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A Theoretical and Experimental Study of 
Shock-Tube Flows 


I. I. GLASS* ann G. N. PATTERSON? 


Institute of Aerophystics, University of Toronto, Canada 


SUMMARY 


During the past 10 vears the shock tube has achieved singular 
recognition as the most important facility for the investigation of 
numerous fundamental problems in nonstationary fluid me 
chanics. The shock tube has been applied successfully to the 
studies of the transition through shock waves and contact frouts; 
wave interactions; flows in the subsonic, transonic, and super 
sonic flow régimes; wave diffraction and refraction; condensation 
effects; sound speed measurements; shock loading of structures; 
relaxation phenomena in gases; flame propagation; chemical 
kinetics; boundary-layer growth and shock-wave attenuation; 
atomic physics of states behind strong shock waves and mag 
neto-hydrodynamics. (See reference 1 for an extensive bibli 
ography.) Thermodynamic states of short duration are avail 
able that have a range of physical quantities from nearly absolute 
zero pressures, temperatures, and densities to high pressures, 
densities, and temperatures, the latter being of the order of tens 
of thousands of degrees Kelvin. This wide utilization of the 
shock tube for basic research work almost tempts one to rename 


it the ‘‘test tube’’ of nonstationary fluid mechanics 

The purpose of the present paper is to give an integrated ac 
count of the shock tube and some of the recent theoretical and 
experimental studies of shock-tube flows 

The reported data for real flows are based mainly on the one 
dimensional, nonstationary research work that was conducted 
at the Institute of Aerophysics during the past 5 years. Included 
are theoretical and experimental investigations of the origin 
problem, total shock-wave attenuation, normal shock reflection, 
simple and double shock wave refraction, head-on collision of 
shock waves, head-on collision of shock and rarefaction waves, 
and the overtaking of two shock waves. In addition, qualita 
tive experiments are presented on normal reflection of rarefaction 
waves, double refraction of rarefaction waves, head-on collision 
of rarefaction waves, head-on collision of shock and wire screens, 
and the emergence of shock and rarefaction waves from the open 
ends of a shock tube into the free atmosphere. 

A comparison is made between theory and experiment, and 
the deviations are discussed. Where a theoretical treatment is 
lacking, the experiments are analyzed qualitatively and the flow 
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quantities, de 


cated 


termined from wave speed measurements, are indi 


NOTATION 


speed of sound 

x/t, a characteristic iu the (x, ¢)-plane 

internal energy per unit mass 

mass flow 

pressure 

time 

particle velocity 

speed of shock wave 

position along shock tube measured from dia- 
phragm 

density 

mean free path 

Viscosity 

kinematic viscosity 

entropy 

coefticient of specific heat at constant pressure 

coeflicient of specific heat at constant volume 

chamber length 

flow Mach Number 

Reynolds Number per foot 


temperature 


forward facing shock wave (particles enter from 
right ) 
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PLATE 1. The wave system produced in a real shock tube from | if | t 
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Vo 0.50 
Fic. 1. Theoretical diagram of the (x, ¢)-plane showing the . 
waves and states in the chamber and channel from the instant a 
the diaphragm ruptures. U = particle path, D = diaphragm, | 2 
H = head of rarefaction wave, and 7 = tail of rarefaction wave | ” 
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SHOCK-TU 


= backward facing rarefaction wave (particles enter 


from left 


( = contact front moving to the right 

He Air = helium and air separated by a diaphragm in a shock 
tube 

He Air = helium and air separated by a contact front 


INTRODUCTION 


Bes SHOCK TUBE was first utilized in France by 
Vieille in 1890, for flame propagation experiments. 
Its application was improved considerably by Payman 
and Shepherd, during the period of 1931 to 1943. 
The most recent and fruitful applications of the shock 
tube were made at universities and government labor- 
Noted among these are Princeton, Michigan, 
Toronto, Cornell, Lehigh, Johns Hopkins, 
NACA and NOL. The Prof. 
Walker Bleakney and his associates at Princeton has 
contributed directly to the development of this field 


of research on this continent. 


atories. 
Brown, 


pioneering work of 


the 
various flow phenomena the 
2- by 7-in. shock tube and the 3- by 3-in. wave inter- 


action tube at the Institute of Aerophysics. 


The photographic records used to illustrate 


were obtained from 


PART (1) THEORETICAL ASPECTS 


1.1) Relevant Wave Theory 


The theory of the shock tube is given in some detail 
in references | and 2. Asa result only the basic con- 
cepts and fundamental equations will be presented here. 
The simple shock tube consists of a straight tube of uni- 
containing a diaphragm which 


form cross section 


separates two compartments (Fig. 1). One compart- 
ment (4) is filled with a gas at a high pressure and is 
designated as the high pressure chamber (or simply 
the chamber), the other (1) contains a gas at a lower 
pressure and is called the low-pressure chamber (or 
Initially, the two gases are in thermal 
equilibrium. When the diaphragm is ruptured in the 
ideal shock tube, a wave system is generated, and its 


the channel). 


representation on the (x, /)-plane is as shown on Fig. 1. 
This wave model may be deduced theoretically from a 
discussion of the transition relations for shock fronts 
and rarefaction waves by utilizing the (p, «)-plane for 
their representation. The postulated wave system 
is also observed in real flows, but it is of a modified 
form. The real flow that appears on the (x, ¢)-plane 
schlieren photograph of Plate | differs in many details 
from the idealized flow shown on Fig. 1. Theoreti- 
cally, it is assumed that the diaphragm is instantane- 
ously and completely removed at time f = 0. Asa re- 
sult a plane shock wave followed by a contact surface 
is instantly transmitted into the channel and a centered 
The 
rarefaction wave is represented as a pencil of straight 
lines (Mach waves or characteristic lines) emanating 


rarefaction wave is propagated into the chamber. 


BE FLOWS vo 


The flow quantities are constant along 
Conse 


from the origin. 
each Mach line but vary from line to line. 
quently state (2), which is compressed by the shock 
wave, and state (3), which is formed from the isentropic 
expansion of the gas in the chamber through the rare- 
faction wave, are uniform regions. The particle ve- 
locity and the pressure in the two uniform states are 
identical. However, because of the different formative 
processes involved, the temperature, density, and en- 
tropy of the two states are different. Essentially, they 
contain dissimilar gases and are therefore separated by 
a contact surface. If the specific heats are assumed to 
remain constant, then the wave speeds and the proper- 
ties of the uniform states are readily determined from 
the initial conditions by matching the unsteady isen- 
tropic rarefaction wave values of the pressure and the 
velocity with those of the steady, irreversible adiabatic 
shock-wave relations. The utilization of these ideal 
conditions is justified from experimental evidence, and 
as a result equations may be derived which determine 
The theoretical basis for the 
It will 


all of the flow properties. 
following relations is discussed in reference 1. 
be noted that only the initial conditions and the shock 
strength are required for the calculation of the physical 
quantities of the flow. 

(a) The Strength of the Shock Wave (P2) as a Function 
of the Initial Conditions. 
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Fic. 2. Transition front profiles through a complete centered 
rarefaction wave showing the various flow parameters. Nott 
The incomplete centered rarefaction wave profiles are identical 
in shape up to the tail of the wave 
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Fic. 4. Representation of the transition through a strong shock 
front 
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Fic. 3(b). Variation of flow properties Re per ft., .W/, u, along L, = shock thickness due to active degrees 
a shock tube. LL» = shock thickness due to inert degrees 
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Variations of the flow parameters in air behind a 


Fic. 5. 
normal reflected shock wave with the incident shock-wave Mach 
Number (W),). 


Eq. (1) indicates that the shock pressure ratio (P»2) is 
a function of the diaphragm pressure ratio (P\,), the 
internal energy ratio across the diaphragm (/j4), and 
the specific heat ratios of the gases in the chamber 
(ys) and in the channel (7;). 

An examination of Eq. (1) shows that there is a limit 
to the shock strength that may be produced even with 
infinite diaphragm pressure ratios, and for the case 
where Ey, < 1, 


[Po Ip, 9A 1+ (Q@ By 44) (2) 


Thus the criteria for strong shock waves produced in a 
simple shock tube are: 
(1) A large diaphragm pressure ratio, Pi; —> 0. 


2) A small energy ratio across the diaphragm, /y4—> 


This immediately suggests the use of such gas com- 
binations as He!|Air, H»|A, or Hs SFs which give a 
small /y4 value. 
channel to further reduce /\, is possible but not prac- 
Subsequently, other methods employed for this 


Heating the chamber and cooling the 


tical. 
purpose will be discussed. 
) The Density Ratios. 

T34 = [PuP ]' ” (3 

Is) = (1 + aP>) (ay a Po) (4) 


} 


FLOWS a 


(c) Speed of Sound and Temperature Ratios. 


Ay = Tu” = Pu™ = [PuPx}* (5) 
Ag = In” = | Pate + = | . (6) 
1+ a Px 

(d) Shock Wave Velocity or Mach Number. 

Wu = [8,01 + aPx)]“* (7) 
(e) Contact Surface or Particle Velocity. 

Un = (Px — 1)/¥1[Bi(a1P2 + 1)] (8) 
Usa = (1/ya8s) [1 — (PrP o1)™] (9) 


(f) Velocity of Head and Tail of Rarefaction Wave. 


Head: Cy = —1 
Tail: Cs, = (1/ysBs) [1 — (PuPa)™] — (10) 
[PisP2 |” 
(g) Local Mach Numbers. 
Mz = (1/Bsys) [(PiusP21) a 1] (11) 
My = (Po — 1) ¥1(81P 21 (a, + P) | (12) 
(h) Local Reynolds Numbers.—If the relation u« 7” 
is utilized, then, for air (k = 0.76), 
T'a4U3 
(Re ‘l)s = = Pa Pi <- (13) 
MsQ4 T3°-* 
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PLATE 2. Normal reflection, refraction, and progression of a 
shock wave through the cold region behind the contact front 
ps 3,068 mm. Hg a 1,136 ft. per sec 
p 20 mm. Hg u 2.34 
P; 153 U2 1.76 
P ae VU 2.7 
VU 1. 4¢ 
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The limiting values of the flow quantities as predicted 
by the above equations are given in Table 1. An ex- 
amination of Table | reveals that the uniform states 
(2) and (3) possess extreme flow conditions. Very high 
temperatures and low Mach Numbers prevail in state 
(2), while very high Mach Numbers and low temper- 
ature are theoretically attainable in state (3). The 
attainability of such states immediately suggests that 
the uniform states may be employed for aerodynamic 
testing. So far, this has not been achieved because of 
the lack of reliable, rapid-response instruments that 
can measure directly the physical quantities of a flow 
of extremely short duration. Furthermore, it has not 
been established whether it is valid to apply nonsta- 
tionary flow test results, which do not possess fully de- 
veloped thermal and velocity boundary layers, to sta- 
tionary flow problems. 

(t) Transition Properties of a Centered Rarefaction 
Wave Moving to the Left (Fig. 1). 


x/t{=u—a (16) 
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u+ [2a/(¥s — 1)] = 2a4/(¥1 — 1 17) 
a/ag = 1 — [(n + 1)/azy] 18 
T/T; = }1 — [(n + 1)/aa]}? 19) 
p/p, = 41 — [(n + 1)/ay]} '”™ 20) 
b/ps = {1 — [(n + 1)/as]}'”? 21) 
M = (n+ 1)/(1 — Byn) 22) 


The profiles are plotted on Fig. 2. It is seen that 
the shape of a rarefaction wave is fixed by the initial 
conditions in the chamber, and its extent depends on 
the diaphragm pressure ratio. When the diaphragm 
pressure ratio is infinite, a complete centered wave re- 
sults, and the extent of the profile is complete—that 
is, the density, temperature, pressure, and speed of 
sound are zero; the particle velocity is 5a,; and the 
Mach Number is infinite at the tail of the wave. For 
incomplete centered rarefaction waves, all the deriv- 
atives are discontinuous at the tail and at the head of 
For a complete wave all the derivatives are 
At the 
tail of the wave they are zero for the density, pressure, 


the wave. 
discontinuous at the head of the wave only. 


and temperature, so that these profiles are tangent to 
the n-axis. The derivatives for the speed of sound and 
particle velocity are finite. The Mach Number deriv- 


ative is infinite. 





th ti 
\ 


= 


(IO) ig) 


croc y 


= 
a Sf 2 é 
VA iz 


5 j 
*R oe oe (IO) \'9 | 
IR 


k 


k 


(6) 








RR 
(4) 
Fic. 8. Interaction of a rarefaction wave and a contact surface. 
U4 t, 
| 
\ | 
\ ‘ 
OA @ a \. 2 
: i 
x . 
rx 
(4) a @ OS 
| | | 
(5) | (I) (5) | (I) 
| x 





The interaction of a rarefaction wave and a stationary 
contact surface 
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A typical set of transition profiles drawn to a scale 
for the shock tube when Py, = 20, case Air Air is shown 
on Fig. 3. The temperature discontinuity at the con- 
tact surface will give rise to a molecular transfer of heat 
between regions (2) and (3). The transition through 
the contact front is analogous to that of a shock front 
and requires the solution of a set of nonlinear equations 
in heat conduction.* 

It should be noted that the greatest entropy change 
exists across the contact surface and arises from the fact 
that state (4) has an entropy that is much lower than 
state (1). 

The matching of the pressure and velocity from the 
steady and unsteady solutions for the shock and rare- 
faction wave, respectively, gives rise to two conserved 
flow parameters on either side of the contact surface 
ie., the energy constant C,7, for states (1) and (2), 
and the momentum constant [2/(y, — 1) ]a, for states 


$) and (3). 


1.2) Deviations in a Shock Tube Due to Variable Specific 

Heats 

The basic equations [Eqs. (1) to (22)] for the flow in 
a shock tube were derived on the assumption that the 
specific heats remained constant. It was shown in 
Table | that for extremely large diaphragm pressure 
ratios (Py,) and small energy ratios (//\4), extremely 
high temperatures could be achieved in state (2) and 
extremely low temperatures in state (3). Under these 
extreme temperature conditions the assumption of con- 
stant specific heat is no longer valid. For the rarefac- 


tion wave the basic equation becomes! 


a4 a ()" p (23 
m4 Ps 


where ® is a correction that is a function of the gas 


u = U4 TT 


properties including the deviation from constant spe 
cific heats. 

The value of ® in general is @ < 1, and for constant 
specific heat @ = 1. Thus the particle velocity (us) is 
increased as a result of the de-excitation of the rota- 
tional degrees of freedom of the molecule due to the 
decreased temperature. In the case of Air Air when 
Py, = 10,000, 7; = 0.16 K 300°K. = 48°K., which is 
below the boiling point of nitrogen and oxygen. Thus 
condensation phenomena can be expected, and constant 
specific heat can no longer be assumed under these con- 
ditions, which are readily obtainable in a shock tube. 

For the case of the shock front, the transition profiles 
for constant specific heat and variable specific heat are 
given graphically in Fig. 4. Immediately at the shock 
front, the temperature jumps almost discontinuously 


toa high value. This is due to the fact that the active 


degrees of freedom (translation and rotation) are ex-. 


cited nearly instantaneously in the distance of a few 
mean free paths (A). Since the previous equations for 
the shock wave were based on constant specific heat, 


the initial temperature jump is accurately predicted 











Fic. 10. Points in the (Y, 7)-plane where the head of the 
reflected rarefaction wave overtakes the tail of the incident 
rarefaction wave (.Y;, 7;), the contact surface (Y,, 7, ind the 
shock wave (X,, r;) 


























¢1,!) 
(4) x 
X Xs= $X, X. 
Fic. 11 Points of maximum flow duration YY; and , for 
regions (3) and (2), respectively Ar; = maximum flow dura- 
tion for state (3 Ar, = maximum flow duration for state (2). 








ta 











Fic. 12. The head-on collision of two unequal shock waves 
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from the Rankine-Hugoniot relation {[Eq. (6)|. After 
the initial jump, thermal equilibrium is reached in an 
exponential manner during the relaxation period when 
energy is transferred from the active to the inert de- 
grees of freedom of the molecule (vibration, dissociation 
electronic excitation, and ionization). This stage is 
marked by a decrease in the temperature which takes 
place in a relatively longer time and over a greater dis- 
tance than that required to excite the active degrees, 
and, therefore, an extended shock wave transition pro- 
file results. If it is assumed that in air the vibrational 
degree is the only inert degree of freedom that is ex- 
cited at temperatures up to about 5,000°K., then it is 
possible to calculate the effects of variable specific heat 
by established methods (reference 4). 

The deviations are as follows. Up to shock wave 
Mach Numbers W}, & 3, the divergence is negligible. 
For Wi > 3 the deviations become monotonically 
greater and the temperature falls, and the density, 
pressure, and particle velocity rise with increasing 
wave Mach Number. For a given diaphragm pressure 
ratio the shock-wave Mach Number is somewhat 
smaller for the variable specific heat case. If the 
equilibrium values are denoted by an asterisk (*), then 
the above results may be expressed symbolically as 


p2* > p2, p.* a pe, Us* > ie, Y Pi < T», and Wr - < Wi 


(1.3) Normal Reflection of Shock and Rarefaction Waves 


In the previous sections the chamber and the channel 
were assumed to be of infinite extent. In actual shock- 
tube work the chamber and the channel are of finite 
lengths with open or closed ends. As a result, after a 
few milliseconds, the flow in a shock tube becomes 
complex because of the ensuing wave interactions. 

Fig. 1 shows a simple shock-tube diagram with both 
ends closed and the resulting normal reflection of the 
rarefaction wave at the end of the chamber and the nor- 
mal reflection of the shock wave at the end of the chan- 
nel. Upon reflection of the shock wave, a new uniform 
state (5) is formed where the particle velocity behind 
the reflected shock is zero, but extremely high tem- 
peratures, densities, and pressures are attained. When 
the rarefaction wave reflects from the chamber wall, it 
also forms a new uniform state (6) of extremely low 
temperatures, densities and pressures, and zero particle 
velocity. It is theoretically possible to achieve a total 
vacuum in state (6) for Air! Air, Ba = 1, when. fy = 
1,620 or Fa = 12:6. 

The following relations apply across a normal re- 
flected shock wave if the specific heats are assumed as 
constant:! 


(a) Pressure Ratio. 
Ps = (ay + 2—-— Py) (1 + a Pp) (24) 


For extremely strong shock waves Pi, > 0 and Ps: > 


Qa) + bi 
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(b) Density Ratio. 
Ts. = (1 + a Poe) /(a1 + Pee) (25) 
(c) Temperature Ratio. 
Ts2 = Ps (a, + Ps2)/(1 + aPs2) (26) 


(d) Wave Speed. 
= = ia — 1) Pre 


WV = 1 (27) 


[(ay + 1) Pyo(ay + Py) 
[2 - (a —_ 1) Py] (Q@y + P12) 


(28) 


Wo Wi = 


For extremely strong shock waves Pj, > 0 and Hx, + 
Wi 2/ a. 

When variable specific heats are considered, Eqs. 
(24) to (28) are no longer applicable. The method that 
was utilized for the solution of incident shock wave 
with variable specific heats is also employed here, with 
the additional boundary condition that the particle 
velocity in state (5) behind the reflected shock wave is 
zero. ‘The final results are shown on Fig. 5 for the nor- 
mal reflection case without and with variable specific 
heats, assuming no dissociation. It will be noted that 
temperatures of many thousands of degrees Kelvin may 
be readily attained in state (5). For example, for an 
initial temperature of 300°K. and a wave Mach Num- 
ber of 5, the variable specific heat temperature behind 
the reflected shock is already 3,000°K. 

This region along with state (2) behind strong shock 
waves offers some promising possibilities in the fields of 
atomic physics and magneto-hydrodynamics. A typi- 
cal (x, ¢)-plane schlieren record of a normal shock re- 
flection appears on Plate 2. 

The following equations give the properties of the 
uniform state (6) behind the normal reflected rarefac- 
tion wave:! 

(a) Pressure Ratio. 


Po = 2Py"* — 1 = 2(PuP2)™ — 1 (29) 
(b) Temperature Ratio. 
Tau = 2P,," — 1 (30) 
(c) Density Ratio. 
Te” = 2P3" — 1 (31) 


In terms of the Mach Number (1/3) in state (3), the 


above becomes 


{(2/(v1 — I] - My 


Ps = * (32) 
‘U2 in — D] + MS 
T., = §(2/(¥4 — 1) — M3)? (33) 
U[2/(¥. — 1)] + MJ 
[2 (¥4 —_ 1)] —_ M; - 
T'64 = (34) 
[2 (¥4 — 1) | + M, 


An (x, ¢)-plane schlieren record of the normal reflec- 
tion of a rarefaction wave is shown on Plate 3. Only 


—- 
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‘ : : , ; : PLATE 4. Normal reflection and refraction of a shock wave at 

\ um- PLATE 3. Normal reflection of a rarefaction wave. Air Air an Air He contact front. 

shind R; = head of incident rarefaction wave a4 1,140 ft. per sec Ss shock wave He) ps 3,080 mm. Hg P 7.2 
R; head of reflected rarefaction wave ts 1,200 mm. Hg ( contact surface (Air) pi 1.9 mm. Hg u 1.8 
Py 12.0 pi 100 mm. Hg a 1,128 ft. per sec 
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(33) | 
Piate 6. The refraction of a shock wave at a layer of helium. 
Case: Air He Air 
S incident shock rs transverse shocks 
(34) S transmitted shock a 1,137 ft. per sec 
: S emergent shock pi,2 747 mm. Hg 
} Plate 5. The refraction of a shock wave at an Air|He S reflected shock Wi 1.82 
flec- | contact front produced by breaking the diaphragm. Ps, = 10, R rarefaction wave u 1.70 
pi = 100 mm. Hg (Air), ps = 1,000 mm. Hg (He), Ps: = 4.60, ( contact surface P 3.70 
nly and = 2.02. ( contact surface a 3.20 
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PLATE 7. The refraction of a rarefaction wave at a layer of 
helium 
R head of rarefaction wave Ds 7.5 mm. Hg 
( contact surface Pis 100 
fi.2, 754 mm. Hg T1,2 ‘ a7 32°C 


the head of the wave shows up in the schlieren photo- 


graph. 
(1.4) Normal Refraction of Shock and Rarefaction Waves 


After normal reflection the shock wave eventually 
refracts at the contact front produced by the ruptured 
diaphragm. The interaction is ‘illustrated on Fig. 6 
and Plates 4 and 5. 

Two cases are possible: a transmitted shock, con- 
tact front and reflected rarefaction wave, or a trans 
mitted shock, contact front and reflected shock wave. 

Mathematically it is simpler to consider this problem 
with a stationary contact surface. It also may be 
set up in this manner experimentally. The interaction 
is shown schematically on Fig. 7. 

The case of the reflected rarefaction wave is given 


by the equation,’ 


Ps” + f( Ps — Pss) _ —g-1=0 (3d) 
a P34 + Piz 
where 
f = [(a5 + Pes) (asPo1 + 1)] (36) 
g = (1 — Pi) [85/(asPs1 + 1)] (37) 
Eys > (a1 + P51)/(as + Pos) (38 


It is also shown that considerable attenuation or ab- 
sorption of the pressure ratio across the refracted shock 
wave results when contact surfaces such as Air) H» or 
Air||He are employed even at room temperatures. 
The case of the reflected shock wave is described by 
the following relation: 
a+ P5sP 43h Px =] P 


E 5 = k Vv Py + ; (39) 
(1 — PsP 43)? Vas t+ Pr 


FEBRUARY, 1955 
where 
h = (asP53 + 1)/(ay + Poss) (40) 
k = (1 — Pss)/V asPss + 1 (41) 
Fis S (ai + Pss)/(a5 + Pee) (42) 


Under these conditions the pressure ratio across the 
refracted shock wave is increased or amplified when 
contact surfaces such as Air|}|A or Air }COs are utilized. 

The problem may be extended to the case of double 
refraction by employing a layer of gas. This inter- 
action is illustrated on Plate 6. A complete analysis 
of normal shock wave refraction may be found in 
reference 5. 

The refraction of a rarefaction wave at a contact 
surface is also produced from the basic shock-tube 
flow as shown on Fig. 8. Here, the rarefaction wave 
generated by the ruptured diaphragm undergoes nor- 
mal reflection at the end of the chamber and refracts 
at the contact surface. Two cases are possible: 
either a transmitted rarefaction wave and a reflected 
rarefaction wave, or a transmitted rarefaction wave and 
a reflected compression wave that steepens into a 
shock wave. Again the problem is simplified by con- 
sidering a stationary contact surface (Fig. 9). 

For a reflected rarefaction wave, 


(G1 Bs) [(PsPs)” — 2(Ps5)” + if 
l — (PoPa)" 


His S (81/85) [1 — Pas™)?/(1 — Pas)?] (44) 

For a reflected shock wave, 
V & (1 — Pos) 
V aP34 + 1 


Pas" 


ees 


| Ds , + 
| fa5 = = 1.) 
% By (P31P45)" — 1] 
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Fic. 13. The head-on collision of shock and rarefaction waves. 
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The overtaking of two similarly facing shock waves 


Fic. 14. 


E., > (B:/Bs)[(1 — Pas™) (1 — Pas” (46) 


IV 


In this case as well the problem may be extended to 
double refraction. The complete analysis of this prob- 
lem will be given in a future UTIA report. 

The double refraction problem for a rarefaction wave 
at a helium layer is illustrated by an actual schlieren 
photograph of the (x, ¢)-plane on Plate 7. The re- 
flected waves did not appear on the records. 

Another problem of interest on this record is the 
stability problem of a contact surface. It is seen that 
a contact surface of the heavy gas (air) accelerated 
(at 6.7 X 10° ft. sec.*) into a light gas remains stable 
(uniformly thick Air) He contact front), whereas for 
the He} Air contact front, a light gas accelerated (at 
2.8 X 10° ft. sec.*) into a heavy gas becomes unstable.® 


(1.5) Optimum Chamber and Channel Lengths for a Shock 

Tube 

The problems of normal reflection and refraction in- 
fluence the choice of optimum chamber and channel 
lengths. Essentially this amounts to a knowledge of 
the points in the flow where the reflected rarefaction 
wave overtakes the tail of the incident 
wave (X3, 73), the contact surface (Y,, 7,), and the 
shock wave (Y,, 7;). The reflected rarefaction wave 
eventually will destroy the flow in states (3) and (2) 
and will also overtake and decay the shock wave. The 
above points are shown on the dimensionless CY, 7)- 
plane in Fig. 10, and they may be determined from the 


rarefaction 


following relations :’ 


XY; = [lag — 1 — ay P34°'|Po sic (47) 
rm = PP, Oo" (48) 

¥. = 2(ay — 1) [1 — Pos] Poys *** (49 
T. = 273 (50) 


U 
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X = WA ur 


ts = 2An73/(Un + An — Wr) (92) 

When the above relations are plotted as functions 
of Py, it is found that shock decay starts at shorter 
distances from the diaphragm for low values of Py. 
For a chamber length of 1 ft., shock decay occurs after 
approximately 4 ft. for He Air and 15 ft. for Air Air. 
It is essential, therefore, to use a chamber of several 
feet for He! Air runs in order to prevent the destruc- 
tion of the flow pattern by the rarefaction wave. 

An examination of Fig. 11 shows that the maximum 
flow duration Az; for region (3), occurs at .Y; for all 
diaphragm pressure ratios (P4,) for which ./; > 1. If 
Ms < The maxi- 


mum flow duration Ar. for region (2) occurs at X,. 


1, the maximum occurs at Y = 0. 


Thus for maximum aerodynamic testing the model 
In order that the re- 


should be placed at Y; and YX, 
flected shock wave does not reduce the value of Aro, 


the channel should be of a length .Y;, such that the 


reflected shock meets the head of the rarefaction wave 
at Y,.. The length -Y; is determined from the following 


equations: 


7 = (X; WA 14) — (X; _ Aa WA 44 (53) 


X54 = [r-A 14 . 2 &.¢ W2) (1 Wa) . te (1 W's) (54) 


If region (3) is under investigation, then it is suffi- 


cient to use a channel length Ys, such that the reflected 
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PLATE 8. Collision of two unequal shock waves in air 
1.5mm. Hg P 11.18 P 5. Bt l 0.7 
1,130 ft. per sec Hi 3.12 WW 2.17 P ( 0 
P ; 10 
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| shock wave strikes the contact surface when the head i 
tf ih ie ay of the reflected rarefaction wave overtakes the tail of — | 3 
R (3) the incident rarefaction wave (X3, 73). This practical 
simplification eliminates the necessity of solving the 
above shock refraction relations. Since 7, = 273, there- 
fore, 
REGION Xs = (1/2)X; (55) | 5 
OF j 
PENETRATION 1.6) Some One-Dimensional Wave Interactions f 
” (2) In addition to the previous interactions directly \ 
connected with the flows in a simple shock tube, the 
following interactions have been investigated by mem-  } 
Rr “R bers of the Institute of Aerophysics. 
(1) \ (a) The Head-On Collision of Two Shock Waves.° 
After the interaction, two receding shock waves result, 
\\ x separated by a contact surface, as illustrated on Fig. 
toa gu dave _——— 12. It is desired to know the strength of the shock 
waves, Py, and P39, after collision, as a function of the j 


initial conditions in state (1) and the two incident 


Fic. 15. The head-on collision of two unequal rarefaction waves. 
shocks P;; and P»,. States (3) and (4), as well as the 




































































contact surface velocity, may be solved in the usual 
manner if the shock strengths are known. Across the 
contact front ps3 = p;and uv; = uy. The solution to this 
problem is given by . 
t 
4 
; Py = V1 + aPssX 
! ) . PrP — Pris sa ot 
! \ 4 — ¢ B+ ( + 1 (90) 
; Vv P3;3(aP yp »P 45 + Ps) 
| (4) 
1 — Pr Py», — 1 
i= p = } 
\ + o,Pis V1l+ aPy a 
= — (4) 
C 5 
k | 1 a (a,P i» 1) — Pis(aa ar P\5) 
(4) | (7) | 6) (5) Pye (ay, Py) (a,P 15 + 1) 
! 
— Since Px, = P4-Pi2-Ps, then for the case of equal 
Fic. 16. The — — = — the shock shocks where P. = P,,. the two roots of Eq. (56) yield 
wave 1rougn < usemann nozzle. e ° . 
= the following solutions: 
For the positive root, 
a P35 = (2 + a= P,) (1 + a P45) (5S) 
t Af | aa ae ae , . , : 
| } “6, | | | | lhis is identical to Eq. (24) and shows that for the 
: | | j ; collision of equal shock waves the problem is equiva- | 
| / | | | lent to normal reflection from a solid boundary, and 
\SR cn WR, 7 no contact surface should be observed. Eq. (56) 
ta (5) [7 | . (5) | (3) F | may be employed to show that for the collision of a 
te e 
XQ / / () shock and a sound wave the strength of the sound wave 
\ j nee” 
(4) WN (as (4) S is increased. 
(1) é P . : 
()_ |x z For the negative root Py; = P35, a linear result, which 
CASE M;<! is applicable only to sonic reflection. 
1 Actual (x, ¢)-schlieren photographs of the head-on 
collision of two unequal and equal shock waves are 
Ay (4) (I) shown on Plates 8 and 9, respectively. ; 
‘ 
| (b) The Head-On Collision of Shock and Rarefaction : ai 





Fic. 17. Flow in a shock tube with a large chamber. Waves.*—In the interaction illustrated in Fig. 13, the 
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Head-on collision of shock and rarefaction waves 
- transverse shock waves, P;; = 1.96, P32. = 2.10, 
746 mm. Hg, and a, = 1,135 ft. per sec 


in air of strength 0.48. 
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The overtaking of two similarly facing shock waves 
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initial conditions are given by state (1) and the strengths 
of the colliding shock and rarefaction waves. As a 
result of the interaction a receding shock and rarefac- 
tion wave are generated. The two waves are separated 
by a contact region of variable entropy that is due to 
the curving shock wave as it passes through the rare- 
The shock-wave pressure ratio is ampli- 
Across the resulting contact 
The solution is given 


faction wave. 
fied after the interaction. 
front p; = ps and uw; = 4. 


by 


Pes = | + ¥Y |] + a P39 xX 


| D —E [(Ps2P15)' TF es (Py»)' ae t 


(59) 


where 


(c) The Overtaking of Two Similarly Facing Shock Waves.' 
State (1) is given along with the strengths of the two overtaking shock waves. 
a transmitted shock and a reflected rarefaction wave separated by a contact surface, for gases of y < 5/3. 


fictitious gas in which y > 5/3, a reflected shock wave is also possible. 


B l | 5) 
il -- ans (Ps — Py) 
(a, + P,;) (aP 15 + 1)P15 


where 


Ps = Py: Px 


A detailed discussion and a graphical analysis of this 
equation is given in reference 10. An actual (x, /)- 
plane schlieren photograph of this interaction is shown 
on Plate 11. 

(d) The Iead-On Collision of Two Rarefaction \Waves."' 
—The interaction appears on Fig. 15, and the problem 
has been solved in a closed form for a monatomic and 
diatomic gas even in the region of penetration. There, 
the characteristic lines are no longer straight, and a 
general solution is obtained by using Riemann’s inte- 


gration method. The following results may be ob- 


tained (see Fig. 15): 
Us = (1 7101) (P,,? - P.,* ) (63) 
Ps," = Py” + Px? — | (64) 


If the waves are of equal strength, then Py, = P2 and 
Us = 0. Similarly, Ps: = 2Py" — 1 and is equiva- 
lent to Eq. (29) for the normal reflection of a rarefaction 
wave from a solid boundary. For this case it is shown 
in reference 11 that the last characteristic line to reach 
x = 0 is x = at for a monatomic gas and x = 2a 
for a diatomic gas. Under those conditions Py,” = 
1/2 and P3, = 0. Theoretically it is possible to 
achieve a complete vacuum in state (3) only after an 


infinitely long time. 
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(Pis + a1) P45 
Ps = a Px, 


SCIENCES FEBRUARY, 1955 
(Ps — 1) } 
Vv ] + P35 


Vor +1 [lL — Py +] 


Da Pie" 


(60) 


= \ aT l N P., (ay + P51) (] + a P3)) | 


It should be noted that Eq. (59) gives a means of 
calculating the final uniform states and wave strengths 
after the collision, but it is not sufficient in itself to plot 
the path of the curving shock through the rarefaction 
The procedure for tracing the shock path is 
Eq. (59) may be employed to 


wave. 
outlined in reference 9. 
study the collision of a sound wave and a rarefaction 
wave. An actual (x, /)-plane schlieren photograph of 
a typical interaction is shown on Plate 10. 


*-The above interaction is illustrated on Fig. 14. 
After the interaction there results 
For a 


The solution for all real gases is given by 


a Pi; 7 l 
~ yo gi ~~ Pets 
l si a P,; 
P3"' -—1=0 (61) 
(62) 
The above condition, Py," = 1/2, is achieved in a 


shock tube for the case Air Air when the diaphragm 
pressure ratio Py A 1,620, or at an incident shock 


strength P., A 12.6. (See Fig. 1 for notation.) 


An (x, f)-plane schlieren photograph of the head-on 


collision of two rarefaction waves appears on Plate 


LZ. 


(e) Expansion of the Uniform States Separated by 
the Contact Front in a Busemann Nozzle for High Mach 
Number Flows.'*—-By expanding the flow behind the 
shock wave in a Busemann nozzle, as shown on Fig. 16, 
it is possible to attain high Mach Number flows in a 
shock tube which simulate flight conditions. It is also 
possible to expand the flow behind the contact front 
in a similar manner. <A theoretical and experimental 
study of this problem is given in detail in reference 12. 


It is interesting to note that the wave configuration 
after expansion is similar to the head-on collision of 
two shock waves. That is, two shock waves facing in 
opposite directions are formed, and they are separated 
In practice these waves are not 


by a contact front. 
However, for a theoreti- 


plane, as shown on Plate 13. 
cal analysis this is a reasonable assumption, and the 
good agreement obtained between theory and experi- 
ment justifies the approach. 
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PLATE 13. The expansion of the uniform region behind the shock : PLATE 2 Condensation elects in the rare ms tion wave 
ora wave through a Busemann nozzle (see Fig. 16). Pu = 10, Pox = 2.8, ps = 746 mm. Hg, 73,1 = 298°K., 7; = 
y 210°K., Wn = 1.61, WM; = 1, chamber length = 4 ft., and 
- channel length = 1 ft. SS. = condensation shock 
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19. The flow in a step-type shock tube 


This procedure combines one-dimensional nonsta- 
tionary and two-dimensional stationary flows advan- 
tageously in order to produce a practical utilization of 
the shock tube as a type of high Mach Number flow 
intermittent wind tunnel. 


(1.7) Production of Strong Shock Waves 


The shock tube has recently been utilized for the 
study of the high-temperature states produced by 
strong incident and reflected shock waves. These 
states become ionized and lend themselves to investi- 


gations on magneto-hydrodynamics, relaxation phe- 
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(b) STRONG SHOCK WAVES 


Fic. 20. Total shock-wave attenuation in a shock tube. 
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nomena, and other problems associated with atomic 


physics. 

Strong shock waves may be produced in a practical 
manner by employing the following methods: 

(a) The chamber is filled with a combustible mixture 
(say a stoichiometric mixture of oxygen and hydrogen 
diluted with 80 per cent helium) at high pressure, 
(400- 1,000 Ibs. per sq.in.) and is then ignited. Or the 
high pressure itself is raised until the metal diaphragm 
The latter 
r-ethod appears preferable from the work conducted at 


Lreaks and results in spontaneous ignition. 


Cornell 
manner high diaphragm pressure ratios and extremely 


Aeronautical Laboratories, Inc.'* In this 


low energy ratios are produced as required by Eq. 
(1), in order to produce strong shocks. 

(b) The gas in the chamber can be heated suddenly 
to a sizable temperature by passing a high voltage dis- 
charge directly through it. This process produces a 
high temperature and pressure in the chamber and 
yields extremely strong shock waves. 

) Two reasonably strong shocks, which may be 
generated without undue difficulties, are made to over- 
take, and after the interaction produce an extremely 
strong transmitted shock wave {see Eq. (61)]. For 
example, from Fig. | of reference 10, it can be seen that, 
when two equal shocks of a pressure ratio of 100 over- 
take, they produce a shock wave of a pressure ratio of 
1,200, which is about half the product of the pressure 
of the two incident shocks (this ratio approaches unity 
for weak shock waves) and represents an enormous 
increase in the transmitted shock strength. 

(d) A worth-while improvement in generating strong 
shock waves is realized by using a variable geometry 
shock tube. Such a tube is shown on Fig. 17. It con- 
sists essentially of a large area chamber and a small 
area channel. The wave system produced when 1/; > 
1 and AJ; < 1 is also sketched. The theory for this 
type of tube was first analyzed in reference 14 and is 
developed in reference 1. The greatest advantage 
A,;—> ©, 


occurs when A, A comparison of the shock 
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VELOCITY PROFILES \\ 


Boundarv-layer profiles with increasing time for the 
equivalent Rayleigh problem in a circular tube. 
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strength produced for a given Py for Ay/A,; > @, 
and A,/A; — 1 (simple shock tube), is given on Fig. 18. 
A further increase in shock strength is obtained by 
combining a large chamber shock tube with several 
diaphragms in steps. A simple example appears on 
Fig. 19. It is assumed that, when A,/A;, is large, 
then the shock produced by breaking the first dia- 
phragm, undergoes normal reflection at the second 
diaphragm, which is ruptured because of the sudden 
increase in pressure and temperature. This increase 
results in a Wy, = 3.9, whereas without the first dia- 
phragm Wy, = 3.3. 

Studies conducted on strong shocks” have already 
shown that a new field has been opened for the investi- 
gation of many problems in magneto-hydrodynamics 


and atomic physics. 
ParRT (2) EXPERIMENTAL ASPECTS 


(2.1) Some Design and Construction Criteria for Shock 

Tubes 

The literature on shock tubes has become so ex- 
tensive in the past 5 years that it would be a rather 
lengthy task to review all the applications of the shock 
tube and the resulting design and construction criteria. 
Instead, some general remarks will be made which are 
applicable to all shock tubes. A particular applica- 
tion, the 3- by 3-in. wave interaction tube, is discussed 
in some detail in reference 1, where the pertinent refer- 
ences are also listed. 

The type of shock tube to be used is mainly deter- 
mined by the funds available and the type of problem 
to be investigated. If shock waves, rarefaction waves, 
and contact surfaces are studied, then a simple cylin- 
drical metal tube plus an optical system is quite ade- 
quate. If magneto-hydrodynamic effects produced by 
strong shock waves are investigated, a cylindrical 
metal chamber and a glass channel is sufficient. IJf 
flow studies over models are conducted in the uni- 
form regions separated by the contact surface, it is 
necessary to use a square or rectangular metal tube 
that will accommodate observation windows for the 
optical systems and will readily permit the installation 
hot-wire 


ot aerodynamic models, piezo-gages, and 


anemometers. The magnitude of the cross-sectional 
area will depend on the size of model to be tested or the 


Reynolds Number to be attained in the shock tube. 


BE FLOWS s8Y 


The wall thickness of the tube will be determined by 
the maximum pressure in the chamber. Welded struc- 
tures should be avoided because of the subsequent 
distortions. A combination of high pressure and large 
cross section will result in a heavy shock tube and may 
require hydraulic or other means for moving the vari- 
ous sections. It is convenient to have a chamber of 
2 or 3 ft. and a channel composed of several sections 
of 2 to 8 ft. 
changeable and easily aligned by means of dowel pins 


The channel sections should be inter- 


so that they may be used in combination as required 
by a particular experiment. 

The viewing windows can be installed in one or two 
sections. These may be coupled with the remaining 
channel sections to give an observation window located 
at a desired distance from the diaphragm and will thus 
give a predetermined flow duration. The windows 
may be of standard plate glass for shadowgraph and 
schlieren observations. It is important that the 
windows be selected and installed in a manner that will 
give uniform refraction due to the slightly wedged sur 
faces that are always present in this type of glass. If 
the wedge surfaces in a window are not uniform and 
monotonic, the window should be replaced. For inter- 
ferometric work the windows should be plane and 
parallel optical flats. The windows may be installed 
with Wood's metal, rubber gaskets, or recently de- 
veloped plastic glass-to-metal bonding cement (for ex- 
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diaphragm pressure ratio (/;,) and distance (x) 


ample, Hysol cement) to give a flush and vacuum- 
tight installation. 

The most satisfactory diaphragm materials are red- 
zip cellophane and cellulose acetate for low-pressure 
work and copper and aluminum sheets for high-pres- 
sure work. The red-zip cellophane has the best ruptur- 
ing properties and leaves the smallest obstruction at 
the diaphragm cross section. Unfortunately, this ma- 
terial is not available in a range of thicknesses and must 
be used in several layers for higher pressure differences 
across the diaphragm. This gives rise to a small vari- 
ation in the time it takes for all the diaphragms to rup- 
ture and results in an overlapping of several wave 
systems at the origin. (See Plate 14.) 

The pressure differential at which a diaphragm will 
rupture depends on the material, cross-sectional shape, 
and area, as well as on the gasket material used to hold 
the diaphragm. If metal diaphragms are used, gaskets 
are not necessary, but the sharp edges of the channel 
section should be rounded to avoid cutting. If cello- 
phane is used, rubber gaskets are bonded to the tube 
sections with (Pliobond or Armstrong) cement. The 
gaskets should have enlarged cross sections so that 
when the tube sections are fastened together, the gasket 
does not squeeze into the tube and cause a discontinuity 
at the origin. 

Shadowgraph, schlieren, and interferometric meth- 
ods have now become standard equipment for shock- 
tube investigations. The twin mirror, single pass, 
horizontal or vertical symmetrically offset schlieren 
system is most useful. First-surface paraboloidal 
mirrors corrected to a quarter wave length of aperture 
F6 to F10 are necessary for good schlieren resolution. 
The types of light sources, continuous or spark, are de- 
scribed in the references. 

Initial conditions in a shock tube are obtained by 
measuring the pressure in the chamber by using mer- 
cury manometers or (Heise type) bourdon gages, and 
in the channel by using mercury or butyl phthalate 
manometers, McLeod, ionization and Pirani gages or 
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(Wallace and Tiernan) absolute pressure gages. The 
temperature of the gases in the chamber and channel 
When gas at high 


pressure is admitted into the chamber, it should be 


is generally room temperature. 


allowed to reach thermal equilibrium (in about a min- 
ute) before the diaphragm is ruptured. Commercial 


evlinders of compressed gas (air, A, He, COs, Ho, No, 


or OQ») are easily available, and they have the advantage 


of a low dew point (about minus 40°C. or lower) which 
eliminates extraneous effects due to water condensa- 
tion. 

Final conditions in the shock tube are obtained by 
measuring the shock wave velocity (IV) continu- 
ously! or by measuring the time interval for the shock 
wave to travel over a given base (accurately measured 
distance between two points). By utilizing the theo- 
retical relations of Part (1) the flow parameters can be 
profiles can be 


calculated. In addition, 


measured directly by using piezo-electric transducers, 


pressure 


temperature profiles may be recorded with hot-wire 
anemometers, density contours are determined with 
an interferometer, and flow Mach Numbers are found 
from shock angles over conical or wedge models. 
When strong shock waves (IV; & 10) are studied, the 
wave front itself becomes luminescent, and its velocity 
may be self-recorded on a drum camera without the 
use of a shadowgraph or schlieren system. 

To date, the wave speed results, along with density 
and Mach Number measurements, have been the chief 
source for verifying the shock-tube theory. Piezo- 
gage data and anemometry data have been confined, 
so far, to weak shock waves. Light-scattering tech- 
niques have also been employed to give a measure of 
the degree of liquefaction in the cold region (3) behind 


the contact surface. 


2.2) Experimental Procedure 


The following data are based mainly on the wave- 
speed measurements of the (x, /)-plane schlieren records 
obtained with a rotating drum camera in conjunction 
with the 3- by 3-in. wave interaction tube. The drum 
rotates at about 4,200 r.p.m. and carries a strip of 
70-mm. Linagraph Pan film externally. The schlieren 
system is equipped with 12-in. parabolic mirrors. The 
arrangement makes it possible to obtain continuous 
(x, t)-plane records over 12 in. and about 14 millisec. 

The photographic records show the propagation of 
discontinuities in the density derivative and, as a re- 
sult, reproduce the paths of the basic flow elements 
found in nonstationary flows—that is, shock waves, 
contact surfaces, rarefaction waves, vortices, eddies, 
and boundary-layer effects. 

The film speed is determined accurately to about 
0.07 per cent. Wave angles are measured to better 
than 0.5 min. by means of a Hilger projector. As a 
result, wave velocities can be found to better than 2 
per cent, even at diaphragm pressure ratios of 10,000. 
To produce diaphragm pressure ratios of this order it 
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is necessary to employ channel pressures of 0.4 mm. 
Hg and experimental errors of about 9 per cent are 
possible. By taking several runs at high pressure 
ratios and utilizing average results, it is possible to 
minimize this effect. 


2.3) The Wave System at the Origin in a Real Shock Tube 


Plate | is typical of the schlieren photographs taken 
at the origin which are presented and analyzed in 
dia- 


Under the cellophane 


When it ruptures, 


reference 1. 
phragm becomes a curved surface. 
three-dimensional effects are introduced in the form 


pressure 


of shock-wave reflection arising from shock curvature. 
Shock diffraction and vortices are also produced by 
the jagged remains of the diaphragm. As a result the 
assumption of one-dimensional inviscid flow at the 
origin is not realized. 

All the records are typical, in that the shock wave 
at the origin appears to be formed within a network of 
converging characteristic lines representative of a com- 
pression wave. However, closer examination of the 
optical records reveals that the method of character- 
istics cannot be applied in order to predict the birth 
point of the shock wave' because the apparent com- 
pression wave is in fact composed of a number of 
These overtake, accelerate, 
The 
wave is initially curved, but the resulting regular and 


coalescing shock fronts. 
and combine to form the primary shock wave. 


Mach reflections form the mechanism that produces 
a plane primary shock wave. The reflections also give 
rise to a train of transverse waves that follow the pri- 
mary shock wave down the tube (even as far as 50 


tube cross-sectional widths). They also propagate 
upstream of the diaphragm. Eventually the trans- 
verse waves become sound waves. However, they 


act as a dissipative mechanism and contribute to the 
total shock-wave attenuation and to the fluctuation of 
physical quantities in the two uniform regions sepa- 
rated by the contact front. 

It is shown in reference | that even for a weak shock 
wave (Ps, < 3) accelerations of the order of 7 X 10° 
ft. sec.? may be produced during the first 10 wu sec., 
and in this respect the rupturing process approximates 
the ideal case of infinite acceleration. 

Plate 1 also shows that the contact front is in fact 
a region which in a matter of 18 in. has grown to an 
apparent thickness of 3 in. 

The region behind the contact front and the rare- 
faction wave is filled with striations that are indicative 
of turbulent and eddying flow. The rough flow is 
caused by the passage of the gas from the chamber 
over the jagged remains of the diaphragm. 

Since the ideal case of an instantaneously plane shock 
Wave is not attained at the origin, the rarefaction wave 
is not of the centered type. The head of the wave 
appears prominently on Plate 1 but not the tail of the 
wave. 
cannot be attributed to the fact that the rarefaction 


The nonappearance of this characteristic line 
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2: Variation of boundary-layer displacement thickness 
6*) with VW. and distance (x) from the diaphragm 


wave is not centered, since this does not affect the 
order of magnitude of the value of the discontinuity in 
the density derivative which responds to a schlieren 
system. 

Summing up then, it may be stated that the flow in a 
real shock tube is a viscous process and deviations from 
an inviscid flow solution may be expected. 


2.4) Some Experimental and Theoretical Aspects of Shock- 

Wave Attenuation 

One of the first quantities to be determined by shock- 
tube investigators was the shock-wave velocity, since 
it was relative easy to measure. It is shown in Part 
(1) that a knowledge of the initial conditions and the 
shock-wave velocity makes it possible to calculate all 
the physical quantities of the flow in a shock tube. 
This approach was followed by many experimentalists, 
and they found that the agreement between theory 
and experiment was good for weak waves but became 
progressively worse for strong shock waves. 

This deviation has been referred to somewhat loosely 
as shock-wave “‘attenuation.’’ As a result of the pres- 
ent set of experiments it is possible to define the term 
‘attenuation’ in a more rigorous manner. It has been 
found that associated with each weak or strong shock 
wave there is a characteristic length (v,;) over which 
the shock forms, accelerates, and achieves a maximum 
velocity.!. Based on the present results a weak shock 
is defined as one for which P2; < 3. The formation dis- 
tance (x,) was found to be somewhat longer for strong 
shocks. 

For weak shock waves the achieved maximum ve- 
locity (IV) is the one predicted theoretically for a 
given diaphragm pressure ratio (Ps Confirmation of 
this point appears on Fig. 25, and it may be noted 
that the wave speed did not change measurably over 
the 142-in. length investigated during the present ex- 
periments. 

For strong shock waves (P2; > 3), it is noticed that 
the shock wave achieves a maximum speed at the end 
of its formation distance (x,), which is lower than the 


theoretical wave speed corresponding to the given 
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diaphragm ratio. Beyond the formation 


length a deceleration occurs, and the shock velocity 


pressure 


decreases monotonically with distance as sketched in 
Fig. 20b. It should be noted that during the present 
experiments, strong shock waves were produced by 
evacuating the channel to pressures as low as 0.4 mm. 
Hg; a procedure which resulted in correspondingly 
high kinematic viscosities (v2) behind the shock wave. 

Recapitulating, it may be concluded that for strong 
shock waves the overall defect or total attenuation in 
shock speed consists of two portions: (a) a decrement 
due to formation and (b) a further attenuation due to 
the distance traversed by the shock wave. These will 
be referred to as a formation decrement and distance 
attenuation, respectively. 

It is reasonable to assume in the case of strong 
shock waves (especially where the kinematic viscosity is 
high) that the formation decrement is due to viscous 
action produced during the shock formation process 
and the distance attenuation results from a boundary- 
layer build-up. The possible growth of the boundary 
layer in a shock tube is illustrated on Fig. 21. It is 
seen that at the position of the shock wave the boun- 
dary layer has zero thickness, since the gas particles in 
the tube have not been affected yet. At the contact 
front the boundary-layer thickness (6) is a maximum, 
since there the particles have been in motion with re- 
spect to the wall for the total length of time (¢;) under 
consideration. Just behind the contact front the 
Reynolds Number increases discontinuously, and for 
an ideal flow without turbulence in region (3) it may 
safely assumed that the boundary-layer thickness will 
decrease in state (3) and through the rarefaction wave 
until it is zero again at the head of the wave. 

Thus the boundary layer in a shock tube at a given 
time (¢t) is a function of position (x). The thickening 
of the boundary layer with time will diminish the in- 
duced mass flow behind the shock wave, and as a re- 
sult its speed will be reduced to balance the lower mass 
flow. No definite mechanism for reducing the shock 
speed is postulated. However, since shock decay is 
caused by a rarefaction wave, it is generally considered 
that the contact front slows up. The deceleration 
produces a rarefaction wave analogous to the classical 
one-dimensional motion of a piston in a cylinder.’ It 
will be shown subsequently that the very opposite 
happens in a real flow—that is, the contact front ac- 
tually speeds up. For this reason, the above concept 
may appear rather artificial. Nevertheless, even if 
the mechanism for slowing up the shock wave is not 
known, it is reasonable to state that the boundary- 
layer growth reduces the mass flow, which in turn de- 
creases the shock speed and is the cause of the distance 
attenuation. 

The present problem may be compared with Ray- 
leigh’s analysis of an infinite flat plate that has been 
started instantaneously from rest to a uniform velocity 
u; at ¢ = 0.'° In the present analysis the equivalent 
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problem is considered—that is, the plate is at rest and 
a flow is generated instantaneously such that its speed 
is zero at¢ = 0,u = u,fort > 0. This is illustrated for 
a cylindrical shock tube in Fig. 22 and was solved in 
reference 17. It is seen that at small times the bound- 
ary layer is thin and uniform over the infinite tube. 
For large times the boundary layer approaches the cen- 
ter of the tube, and the motion becomes a pipe flow. 
In a real shock tube the boundary layer does not grow 
uniformly, and it occupies only a finite length of the 
tube. Furthermore, the end boundary conditions are 
complex because of the presence of the shock wave and 
the contact region, which extends right to the dia- 
phragm station due to the no-slip condition at the wall. 
However, it will be shown subsequently that some useful 
results may be obtained by considering the two prob- 
lems as equivalent. 

A dimensional analysis of the instantaneously started 
smooth tube (otherwise a roughness parameter must be 
included) indicates that the boundary-layer thickness 
(6) or displacement thickness (6*) (see Fig. 21), neglect- 
ing heat transfer, is a function of the following quanti- 
ties (for a real flow the tube walls are cold and heat 
transfer exists between the gas and the walls, and the 
Prandtl Number is an important parameter) : 


6 = d(p, p, uw, ¥, u, d, t) (1) 


(2) 


F(6, p, p, u, ¥, u, ad, t) = 0 


Therefore, there are five dimensionless z-functions, 
such that Eq. (2) may be expressed in terms of the fol- 
lowing dimensionless parameters as:'* 
6 u> = pud_s ut 
d yp/p pw @& 
6/d = g[M, Re, (vt/d*), y] (4) 


Therefore, if tube roughness and heat transfer are 
neglected, then for a given gas (vy) the distance attenu- 
ation that depends on the boundary layer may be ex- 
pected to vary with the Mach Number (J/), the Rey- 
nolds Number (Re), and the time parameter (yf/d*). 

In the case of the shock tube the gas particles at the 
contact front are in motion with respect to the wall for 
the greatest length of time, and it is the time that ap- 
parently would correspond to the time considered in the 
equivalent Rayleigh problem. If the total shock wave 
attenuation is neglected, this time (f,;) is computed as 


(Fig. 21), 


since they are all identical. It is convenient to meas- 


ure /; in terms of the shock position from the diaphragm 
(t; = x,/w,), and the time parameter may be expressed 
as 


vot /d? = wx/Aw (6) 


where A is the cross-sectional area proportional to d’. 
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The above type of analysis cannot be applied directly 
to the formation decrement, since it would also depend 
on the diaphragm rupturing process, which is a complex 
three-dimensional viscous problem. 

In order to check some of the factors contributing to 
the formation decrement and the distance attenuation, 
the following experiments were conducted: To investi- 
gate the formation decrement, different types and com- 
binations of cellophane diaphragms were employed. 
The diaphragms were maintained at a constant pres 
sure difference Ab = (ps — pi), of 100, 500, and 1,000 
mm. Hg over a wide range of diaphragm pressure ratios 
(P;;) in order to study influence of static loading. 

The chamber pressure (ps) was limited to about 7 
atm. It was not possible, therefore, to maintain a con- 
stant channel pressure (p;) over a large Py, range. In 
stead /; was varied to obtain a wide range of P4. 

The influence of extraneous disturbances on the for- 
mation decrement by the shredded dia- 
phragm or the diaphragm breaker plunger were also 
The broken diaphragm introduces a discon- 
This effect was emphasized by 


produced 


observe. 
tinuity at the origin. 
placing protruding metal frames at the diaphragm 
station and noting the differences produced in the 
formation decrement. By placing the diaphragm 
breakers upstream or downstream of the diaphragm 
under identical conditions, similar effects could be in- 
vestigated. 

Finally, to study the distance attenuation, the shock 
velocities were measured with a wave speed camera and 
schlieren system continuously over a 12-in. length at 
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The experimental results may be summarized as fol- 
lows: Cellophane diaphragms, CIL MSTL 600, used 
singly or in combination, gave the closest approxima- 
tion to an ideal rupture when they were loaded to about 
half or greater of their static rupturing load. The 
energy contributions from the breaking cellophane 
diaphragms may be neglected. Introducing discon 
tinuities at the origin or small protrusions has no meas 
urable effect on the shock-wave velocity. The forma 
tion distance where the shock becomes plane and 
achieves a straight path on the (x, /)-plane varies. For 
weak waves (P2; < 3) about 4-tube widths are required. 
For strong shocks where P2; & 17, the formation dis 
tance extends over about S-tube widths. 

The effects of distance attenuation and the formation 
decrement for 1 < Ps, < 10,000 and 33 in. < x < 142 


in. are summarized on Fig. 23. 


The theoretical relations [Eqs. (1) and (7), Part (1 
may be approximated by the equation, 
log Pa = 1.34 (Wi; — | - 0.176 (7 


Whereas, all the experimental results are 


described quite well by the expression 


(Wy 


present 


log Py = 1.33x°-% 1) + 0.08 (S 


It was desirable to have a relation that would be 
applicable to all shock tubes. As a possible solution 
the viscous incompressible equations of motion were 
solved for a square tube having the following initial and 


boundary conditions: 


Ou/Ot = v|(O-u/OZ-) + (O°u/OV 


distance intervals starting at the origin and beyond, up (O<2<z 
to 142 in. (O- to 50-tube width) from the diaphragm u(z, y, O) u; lor W<cve¥ 
° a e ° e a -0 
station. The range of diaphragm pressure ratio (P4;) (9 
was simultaneously varied over the range 1 < Py < u(O, v, ¢) | 
10,000. (On the graphs that follow x is the distance u(z, O, ¢) 0 fort > 0 
; : ; . or | 
from the diaphragm station to the center of the 12- u(Zo, ¥, t) 
in. schlieren beam.) u(Z, Vo, f 
The solutions for the velocity distribution and the mass flow are given by the following: 
ut = > | h 2 ° ans ° b TV 
ll Lo - (1 = cos ae) (1 — cos barje - sin - sin (10 
Uu 1=1b=1 abr? Zi Vv 
m $ ~~ wo (1 — Cos am)? (1 — cos br)? A) (a2 +b 
= . Ya iw aie += 
m; wia=1b=1 ab? 


If the shock attenuation is due to a defect in the mass 


flow, then this may be expressed for air as 


ml (Px a ') 6P + 1 (P2)i + ") 
= (12 
) 7 } (P2;) — | 


m \P>o, + 6 6(P.2); + 1 
By assuming that both problems were equivalent [{i.e., 
Eqs. (11) and (12)], it was possible to plot Fig. 24. 
The results indicate that the viscous flow analysis gives 
a good description of the real flow up to diaphragm 


Beyond 200 the attenuation 
This 


confirms the previous statement that, when the time 


pressure ratios of 200. 
predicted by the present analysis is too large. 


parameter (v//A) is large, the viscous effects may be ex 
pected to be too great, since the flow becomes a non 
stationary pipe flow analogous to the Poiseuille flow for 
the stationary case 

A more realistic approach is to assume that the flow 
consists of an inviscid core and a viscous boundary layer. 








94 JOURNAL OF THE 


The boundary-layer thickness varies with x at a given 
t, between the shock and the contact front. But it is 
assumed that the boundary-layer displacement thick- 
ness (6*), at the contact front, controls the defect in 
mass flow which in turn causes the shock-wave attenu- 
ation. The mass flow ratio depends on the initial and 
final strengths of the shock wave and is found from Eq. 
(12). 
The mass defect at the contact front, 
Am = p2tt.6*D 


perimeter of the shock tube and p», 2, and 
For a square 


where D = 
6* are measured at the contact front. 


tube then, 


56* = (yo/4) [1 — (m/m;j) (13) 
The above definition of 6* may be questioned. How- 


ever, it can be tested experimentally, perhaps by inter- 
ferometric means. Some work in this direction is pres- 
ently under way at the Institute of Aerophysics. 

By utilizing the results of Fig. 24 in conjunction with 
Eqs. (12) and (13), it is possible to plot a family of 
curves of 6* vs. x for varying shock strengths or flow 
Mach Numbers as shown on Fig. 25. During the pres- 
ent experiments the ratio v2./Aw, was held constant 
for a given J\/, and x. The results of the above figure 
indicate that 


dé*/dx @ 0.0004 


over the Mach Number range 1.3 < J/. < 1.5. The 
equivalent value for stationary flow is about tenfold. 
When the lines are extended toward the origin, it may 
be seen that 6* is composed of a sizable formation dec- 
rement that increases with J/. or v» and a distance 
attenuation that exhibits an almost constant increase. 
The value of 6* & 0.1 in. at the 150-in. station for 1/, = 
The viscous flow result of 
6* & 0.3 in. appears too large. A direct check of the 
above analysis by measuring 6* at the contact front 
would verify the usefulness of the present approach to 


1.5 is a reasonable figure. 


the problem. 

It should be emphasized that at the present time 
there is no theoretical or empirical analysis that would 
account for, or unify, the various wave attenuation re- 
sults obtained by different investigators. 

For low diaphragm pressure ratios it is possible to 
discern the rear and front boundaries of the contact re- 
gion, as shown on Plate 1. As the diaphragm pressure 
ratio increases, the rear boundary merges with the 
rough flow of state (3) and, finally, at low channel pres- 
sures even the front boundary is difficult to distinguish. 
(For channel pressures below 0.3 mm. Hg the shock 
wave itself becomes indistinguishable.) Typical re- 
sults for the speed of the front of the contact region 
(U’,) as a function of Py; and x are plotted on Figs. 26 
Briefly, the behavior of the contact front is 
For 


and 27. 
opposite to that of the shock front (see Fig. 28). 
weak and strong contact fronts there is a deceleration 


AERONAUTICAL 


SCIENCES FEBRUARY, 1955 

of the initial jet of gas following the rupture of the 
diaphragm.' At the end of the formation distance an 
optimum velocity is reached. 
the optimum velocity is the one predicted by Eq. (S) 


For weak contact fronts, 


(Part (1)], and appears to be independent of distance 
(x). For strong contact fronts the velocity at the end 
of the formation distance is greater than the theoretical 
value, and the velocity further increases with the dis 
tance traversed by the contact front. 

To summarize, beyond the formation distance strong 
shock waves decelerate, whereas strong contact fronts 
accelerate concurrently. The increase in contact front 
velocity cannot be explained on the basis of the heat 
conduction solution for a “laminar’’ contact region.’ 
It is not known what effect this speed-up has on the 
flow quantities in the shock tube. 

In a real flow the boundary-layer growth is a function 
of position and time. The continuity equation indi- 
cates that the gas properties and the effective tube area 
will also be a function of the above independent vari- 
ables, and what ini‘luence this has on the contact front 
velocity is not known. 

Nevertheless, any new approach to a unifying solu- 
tion of shock-wave attenuation must consider the con- 
tact front acceleration simultaneously if it is to describe 
the real physical processes. An experimental deter- 
mination of the flow quantities between the shock 
wave and the rarefaction wave appears as a necessary 
step to aid in the formulation of a realistic mathematical 


model 


2.5) Real Flows in the Two Uniform States Separated by 
the Contact Region 

Since shock-wave attenuation, contact front acceler- 
ation, transverse waves, and eddying occurs in a real 
flow, it may be expected that deviations will occur 
from an inviscid analysis. So far, a systematic investi- 
gation of the two uniform states by measuring density, 
temperature, pressure (and particle velocity) over a 
wide range of shock Mach Numbers has not been done. 

Most of the available data on the measurement of 
physical quantities deal with flows induced by extremely 
weak shock waves. This includes pressure, temper- 
ature, and density profiles of the entire flow. 

For strong shock waves the data consist almost en- 
tirely of Mach Number measurements determined from 
cone- or wedge-type models. The results are not use- 
ful for predicting the actual values of the flow quantities 
or the flow uniformity. (Some piezo pressure profiles 
in reference 3 show that deviations of about 30 per cent 
are not uncommon for diaphragm pressure ratios of 
about 1,000.) Typical wedge measurements of the 
Mach Numbers J/. and JJ; as a function of Py, (from 
reference 19) are shown on Figs. 29 and 30. It is seen 
that the experimental results for /, agree quite well 
with the theory up to Mz 81. For A/, > 1 the Mach 
Number starts to decrease and is indicative of the effects 
the Mach 


of stronger shock waves. In region (3 
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29. The variation of Mach Number in state (2) (.V/2) with 


Velocity of the front of contact region (Us) as a FIG 
diaphragm pressure ratio (), Case: Air Air 


distance (x). 


Fic. 27 
function of diaphragm pressure ratio (Pa) and 
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Number appears to approach a limit of about J/; & 4.5 
for diaphragm pressure ratios (P4,) > 10,000. This 
may be due to the fact that the predicted low temper- 
The 


as direct 


atures in state (3) are not achieved in practice. 
real reasons may become evident as soon 
measurements of the flow quantities become available. 

A check on the temperature would have been pos- 
sible if the tail of the rarefaction wave had appeared 
on the (x, /)-schlieren records. Unfortunately, this 
characteristic line does not appear in a schlieren record 
of the rarefaction wave produced by the ruptured dia- 
phragm. However, it was possible to utilize the head 
of the rarefaction wave (Plate 1) in order to determine 
the speed of sound from a gas dynamic point of view. 
The values of the sound speed at 0°C. for air = 1,087 + 
2 ft. per sec., argon = 1,009 + 1 ft. per sec., and car- 
bon dioxide = S44 + | ft. per sec. are in good agreement 
with acoustic results.*? This also bears out the previous 
observation that for weak shocks, starting with a sound 
wave, the wave speeds are in excellent agreement with 
theory. 

It is of interest to point out that condensation effects 
may readily be observed in the rarefaction wave region. 
Plate 15 is typical of the flow at the origin for wet and 
dry air. The condensation shock path is straight, and 
its velocity is uniform at about 630 ft. per sec. relative 
to the tube. Theoretically, the path should be curved. 
The droplet precipitation as a function of time is strik- 
ingly demonstrated and suggests the possible study of 
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The Variation of Mach Number in state (3) (./;) with 
diaphragm pressure ratio (P?3,).. Case: Air !Air 
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Fic. 31. Variation of the reflected wave speed (Wai), the 


refracted wave speed (W’;,), and the shock-amplified wave speed 
(W’s:), with the incident wave speed (W, Case: HE Air 


the rate of nuclei formation from this type of record by 
using a photodensitometer. 
2.6) Real Flows in the Vicinity of an Open or Closed 
Chamber or Channel 
Plate 2 illustrated the normal shock reflection at the 
Plate 4 for 


refraction at a 


closed end of a channel for Air Air and 
He) Air. It that shock 


contact surface may also be studied from such records. 


is seen wave 
Typical results from reference 21 are shown on Fig. 31. 
The graphs indicate that the normal reflected shock 
speed approaches the variable specific heat results for 
shock Mach Numbers W2,; > 3. For incident wave 
Mach Numbers |; & 
a few mm. Hg, luminescence produced by ionization 


©, at channel pressures (/;) of 


may already be observed in state (5). 

Few direct measurements of the physical quantities 
in state (5) are available. This region will probably 
be utilized extensively for the study of atomic physics 
and The actual refracted 


wave speeds (I1’3;) deviate considerably from the theo- 


magneto-hydrodynamics. 
retical curve. The deviation is similar to that for J/ 
for the same region (Fig. 30) and indicates that the 
low temperatures may not be achieved in state (3). 
It is worth noting that in region (2), (Plate 2), which 
has the lower Reynolds Number, the shock wave does 
not appear to undergo too much of a change, whereas 
in region (3) the shock appears as broken up into multi- 
ple shock waves. The flow is typical of the interaction 
of a shock and a jet and is analogous to a flow with a 
thick boundary layer in stationary, supersonic pipe 
flow. Since state (3) has a rough and eddying type of 
flow, it produces a similar result. An additional effect 
due to the boundary layer, which is more pronounced 
in state (3) than in (2), is indicated by the refracted 
shock wave racing ahead inside the boundary layer. 
The normal reflection of a rarefaction wave at a closed 
channel is shown on Plate 3. Limited information can 
be obtained from an (x, /)-plane schlieren record, since 
only the head of the waves appears. The uniform 
region (6) has been neglected by experimentalists to 
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Fic. 32. Variation of the transmitted wave speed (W,4;) with 
the incident wave speed (W2,) for the head-on collision of two 
shock waves in air. WW = 2.18. 


date, and yet it offers excellent possibilities for research 
in a low-temperature state. 

When the channel is open the shock wave emerges 
into the free atmosphere where it diffracts and ap- 
proaches a spherical wave. It is followed by an annu- 
lar vortex. A rarefaction wave followed by a shock 
wave are also propagated back into the chamber. 
Large viscous effects are indicated on Plate 16. A 
theoretical analysis that predicts the existing wave 
configuration and the wave motion at the open channel 
is presently not available. 

Plate 17 shows the emergence of a rarefaction wave 
into the free atmosphere. A compression wave does 
not appear to be reflected back into the chamber as 
About 500 


usec. after the rarefaction has emerged from the cham- 


predicted from a perfect flow analysis.*° 


ber, an almost stationary vortex is formed within the 
tube by the inflow from the open end. (The vortices 
from the ruptured diaphragm soon merge with this 
vortex.) The flow appears to be rough, judging from 
the large number of striations, and is typical of unfaired 
subsonic entries. An adequate theory to cover this 


type of flow has not been developed. 


2.7) Experiments on One-Dimensional Wave Interactions 


The pertinent equations dealing with wave refrac- 
tion and collisions were discussed in Part (1). The 
present section will supply some basic experimental data 
on the interactions. 

Fig. 32 shows a typical result for the head-on collision 
of two shock waves. It is seen that the wave speeds 
tend toward the variable specific heat values (lower 
wave speed) as in the normal shock reflection case (Fig. 
31), which is the symmetric case of two equal colliding 
shocks. 

A typical experimental result* for the head-on colli- 
sion of a shock and a rarefaction wave appears on Fig. 
33. It is seen that the agreement between theory and 


experiment is quite good. 
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The case of shock refraction at an Air|/A contact 
surface formed by a nitrate dope film of several thousand 
ngstrom units is shown on Fig. 34 taken from reference 
sy 

On Fig. 31 it may be noted that the primary shock 
wave undergoes normal reflection at the end of the tube 


a 
2 


and subsequently refracts at the Air||He contact front 
that was produced at the origin. This problem is con 
sidered in detail in reference 21. It may be stated 
that the results for stationary contact fronts are in good 
agreement with the theory even for double refraction. 

The theory for the refraction of a rarefaction wave, 
which was outlined in Part (1), is presently being tested 
experimentally at the University of Toronto, Institute 
of Aerophysics. 

Plates 18 and 19 illustrate the interaction of a shock 
After the 


interaction the wave system consists of a reflected and 


wave with a No. 40 and No. 30 mesh screen. 


transmitted shock wave separated by a contact front 
and a vortex. The vortex is on the downstream side 
of the screens and appears on the records for higher wave 
Mach Numbers when the particle velocity is sizable 
enough to move the vortex downstream as on Plate 
18. It has been tentatively established that the Mach 
Number ./; remains quite constant for a given mesh 
size (see Plate 18 for notation). 

The schlieren photographs show that region (3) is 


at a lower density than region (2). Since the pressures 
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PLATE 17. Emergence of a rarefaction wave from the open a; 
end of a 3- by 3-in. chamber. R = head of rarefaction wave, 
V = vortices, O = open end of chamber, pj = 740 mm. Hg, pi 
Py = 370, Ty, = 298°K., and a, = 1,135 ft. per sec. Pu 
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SHOCK-TUBE 


p p», the temperature 7; > 7». This is reasonable 
since region (5) has undergone two shock compressions, 
which also results in a low particle velocity u;. There- 
fore, a rarefaction wave must exist which takes the 
high-pressure, low-velocity region (5) to the low-pres- 
sure, high-velocity region (3). The rarefaction wave is 
not seen on the records. However, along with the 
shock wave diffraction it may account for the formation 
of the vortex, which makes it impossible to determine 
state (4) from a wave speed analysis alone. 

In reference 24 a solution is given for an analogous 
interaction that occurs for the magneto-hydrodynamic 
problem when a strong shock wave collides with a 
field. 


state (4) it was not pe yssible to compare the two cases. 


magnetic Because of the indeterminancy of 

However, the excellent reproducibility of the shock 
and screen interaction may lead to its utilization in the 
study of the “turbulent” contact front and the proper- 
ties of wire gauzes in compressible flow turbulence re- 


search. 
CONCLUSIONS 


The theory of the shock tube is outlined in Part (1) 
The 
tions that are usually employed for flows with variable 


from one-dimensional considerations. modifica- 
specific heats for rarefaction waves and shock waves are 
indicated. Some of the basic results for wave reflec- 


tions, refractions, and collisions are included. For 
flows with large viscous effects such as in the diaphragm 
rupturing process, the transverse train of shock waves 
which follow the primary shock, the emergence of a 
shock wave and a rarefaction wave into free air, and 
the collision of a shock with a wire screen, a simple one- 
dimensional inviscid analysis is not possible. These 
problems still await adequate theoretical solutions. 

In the study of the origin problem a real (xv, ¢)-plane 
(see Plate 1) gives the appearance of a compression net 
of Mach lines which would result from the motion of a 
piston. Closer examination of the photographs shows 
that the head of the compression wave is supersonic 
and the region really consists of overtaking shock fronts 
which coalesce to form the primary shock wave. 

Shock-wave attenuation is a defect in velocity which 
consists of a formation decrement and a distance atten- 
Both 


viscous effects are small. 


uation. are absent for weak shock waves when 
The formation decrement 
increases with shock strength and the distance atten- 
uation rate remains reasonably constant. 

Concurrent with the shock attenuation phenomenon 
Like 
the shock front the contact front velocity is in agree- 
It is noted that a 
unifying theory does not exist which will account for 


is associated the acceleration of the contact front. 
ment with theory for weak waves. 


the shock attenuation measurements of different ob- 
servers. An adequate theory must take into account 
the apparently paradoxical behavior of the contact 
front since it occurs simultaneously with shock-wave 
attenuation. 
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The flow uniformity of the constant state regions 
separated by the contact front has not been established 
experimentally over a wide range of initial conditions. 
Neither has it been confirmed that it is justified to de 
rive the flow quantities from a measurement of the 
shock speed. These factors will have an important 
bearing on the use of the shock tube as an intermittent 
wind tunnel. 

The attainability of high Mach Number in the state 
between the contact front and the rarefaction wave 
Little data 
are available on the physical quantities behind re 


also requires an experimental decision. 


flected shock and rarefaction waves. 


One-dimensional wave interactions may be ade 
quately analyzed by the simple wave-element theory 
in conjunction with the (p, w)-plane. This has been 
adequately substantiated by the wave interaction tube 
results. 


Although emphasis has been placed on one-dimen 
sional studies at the Institute of Aerophysics, some 
notable two- and three-dimensional work has also been 
conducted at Princeton University. Yet, it is diflicult 
to state that the shock tube has proved itself as a type 
of intermittent wind tunnel on the basis of the work 


that has been done to date. The difficulties are due to 
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Fic. 34. Variation of the transmitted shock pressure ratio 
P2,) with the incident shock pressure ratio (/’;;) for the refraction 


of a shock wave and a contact surface. Cas« AirA, Ai; = 
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the lack of direct rapid response instruments for the 
measurement of the transient flow quantities. Here, 
the interferometer has already achieved wide promi- 
nence, and other instruments are being developed. 
Little is known at present about the validity of applying 
nonstationary flow results to stationary flow which 
have associated boundary-layer and_heat-transfer 
problems. 

What certainly may be said without qualifications, 
judging from the wide utilization of the shock tube for 
numerous fundamental researches, is that the shock 
tube has proved itself to be a unique instrument for the 
investigation of basic problems in nonstationary fluid 


mechanics and possibly in the field of molecular physics. 
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ABSTRACT 


A method is presented for investigating the influence of flex- 
ibility in the mounting structure of hydraulic servo systems on 
the dynamic stability quality of contro! systems for missiles and 


airplanes 

It is shown that the compliance effects can be interpreted as a 
fictitious bulk modulus of the hydraulic fluid and a fictitious 
actuator leakage coefficient Introduction of these fictitious 
terms permits treating the system as though the actuator were 
rigidly mounted 

Formulas and curves are presented which permit obtaining the 
fictitious leakage coeflicient and the fictitious bulk modulus as 
functions of the real bulk modulus and the real leakage co- 
efficient for any given overall servo loop system gain and ratio 
of true natural frequency of oil compressibility resonance to natu- 
ral frequency of the mounting structure 


INTRODUCTION 


és OBTAIN RAPIDLY FIRST DESIGN SPECIFICATIONS 
in missile and airplane control problems, struc- 
tural elasticity between the final control member and 
the actuating servomotor is often disregarded. 

With these specifications the design of the hydraulic 
valve and actuator system is then optimized. How- 
ever, it is not realized that the mount compliance 
might be so detrimental that it deteriorates the other- 
wise high performance quality of the servo system. 

This paper includes the mount compliance into the 
control investigations. It shows that the compliance 
effects can be interpreted as a fictitious bulk modulus 
of the hydraulic fluid and a fictitious actuator leakage 
coefficient. Introduction of these fictitious terms 
permits treating the system as though the actuator were 
rigidly mounted. 


Basic EQUATIONS 


The basic physical features of a closed electrohydrau- 
lic servo loop are shown in Fig. 1. In deriving the basic 
linearized equations the following assumptions are 
made. 

(1) The various elements within the servo loop can 
be represented by lumped constants. 
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(2) Only small perturbations about steady-state oper- 
ating conditions are considered. 

(3) The actuator piston is balanced. 

(4) The small masses of actuator cylinder, piston, 
and rod are neglected compared to the considerably 
larger mass of the control member. 

From Fig. 1 then, the basic equations become as fol- 


lows. 


Amplifier: 


tg = K,Veg (1) 
where 
is = differential current, ma. 
Ve = error signal input, volt 
AK, = amplifier gain, ma./volt 
Solenoid: 
fs = Ksis (2) 
where 
fs = magnetic force, lb. 
Ks = solenoid gain, lb. ma. 
Valve 


(Xp /wr®) + [(2év/ov) Xv] + Xv = fs/hi (3) 


where 

we = Vky/My = natural frequency of valve, 
1/sec. 

f& = (1/2) V by*/(kyMy) = damping ratio’ of 
valve, | 

ky = spring constant, lb. /ft. 

by = damping factor, (Ib./ft.)/sec. 

My = mass of solenoid-valve system, slug 

Xy = valve spool displacement from initial equi- 


librium position, ft. 
Hydraulic Flow: 
Ky-Xy = V 


Total _ {Compressibility 
flow rate} flow rate 


Leakage 4 Displacement 
flow rate flow rate 


B)p| + lap + a(Xup — Xac) (4a) 
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kK ry 
Kc wal S|= 
ACTUA 
tie 
FRAME |” 
3 ; ys 
Fic. 1. Physical diagram of a hydraulic servo system ) ( 
Q. = compressibility flow 
OQ. = leakage flow 
Qo = displacement flow 
Or = total flow res 
j 
From Eq. (4a) follows (r2a/Tom)p = —wmuc*Xac = wup*(Xap — 76) (5c) 
where 
lem Sache p 2&4 Tou 1 ss és 
; KyXy = Z. + + - (X ap ~~ X ac) (4b) 9 . > : 
(ra)- Pa Wy cr a wup = V (kypr?)/Ieu = natural frequency of mount Fr 
| 3 
between piston and control member, 1/sec. 
where wuc = V (Rucr*)/Icom = natural frequency of mount | 
: ; ; : between actuator cylinder and frame, | 
o, = V [B(ra)*]/( Vite) = natural frequency of : 
° og ege Sec. ‘ 
oil compressibility resonance, | /sec. : : Ce 
] BI y : : kyp = spring constant of mount between piston and ati 
& = (1/2)1,W ( Mw) (ra)*| = damping ratio A ‘ 
Sa is ay sd : As ra)*] lamping rat control member, lb. /ft. 
of oil compressi y resonance, | . : P 
1" | : oe ; get kuc = spring constant of mount between actuator 0| 
= volume der compression, [t.” : f P 
C : : € - : P , it cylinder and frame, lb. /ft. 6 
3 == bulk moc s, Ib. /it.* P : . . ' 
I - — ‘ ; ae 6 = control member deflection angle, radians, | 
p = pressure across actuator piston, lb./ft.° 
li, = leakage coefficient, (ft.*/sec.) /(Ib./ft.”) Control Member: 
e ° ° 9 | 
a = actuator piston cross sectional area, ft.” "2 : 
: . : = ee 0 2EoM * S l warp” | wh 
X,» = actuator piston displacement from initial - 6+6= - (Xap — 18) (6) 
equilibrium condition, ft. wom cm Y cM” 
X,- = actuator cylinder displacement from initial where 
equilibrium, ft. 
r = arm from actuator piston to hinge point of wea = WVkhem/Icm = natural frequency of control | Eq 
final control member, ft. member, | /sec. 
Icey = control member moment of inertia, slug ft.° few = (1/2) V bew?/(Remlem) = damping ratio of ’ 
‘ . as eli 0 
Ky = valve flow gain, ft.*/sec. control member, 1 
; 6, 
kom = any spring-type load torque on control mem- 
Actuator: : ie 
ber, lb. ft. 
Cylinder: ap = —kucX a (5a) bem = any viscous type damping load torque on 
, ; , control member, lb. ft. sec. 
Piston: ap = kup(Xap — 76) (5b) 
Feedback: 


Eqs. (5a) and (5b) can be combined into the follow- 


ing relation: 
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ny 
| where Eqs. (S) and (10) can be written as 
| feedback signal, volt | | PE, lew KX lou; 
; - a at ; . + -}pt+ p ApeAy ) 12) 
A feedback gain, volt ft. = wry’ eo ra ra 
Eqs. (4) and (5) can be combined into the following — | rK 1 wry?) [6 + 2Eeyweyd 4 
equation: wou + wy7)d 13) 
| ) : ) DE 
( 1 WiuP T OMe —Sa 
T . pt P OPERATIONAL BLOCK DIAGRAM 
\w WP OME Wy 
I ws | ane Application of Laplace transforms to Eqs. (1)—(5), (9), 
(ra)* nae ra = (12), and (13) with all initial conditions zero results in 
the following set of equations: 
Eqs. (5) and (6) can be combined into the following 
. bal 4 Vel 14) 
equation : 1s(S K,Ve(s | 
< : =e . P f etal 15) 
(0 Wey + (2Eey wery)d + 6 = (ra ky uw)P (9) . A sts\s ‘ 
{ ay i 
: - ai ~ . E 2 (S- @ + ty Wyv)S 7 1¢ Ay § 
Eqs. (5)—(7) can be combined into the following equa . E 
Io(S ky 16) 
tion 7 
.) 2bem " ra a 
| | + s+ 1) ds p(s lv) 
| rh tes 2 : ; 7" p 4 Woy” WM ky vi 
7 [wup rs (wup? + wc?) | 7 
ee l DE 
j r or 3 : WyPp we x | 10 1 o p(s 
a) 4 : n ; 
) } -tc¢ uwcyo WOM ’ yu : > a Wy” Wy Wy 
WP Wilt 
I, VW } 7 I, Vu A 
, ° ° - ° ; ° ) cAl ¢ Q 
Introducing, for simplification and abbreviation, the aa Xy(s ia" 1S) 
resultant natural frequency of the entire mount 
j V 4 S = rk r \ l Wy [s T 2E, wWwoMs T 
= Wy = WyPWiK0 V wyp- TT WV 1] (ay vie T Wy 16 ‘ 19) 
(5c) 
punt From Eqs. (14)~-(1S) the following forward transfer function is obtained 
sec, 
6(s K iA shy j l l 2E., ; , | .) DE, 
unt a ee - "Y. (S T 2& ywWous T Wom” T : T Ss oT l 20) 
me, lp(s ky(ra { Wa wy" We { wy ww 
' Combining this forward transfer function [Eq. (20)] with the feedback transfer function [Eq. (19) | in a degener- 
an . . ccs ae 
ative closed feedback loop, the following closed-loop transfer function 1s obtained. 
ator 8(s : j l | (DE, ies 2§ 
Ke /\- oe SF ie 5 (s° + 2eyweus + ecm 7, = ( ee S 
! G(s / {| \W, Wy” Wa _} Wy” Ww) 
| ne 
. o <cCOyoOuV WoM 
K, ee | 1) 
| Wy Wy Wy 
where 
(6) 
Ky (A,AsAyAr)/(Cya) = overall servo system gain, 1 /sec. 
4 ‘ Z 5 
6 = Ve/(rKpr) = reference input, radians, | 
trol Eq. (21) can be rewritten in the following form 
| 
<4 [ ; 1 
» of ' ky 1 / tn 2 s+ 1 
O\S ‘ | / Wy” Wy 
Ges ¢ 7 : c 
om- ' ly/ 1 | > &a Ky ' hr... » a 
= — 7 zs 3 = 2 = rs / -T = ie aw l { bs) T /|-tcwWcus T Wev 
As) { Wa” Wy~ Wa Wy” / wy” wy 
on | . 
5° rea 
ky | 4+ 2 + | 
S / wy wy 
+ ») 
Lyf | _— A, /[s > § - rs 
1+ iz ot ie wy 2 + ia d lL, — 2 s+ 1 ( i + 2E, uwWoemusS +t We”) 
(7) § { Ww Wy Ww Wy / wy wy 
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Fic. 2. 


The closed-loop transfer function in the form given 
by Eq. (22) allows an extremely instructive block dia- 
gram representation, as shown in Fig. 2. 


MopIFIED OPERATIONAL BLOCK DIAGRAM 


It is generally permissible to assume that the nat- 
ural frequency wy of the valve is high compared to the 
natural frequency of the structural mounts. Therefore 
it is possible to approximate the block ‘‘Mount Dy- 
namics” in Fig. 2 by the following simplified transfer 
function. 


(s/w?) + 2} [(1/2) (Ke/wm) |/ou} (23) 


This in turn permits combining the two parallel 
blocks, ““Oil Compressibility Dynamics’? and ‘Mount 


Operational block diagram of servo system. 


Dynamics,’ into one block with the following transfer 


function. 


l l E, (1/2)(Ke/wy) 
S ie ne 2 + = 
Wa” Wu Wa WM 
“Jlatali 
<s/ ei: 
( / a oa (Wy / Wy )2 | 


x. see *) : 
t j Wy 
+s Vit {oe 

“Saa \OM J/ wy ( 
[we Vi1+ (a, wu)? | 
This transfer function can be interpreted as the transfer 
function describing the oil compressibility dynamics of a 
fictitious but rigidly mounted actuator with the follow- 
ing natural frequency and damping ratio. 
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~ Ke s2 €v -4} — 
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OlL COMPRESSIBILITY 
DYNAMICS OF FICTITIOUS 
BUT RIGIDLY MOUNTED 
ACTUATOR CONTROL MEMBER 
DYNAMICS 
5. sofe 
2 +22 s? + 2€cy womSt wey os 
Fic. 3. Modified operational block diagram of servo system. 
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ae ; + 
; 7 | | | 
{ w | vit (“) W, OF LO bg—}—}_} ++} + {4} +—140 
Wy 0.9) or ee ee ae 
] ] = 0.8}— = ee ee 4 4 
= ( ) : 1a] 
a 1 + 2 + =. Oh jp 
6 Wa Wa Wg w| 
06+ see! oS a ae oe a ee Vue? 
K ) and | ‘a 
j Wa 4, , 05+-—+——+\+— +—+_+—_+. —+ + +420 ifw2 
5 ) . | \' & or 8 8 a 
j =SaWa Wy { B 0.4 + + Br + + + + + + + + 6 
ye? ye aon ae ee ae ae Z 
-—Sa —oa oF } 
* —_ | (260) 0.2+—++ + +—__+ <—t- + } + = + 8 
W, - E,W, a W, 
1 E, a a 0.1 4 ) a oe i 4 + + + + 3 
The block diagram of Fig. 2 can, therefore, be reduced 0 . : z : 
to the simplified block diagram of Fig. 3. @ 
a . 
| = Mopiriep BULK MopuLus AND LEAKAGE COEFFICIENT FIG. 4. Ratios B'/B; wa'/wa; and (1/w." wa*) versus the 
ratlO Wa/wmM 
Assuming for this fictitious but rigidily mounted BB’ = realiand Sctitions bulk modulus 
actuator the same mechanical and geometrical prop- Wa,Wa = real and fictitious natural frequency of the oil com- 
. ° ° ee pressibility resonance 
erties as for the real actuator, the natural frequencies wy = resultant natural frequency of mount 
w,’ and the damping ratio ~,’ can be obtained as a result 
1 ping g 
of a fictitious bulk modulus B’ and a fictitious leakage actuator has been increased. The resultant natural fre- 
j ireeig . ee ‘ gree 
J coefficient /,,’ quency, w,’, of the fictitious oil compressibility reso- 
The following equations hold true for the real and nance of the actuator is smaller than w, as well as smaller 
fictitious actuator. than wy, as can be seen by rewriting Eq. (25) in the 
following form. 
nsfer IB(ra)? ; B'(ra)? ~ 
ma = | Oy = (<4) , i > 2 
4 4 Wy = Wu Vi+ (Wy Wa (31) 
J al vi I al Vv 
| | | BI | B’] Fig. 4 is a graphical representation of Eqs. (25) and 
i i CM — / CM 9 : me ge : . aoe 
ie’ la 77 = Sf = 5 la Y . (=) (29), while Fig. 5 is a — representation of Eqs. 
- a\’a)- - wate)* aia . . > 
' : (26) and (30) with K = Ke/2é,w, asa parameter. 
From a comparison of Eqs. (25) and (27) it follows that A comparison of Figs. (4) and (5) indicates that the 
for the fictitious bulk modulus B’ following relations exist for any ratio w,/wy 2 O 
; a Z : 
. = = = (29) (&," Wy’) (€_/ Wa) S (I W,’”) (1/ wa" (32) 
(24) B 1 + (w,/wy)? Wa" 
Z Z 
; _ ; ons int ton KS co § 2 tw 
From a comparison of Eqs. (26) and (28) it follows that for A | or Ae = SaWa- 
. 
ister for the fictitious leakage coefficient It can now be ieee that for practical purposes only 
ofa | ; K : ver ‘an parameter values K < 1—i.e., Ac < 2£,w, should be 
low- | ln i \¢ Wa \~ Sa Sa o fs 
_ eS = . (30) considered. 
=EaWa \Wa W, | . - ———— . : 
_ . <= For wy = © Fig. 3 represents the block diagram of a 
; These two equations indicate that the structural com- rigidly mounted servo system, wg’ = wa; & = 4. The 
pliance of the servo mount acts as if the compressibility closed-loop transfer function of such a servo system 
of the hydraulic fluid as well as the leakage flow of the can be taken from Eq. (21) for ay = « 
e \ / : 9e 2 yi 
O\S) ‘ A) <$a ze - S) =—s} : ‘ 
Ke : te (s? + 2Ecuweus + &cm”) s +: adie s+ 1)+ Ke (33 
6c(s) | Wa" Wy wy wy 
If the valve is idealized by assuming wy = ©, Eq. (33) reduces to 
: / 
6(s) § 2&, val 
= Ke/ 3 a i (s? + 2ZEcmwcus + wom") + S| + Ky ( (33a 
Ac(s) / \ Wa" Wa 
i If additionally the actuator is idealized by assuming w, = ©, Eq. (33a) reduces to 
' al . 
6(s)/Oc(s) = Ke/(s + Ke) (33b) 
Application of the servo loop as part of the forward loop in a closed control loop of an airplane or missile, as 





shown in Fig. 6, gives the following closed-loop transfer function for wy = ©. 
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A(S) j N; wys ky Wy "WI , ' 
Ap S) | dD, w(S) [( a 2E Wa ) (s? r 2&; ywwows On uw") T Was (s* 2E, Wy s Wy “) + Ky Wa “wy? f 
j N; B(S) Ni w(s) Ky Wy “Wy P | , 
| D; B\S) dD, w(s) [( s+ 2E Wa) (s? + 2E, wwousS + Wem) + Wa" | is” + YE, wvsS + wy + A, Wy Wy” f 
For the idealized valve it follows that with a, 
Ay(s) j N; ws ky Wa ; ] 
Or(s | Dew(s) [(s + 2a) (8s? + 2eyweus + we") + ay’ + Kew," f 
j Nrp( S) N; wl Ky Wa | 
‘i = - : ta 
( Dep( ;) D; nr(s) [(s + 2E,wWq) (Ss? + 2&, wwous + wey + w,7s] + Kew,’ { 


For the idealized actuator and valve it follows that wit 


hwy = and w, = 


A(s) New ( S) Ky | N; B\S) N; W\s) ie 244 
—_ : + = ce) ) 
Ors) Drw(s) kh, 7 S / Dep(s) D; wis c ge 
where 

A = output variable of the controlled system 

Or = reference variable of the controlled system 

[New(s)] [Dew(s)] = net transfer function of the forward loop outside the servo loop 

|Nra(s)|/[Dera(s)] = net transfer function of the feedback loop 




















































































































Fic. 5. Ratios Ja’/la and (£'a/wa’)/(ta/wa) versus the ratio 
wa/wy With A = A,/2twa as parameter. 
la, la’ = real and fictitious leakage coefficient 
fa, 2’ = real and fictitious damping ratio of the oil compressi- 
bility resonance 
®@a, @.’ = real and fictitious natural frequency of the oil com- 
pressibility resonance 
WM = resultant natural frequency of mount 
K. = overall servo system gain 
NET FORWARD LOOP HYDRAULIC 
Gals) + OUTSIDE SERVO LOOP SERVO LOOP 
R 85 (s) 
5 New (8) TRANSFER FUNCTIONS, 
D Fw (s (33) OR(33a) OR( 3 3b) 
NET FEEDBACK LOOP 
N FB(s) 
D FB(s) 
Fic. 6. Block diagram for complete control system. 


For physical systems, the following nonrestricting 
assumption can be made for any control loop concerning 
missile and airplane control systems. 

If asa reference input a step 


Or(s) = |Op\/s (39 


is applied to the system at / 0 sec. with zero initial 
conditions for the output variable %, the true initial 


error lim s[@e(s) — 6 {s)] is equal to the initial actu- 


ating error lim s | On(s) — [Nep(s)/Dep(s) ]0o(s) { which 


in turn can be assumed to be equal to the step height 
On. This means that the product of the total for- 
ward transfer function times the total feedback trans- 
fer function of the loop has the limiting value zero for 
s— o, This, in turn, means that the degree of the 
numerator polynomial of the product is at least one 
lower than the degree of the denominator polynomial. 
In order that even the most idealized case of Eq. (384b 
reasonable control system it 
is necessary that the 
Nrp(s)New(s) is less than or equal to the degree of the 


describes a physical 


degree of the polynomial 
polynomial Drg(s)Drw(s). 
With this result and assuming for the two _ poly- 


nomials the expressions 


Nerr(s) New(s) = N,s* 4. N, fi 1 N, 9 
sm" -~?F?+,,..+ Ms?+ Ms +N N, 20 
(36 
st PD. ss" —? + 
+ Dos? + Dis + Do 


Drp(s)Drw(s) = 
Dp, as" 2 + 

the characteristic equations can be found 

They are 


from the 
denominators of Eqs. (34)—(34b). 
| 
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On Some Solutions 


of the Hodograph 


Equation Which Yield Transonic 
Flows Through a Laval Nozzle 


F. EDWARD EHLERS* 
Boeing Airplane Company and Oregon State College 


SUMMARY 


rhe properties of a one parameter family of solutions of the 
Tricomi equation are studied which have the sonic singularity 


l ~ 
(Yoly=0 ~ (1/s)* 


— () where y is the stream function, @ is the direction of flow, 


as 

and s is a function of the velocity such that s = 0 for the sonic 
velocity. These solutions represent the stream functions for 
transonic flows through a Laval nozzle. For a value of the 
parameter « = 1, the solution yields a sonic line that is concave 
to the incoming flow and has positive, nonzero curvature. The 


acceleration across the sonic line is positive and nonzero along all 
streamlines. A family of solutions yielding nozzles with dis 
continuous wall curvature is also obtained and discussed. For 
1< yp < 3/2, the flow pattern is similar to that for 1, but the 
sonic line has zero curvature at the axis of the nozzle, and the 
acceleration at the sonic point is zero for the center streamline 
but positive for all other streamlines. For 3/2 <u < 1/4, the 
sonic line is straight with all streamlines having zero curvature 
and zero acceleration at the sonic point 

For values of «» = 1 and 3/2, the hodograph sclutions are 
simple, and the computation of the throat of a nozzle is practical. 
3y means of Loewner’s transformation of the hogodgraph equa- 
tions, good agreement with the exact equations is obtained in the 
Mach Number range from 0.8 to 1.35. A method for designing 


perfect nozzles by means of these solutions is also given 


SYMBOLS 


1,38, C, D, E = constants depending on the parameter yu 

: = linearization for 1/p, defined by f’(s*)/f 

d = element of arc length 

/ = Mach Number 

t = pressure (equals unity for stagnation condi 
tion ) 


= velocity magnitude normalized to the stag 
nation sound velocity 
q* = critical velocity or the flow velocity for which 
the local velocity equals the sound ve- 
locity 
r = (6? + 45%/9 


3 pty (p/q)dq 
7 ] aq 


= 2.1120 s*/(1 + 1.2304 s* 


t = 2r/(@ — r), also used ast = 2s 3 
u, = Cartesian components of the velocity 
| = specific volume, 1/p 
y = Cartesian coordinates 
z =x-+1y 
a = }[4(—s)'/?/9] — 0} 
) = iffy — t) 


Received March 25, 1954. 
*Aerodynamicist, Boeing Airplane Co., on leave from the 
Mathematics Department, Oregon State College. 


ratio of specific heats ( 1.4 

a and 8 also used in transformation (@ t t A — 
rt =_ e@ rip 

6 = cos! 6/|(4/9 — 

: = {1.6118 s*/(1 1.2304 s* in l/p rela 
tion; also : Vn 

‘ —1.6118 s*/(1 + 1.2304 s* ‘<0 

Y —(9/4) (@" 

p density (unity for stagnation 

p critical density 

w} = 1.453820 — 0.84664 

a (1 — 1J?)/p? 9.4204s*/(1 1.23045 *)° 

iv = singularity parameter 

6 angle of flow 

¢ velocity potential 

¢’ = related potential [Eq. (5 

y streamfunction 

y’ = wiv 


(I) INTRODUCTION 


6 ew DISCOVERY THAT A GAS could be accelerated to 
supersonic velocities by forcing it through a con- 
striction was made by a Swedish engineer, De Laval. 
The first attempt to construct analytically the stream- 
lines and lines of constant velocity in the neighborhood 
of the throat of such a constriction was made by Meyer’ 
in his 1908S dissertation. Meyer expanded the velocity 
potential in a power series in x and y, the Cartesian 
coordinates about the critical point of the nozzle. 
Later, Taylor,? using a method similar to Meyer's, 
computed the transonic region in a nozzle whose walls 
consisted of ares of circles. Both of these methods 
employed series expansions, and, since there is no 
knowledge concerning the convergence of these power 
series, conclusions drawn from them cannot be relied 
upon with certainty. 

The nature of the hodograph of the continuously 
accelerated flow through a nozzle was investigated by 
Lighthill® and Frankl.* 
the derivative of the streamfunction with respect to 


Their analysis showed that 


6, the angle of flow, was singular at the critical point of 
the straight streamline and that the streamfunction is 
triple-valued in a region of the hodograph plane 
bounded by one pair of characteristics which pass 
through the sonic velocity and 6 = 0. 

Using a transformation of the hodograph equations, 
Loewner’ developed an approximate hodograph method 
that is based upon a hypothetical pressure density rela- 
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tion that gives good agreement with the adiabatic, isen- 
tropic relation in the Mach Number range from 0.8 
to 1.35. At the sonic point the values of pressure, /, 
and the derivatives with respect to the density p, 
dp/dp, d*p/dp*, and d*p/dp*, agree with the values ob- 
tained from the adiabatic isentropic relation: p = 
p’. The resulting differential equation to be solved 
for the streamfunction is the Tricomi equation in the 
transformed hodograph variables. 

In this paper we shall study some solutions of the 
Tricomi equation which can be shown to yield the 
hodographs of flows through a Laval nozzle. We shall 
investigate some of the properties of these flows as ex- 
hibited by the singularities of their hodograph solu- 
tions. 

(II) THe APPROXIMATE HopOGRAPH METHOD FOR 
COMPRESSIBLE FLow 


The equations for the isentropic, nonviscous, irro- 
tational flow of a compressible fluid are 


O¢/Ox = (1/p) (Ow/Oy) 


I 
dy Oy = —(1/p) (OW/dx) | My 


where p is the dimensionless density and equals unity 
for zero velocity. The quantities y and ¢ are the 
streamfunction and velocity potential, respectively; 
and x and y are the Cartesian coordinates of the flow. 
Chaplygin'! and Molenbroek" have shown that these 
equations can be linearized by transformation to the 
variables of the hodograph plane. With u = Owy/Ox, 
v = O¢ Oy, and Eqs. (1) we have 


dg + (t/p)dy (udx + vdy) + 1(—vdx + udy) 

= ge " dz 

The quantity g is the velocity magnitude, and @ is the 
Solving for dz, we obtain 


dz = (lIq)e” {dg + [(i/p)dy]! (2) 


flow angle. 


The elements dx and dy are easily found by separating 
Eq. (2) into its real and imaginary parts. 

Let us now consider g and 6 as independent variables. 
We then can write Eq. (2) as 


Oz Dew ( Oy i ot ) 
eq 00 | p O80 


Oz Lg ( Og t ot) 
= e + a 
Og g Og p Og 
Assuming that p is a function of the variable g only, we 


differentiate the first expression with respect to g and 
Since 0°2/060q = 


I 


and 


the second with respect to @. 
0°z/0g08, it follows that 


re) | o 
e*( - 4 = 
q Og p Og/ 


al | dy és) 
e —— - + 4 
q’ O80 dq \pq/ 086 
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Equating real and imaginary parts yields 


Og 06 = (q/p) (OW/0gq) 


O¢ fo? (2) | 2 
q =|q “— 
Og dq \p pj O86 
These are the hodograph equations first obtained by 
Molenbroek.'’ Since p is a function of g, these equa- 
tions are linear. Therefore, the principle of super- 
position of solutions may be ap; lied to find the velocity 
potential and the streamfunction in terms of the hodo- 
graph variables. Since the resulting solutions for the 
streamfunction and velocity potential of a flow are ex- 
pressed in the hodograph variables, the mapping to the 
physical plane must be obtained by integrating Eq. (2). 
For our own purposes it is convenient to introduce 
another variable defined by 0/0s* = —(q p) (0/09). 
This results in the following pair of equations. 


de/00 = — (d¥/ds*) 
0¢/Os* = w*(Op 00), w* = 3 
(d/ds*) (1 p) + (1, o°)| 


We shall approximate the quantity w* to obtain a more 
simple pair of equations valid in the neighborhood of 
the sonic velocity. 

Loewner has developed the general forms of two si- 
multaneous differential equations in two dependent 
variables which reduce to the Tricomi equation for 
one of these variables by means of a Baecklund trans- 
formation. We shall consider one such form. Let 
y and ¢ be defined by the relations 


w (Ow /Os) = —(1/a;) (O¢/08) | 
(4) 
(1 a) (O¢g Os) = SG) (Ow 06) \ 


Now consider a transformation of the dependent vari- 


ables given by 


y’ = wip | 
¢;) = —(1/a1)¢s , (5 
¢o. = —(1/a, ) Qo 4 (da ds wh 


By eliminating g’ from the preceding equations, we see 
that dw,*/ds* must be zero if the resulting expression is 


to be compatible with Eqs. (4). Hence, we must have 


w =ast+b 
Substituting Eqs. (5) into Eqs. (4) we obtain 
Vs’ = G9 ( (G 
g.’ = —sype' I 
These relations lead to the Tricomi equation 
Vs5’ + soo’ = 0 (7 


by eliminating ¢’. Now Eqs. (4) may be written in 


the same form as Eqs. (3) by the introduction of a new 
* 


independent variable s Thus 
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SOME SOLUTIONS OF THE 
oy l ds O¢ 
Os* 7 w,? ds* O¢ 
(S) 
O¢g , ds Oy _ Oy 
Oost  — ds* 08 ~——00 


By comparing these equations with Eqs. (3), we obtain 


ds Ss 
s*¥ = - + constant = : 


Ww b(as — b) 


and we have 


In terms of s*, 


from which 
— ab s*5) 


The function @ can be made to approximate w* best 
by determining the two constants so that the first two 
derivatives of w* and @ coincide at the sonic velocity. 


Thus 


dq 7 


(= ) = | (4) 
ds" J aca -_. dq ds* | 28/(28+1 


(d2a/ds?*) «-9 = —2(2y + 5) [(v + 1)/2 

From this we obtain, for y = 1.4, 
@ = 9.4204 s*/(1 + 1.2304 s*)5 (9) 
s = 2.1120 s*/(1 + 1.2304 6*) 
w, = (1.45320 — 0.84664 s) 

and Wy = (1.45320 — 0.84664 s)~! y’. The function 


¥’ for our purpose will have the same singularities as 
y since 1/a, is not singular in the region in which we 


are interested. 


(III) THrk EQuaTION OF STATE OF THE GAS LEADING 
TO THE TRICOMI EQUATION AS THE ExAcT HoODOGRAPH 
EQUATION 


In the above analysis we have approximated the 
differential equation for the streamfunction y in the 
neighborhood of the sonic point by the Tricomi equa- 
tion in a related function y’. Now from the stand- 
point of understanding the range of validity of our 
approximation, it is better to find the pressure-density 
relation that yields the Tricomi equation as the exact 
differential equation for the related function y’. We 
can then see how well this pressure-density relation 
approximates the adiabatic relation for the ratio of 
specific heats y = 1.4. To this end, we shall write the 
quantity w* in Eq. (3) in the form 
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Ww = (] p~) _— (] p- (dp ds” 


Ww = | + (dV ds*) 


Here we have replaced the reciprocal of the density by 
the specific volume |’. After substituting the approxi- 
mation for w*(s), we obtain the following differential 


equation for |’. 


(dV ds*) + V? = b8s*(1 — abs*)— (10) 
Since we require our solution in the neighborhood of s 
= 0, we shall obtain a power series for |". Accordingly, 


we substitute 


into Eq. (10). After equating coefficients of like powers 


of s* and evaluating the constants, we obtain 


1.5774 — 2.4883 s* + 8.6354 s*? 


30.175 s*3 + 88.024 s*¥4 — 223.32 s® +... (11) 


lp 


The first four terms agree with the expansion of 1p in 
powers of s* for the adiabatic pressure-density relation 
for y = 1.4. 
It is now possible to check the agreement of the 
with the actual adiabatic 
First of all we must have 


pressure-density relation 

relation at the sonic point. 
an expression for g, the velocity magnitude for our new 
gas, from which we calculate the pressure p. For this 
we make use of the differential relation between our 
variable s* and the velocity magnitude g—namely, 


—p(dq/q) = ds* 
or 
—(1/q*) (dq/ds*) = 1/pq 
Finally 
(d ds*) (1/q) = (1/p) (1-q) 


The mth derivative is 


d” l > n— | d l d” l 

wo (a) SaC mn )asen (Daven (9) 
derivatives may be evaluated successively at 
0; and a MacLaurin expansion for 1/q then may 
be obtained. The final expansion for 1/q with the 
proper numerical values substituted for p and g at s* = 


These 


s* = 


Ois 


1/g = 1.0954 + 1.7280 s* + 1.7199 s** — 

3.9340 s*4 + 6.7097 s*§ — 9.5639 s**° +... (12 
The coefficients in the expansion for |/g agree to the 
fourth power in s* with the actual expansion of 1/q for 
y = 1.4 and the series for g, of course, will also agree 
to the fourth order in s*. Now an expansion for p 
may be calculated from Bernoulli's equation—namely, 








110 JOURNAL OF THE 


y(dp/p) + qdq = 0 


y(dp/ds*) = q¢ 


The quantity y appears from normalizing p, p with their 
stagnation quantities and g with the stagnation sound 
Now 


velocity. 


= qg°(d s*/dp) 


dy on s*" 


n=O 


y(dp/dp) 


Il 


Since dp/ds* is correct to the second order in s* then 
Co, (, and Cs will be the exact values of the coefficients 
for the adiabatic equation of state. Thus at s* = 0, 
dp/dp for our new gas will coincide with the adiabatic 
relation. By successive differentiations it can be seen 
that d°p/dp? and d*p/dp* agree with the adiabatic rela- 
tions at s* = 0. If a = O in the expression for «, fol- 
lowing Eq. (5), however, the /, p relation will agree only 
up to the second derivative d*p/dp? at s* = 0. 
Knowing the expansion for g*, we are then able to 
obtain the pressure-density relation in parametric 
Now from Bernoulli's equation in differential 


i/ (5 aa) - 


where d, is the coefficient of s*” 


1/q. Evaluating the coefficients e,, and integrating 


form. 
form we have 


dp 


" ds* ss T 7 


2, ta" 
n=0 
in the expansion for 


gives us the following expansion for p — /p,. 


v(p — ps) = 0.8333 s* — 1.3145 s** + 
0.6912 s** — 0.3816 s*4 + 1.0939 s** — 


1.3948 s¥®+ ... (14) 


This relation in conjunction with Eq. (11) constitutes 
the pressure-density relation leading to the Tricomi 


equation. 
n le precec ing yaTagrapns we used a power se ies 
In the 7 ling paragraphs we used a {| r ser 
for 1/p. Unfortunately the coefficient of s*” in the 


power series increases rapidly with 7 so that the series 
converges extremely slowly. It is possible to obtain 
a solution of the Riccati equation in closed form. If 
fin Eq. (10), we obtain 


(1 + cs*)5 f” 


we let 1/p = f’ 


= ds*f 


where c = 1.2304 and d = 9.4204 are substituted for 
convenience. This equation may be simplified by let- 
ting 

f = (1 + cs*)W(x) 
where x = 1/(1 + cs*). In terms of the new variable 


x, the differential equation becomes 
W.2 = (d/c*) (1 — x) W 
Another transformation, 7 = (d“*/c) (1 — x), vields 


Win — 1W = 0 
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This is the well-known Airy equation whose solutions 


are 


A linear combination of these two functions gives us 
our general solution; thus 

f = (1 + cs*) [Aé ‘J.,,(&) + BET 
Since | /p = f’/f, we are only at lib 
A,B by boundary condi 


where & = ?/3n 
erty to determine the ratio 
tionson1/p. This is natural since the Riccati equation 
is a first-order differential equation. The specific vol- 


ume v is given by 


‘i C 
aie ‘d - 1 + cs* 
((3/2)d]“* [AE “T_2,(€) + BET, (E 
(1 + cs*)? [AE (&) + BET (é 
For s* = 0, v = 1.57744 = ad. Thus we obtain for 
the initial condition 
c + (3d) “* ([T(2/3)/T(1/3)] (A/B) = a 
or 
A/B = —[(a — c)T(1/3)]/[(3d) “T(2/3)] 
Replacing the constants by their numerical values we 
get 
] 1.2304 4 
po 1+ 1.2304 s* 
2.4176[€ 7 _2/,(€) — 0.19638 *T2/,(E) ] ” 
: (12) 
(1 + 1.2304 s*)? [E°71,,(&) — 0.1963 “7 _1/,(€) 
where 
€ = [1.6118 s*/(1 + 1.2304 s*) 


The variable & becomes imaginary for supersonic values 
of the functions 
E77 .17,(€) 


expansions in the argument 1.6118 s 


speed (s* < 0). However the 
regular power series 
*/(1 + 1.2304 s*) 
— 1.2304) 


For super- 


and é’*J:,(&) are 


which converge for any finite s* (except s* = 
and are continuous across the sonic point. 

sonic values of s*, the functions /.., must be replaced 
by the Bessel functions J , Letting = [—1.6118s*+ 
(1 + 1.2304 s*)]’” the relation for 1 p in the supersonic 


range 1S 


I 1.2304 
p ~ | + 1.2304 s* 
2.4176 [E”*J_2/,(£) — 0.1963 £7 *J2,,(é)] 
(1 + 1.2304 s*)? [—2'J.,.(2) — 0.1963 £7 J_1,(é 


(16) 


We shall now compute the relation between s* and the 


velocity magnitude g. Thus 


= —ds* 


p(dq/q) 


dq/q = —(f’/f) ds* 


th 
va 
by 
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Th 
ph 
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Finally 


q q* = f(Q))/T(s 


— 0.1667 
1 + 1.2304 s*) [€°*S1,,(€) — 0.1963E / 7 (é 


To obtain the pressure in terms of s* from Eq. (13) we 
would have to integrate numerically. However, the 
series expansion of Eq. (14) is sufficiently good for the 
range of values that we shall use. Knowing the spe 


cific volume in terms of s* we may also compute the 


Mach Number from the relation 


if? d ( ) . 9.4204 s* 
_ ak — 
p ds* \ p p- (1 + 1.2304 s*)5 (18 


M? = 1 — [9.4204 s*p?/(1 + 1.2304 s*)5] 

Thus the state quantities g, p, p, and .\/ for our hypo- 
thetical gas are given by Eas. (14)-(18) in terms of the 
variable s*. The pressure-density relation represented 
by Eqs. (14) (16) agrees closely with the adiabatic re- 
lation for y = 1.4 in the range VJ = 0.85 to VW = 1.35. 


IV) THE MAPPING TO THE PHYSICAL PLANE 


The mapping of the streamlines in the physical plane 
Along the 
Replacing 


is affected by the integration of Eq. (2). 
streamlines, dy = O and d6/ds = —yYy Wp. 
the derivatives of ¢ by those of y from Eq. (4), we ob- 
tain the following relations for the mapping of the 


streamlines. 
(Y.°/ Po) + sWe 
— (w2/g)e[vs + (so? Vs) 


9 id 
(dz ‘ds y=K = (o@i°/qgje 


(19) 
dz/d0)y=xk = 


constant) in the 


The lines of constant velocity (s = 


physical plane can be found by integrating 


dz/d@ = (1/q)e” | —ar2y. + [(2/p)Wol (20) 


Here we made use of Eq. (3) and the fact that ds/ds* = 
1°. 

From Eqs. (19) we can obtain expressions for the 
curvature of the streamlines and velocity gradient 


along the streamlines. These are 


(‘) = qv 
dl y=K wr "We" [(Ys/ Wo)? + 5 


(‘7) q 

dl ¥=K wtp Wal (Wy; We)? + s] 

dl being the element of arc length along the streamlines. 
The inspection of the above equations reveals that dis- 
continuities in the derivatives of the streamfunction 
with respect to the hodograph variables will introduce 


(21) 


discontinuities in the curvature of the streamlines and 
in the acceleration of the fluid particles. These dis- 
continuities are, of course, only admissible along the 


rIONS OF THE 


HODOGRAPH EQUATION 11] 


characteristics in the hodograph plane and therefore 


along the Mach lines in the physical plane 


E TRICOMI EQUATION WHICH 
\ NOZZLE 


(V) Tue 
YIELD TRANSONIC FLOWS THROUGH 


SOLUTIONS OF Ti! 


Some of the solutions of the Tricomi equations which 
we shall investigate are those given by Tricomi* in his 
famous paper on partial differential equations of mixed 


type. For the streamfunction we shall use 
Ls r—6 fy 6\ 

y = } a FP, ” F, bs 
wy 2r 27 J 
where r = V6? + (4s* 9), /, is the hypergeometric 
function Fyi(x) = Fl[(2u/3) — 1, (4/3 (24/3); 2/3; 


x], and yw is an arbitrary parameter. For r real, the 
line @ = 0 in the hodograph plane is the straight stream 
line Y = O, and the streamfunction y is an odd function 
of 6 and yields a flow that is symmetrical about the 


streamline Y = 0. Since the hypergeometric function 


converges only when the absolute value of the variable 
is less than unity, we need the following analytic con- 


tinuations. 


d a r—@ 
ay=r"" | -4t F,( ) T 
2r 
@—Fr _[(r-6\] 
Bu Fy | 
2r 2r 
where 


aw = {[r(5)° 6) 


B(p) = 
and 


Fx or 2 r NSN 


l—(2p/3 


wow = —(A(wE(u) — Biu)D(u) Ir x 


9 2u/3)—1 Dy 
( ) r:( )- (A (yu) F(u we 
6é-—r r—0 
Dp 1—2y)/3 9 
B(u) C(u) > #3 ( <! ) F,( 2r ) (24) 
Gg-—r rF-¢ 


where 
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Fic. 1. Lines of constant Mach Number ./ and constant angle 


6 for » = 1. 
Diy) = 
| 7-4 | /5 — 2 6—2 
PG) (Sr Cs") °C") 
a) oO / 2) 2) 
E(u) = 
2 i-4 / 4-2 > — 2p 
Gs) WA) 
} 3 } 3 
F(u) = 
2 4u — 7 / 24 — 3 2 2 
(a) a Ws) OS) 
2 ) ) ) 
and 
2yu—-2 nw-—-3 44-4 
F3(x) = F( . =, oa - = x) 
3 3 3 
o9— 2u 4— 2 10u — 4 
F(x) = F( = . = ; = *; ) 
é >) >) 
The essential triple valued property is apparent from 
Eqs. (23). If the line 6 = O for r imaginary is ap- 
proached on the correct branch of [(@ — r)/2r]'', the 


streamfunction vanishes, and we obtain the analytic 
continuation of the straight streamline into the super- 
sonic region, s < 0. On other branches, we obtain 
streamlines which cut the @ = 0 axis in the hodograph 
plane for s < 0. 

For values of u > 3/2, we shall use the following solu- 
tions found by the author instead of those given above. 


1+ 90? 
wy = 6s” F ( 4 ey 3/2; — :) (25) 
3 3 ts3 


The singularity for (¥)y-9 at s = 0 is the same as the 
solutions in Eq. (22). The reason for this choice will be 
apparent later in the discussion. 
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Let us now investigate the nature of the mapping of 


the straight streamline in the neighborhood of the criti 
cal velocity. From Eq. (19) we have, since y, = 0 
for y = 0, 


(dx/ds)y=9 = (ar7/q)sWo 


Now the solutions in Eqs. (23) and (25) have the singu- 


larity 
(Woy ges” 
Thus 
(dx/ds)yo9 = s'-* 
If wu < 2, the sonic point s = 0 maps into some definite 
I I 
point in the physical plane. If uw < 1, then 
(‘*) dq ds ' 
= es > wass—>O0 
dx y=0 ds dx 


We obtain infinite acceleration, and the solution yields 
a limiting line and cannot be continued into the super- 
sonic domain. If « > 2, then the sonic point maps into 


x = ©, Therefore, in order to obtain a transonic nozzle 
flow we must restrict uw to the range 1 <u <2. We 


shall see that we must further restrict u to less than 7/4 
in order that ¥ remain finite on two characteristics in 
the s, @ plane. 

When s is set equal to 0 in Eq. (24), the streamfunc- 


tion becomes 
r — (2/3) +1 yr 9 | 
= 1-4 T r 
v 1.4532 ( iL (*) (3) 
aus ‘ oo 
fr (2= 24) 2 (= 21) on 
3 3 f 


For uw in the range | to 3/2, the magnitude of W increases 
with increasing @. For values of uw larger than 3/2, 
y—> © as#@— 0. In the ensuing discussion it will 
be shown that, when 1 > uw < 3/2, then the line of 
constant sonic velocity in the nozzle is concave to the 
incoming flow. The choice of u equal to or greater than 
3/2 leads to flows where the sonic line is straight. It is 


) 


apparent from Eq. (26) that when uw = 3/2 the stream- 
function reduces to the trivial solution y = 0. How- 
) 


ever, wu = 3/2 results in a definite nontrivial solution 
from Eq. (24). 


(VI) PROPERTIES OF THE FLOW FOR yu = | 


When uw = 1, the hypergeometric series in Eq. (11 
can be summed, and, except for a multiplying constant, 
the streamfunction becomes 


yY = (l/a,) [(7 + 8) — (r — 6) (27) 


A section of this flow was computed by the author in 
reference 5. Most of the integration was performed 
numerically; however, the following parametric repre- 
sentation of the sonic line in the physical plane was ob- 
tained. 
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) (3/2) + [(3/16)e2] + [U1 


x = 0.964225 6° 


y = 6”* [1.4981 — 0.6078 6° + 0.0497 64 +... 


THE 
112) ) +... 


] + 0.964225 6 
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( 


‘is — 1.2509 6 + 0.2002 6? — 


0.0099 @ + 


(30/5) — (03/22) + (@/680) — 


Similarly the velocity distribution on the straight streamline, Y = 0, can be expressed by the following power series 


in s*™ 


xX = 


Because of its simplicity, the solution of Eq. (27) is extremely practical for designing nozzles. 
It can be shown by differentiating Eq. (27) that 


stant velocity are easily computed by using Eq. (20). 


vs = 


—{(1 @)) (da, ds)wy] — 1(3/2) 


2.9372 s* — 3.1042 s*2 + 3.1924 s*¥3 — 2.0035 s*4 — 2.5798 s*° + 14.0267 s**§ + 


The lines of con- 


‘wy Vo 


With this substituted in Eq. (20), the coordinates of the lines s = constant are given by 


10) dw, i@ " @y 
3—- X= — ewvt 
gq ds 2 gq 


The results of some computations of the lines of constant 
velocity are shown in Fig. 1. 

The streamfunction inside the region (4.9) (—s)* > 
#2 of the hodograph plane is triple valued, because the 
quantity 7 becomes pure imaginary. In this region 
Eq. (27) can be written 


6 2an 
(—s) ’* cos 3 > — , «e=61,3 


where 


V (4/9) (—s)* — @ 6 
6 = tan“ = Cos - 

0 (4/9) (—s)8 
A number of the streamlines are plotted in a small 
section of the hodograph plane in Fig. 2. A positive 
streamline first cuts the line (2/3) (—s) = @, then 
crosses the axis 6 = 0 and becomes tangent to the char- 
acteristic (2/3) (—s)’* = —@. Thus this latter char- 
acteristic is the envelope of the streamlines Y = A, A 
being an arbitrary positive constant; and along this 
characteristic the Jacobian J = [O(x, y)] [O(s, 4)] 
becomes infinite. The segments of the positive stream- 
lines cut off by the characteristics are mapped into the 
triangular shaped region COA in Fig. 3 containing the 
line 6 = 0 for y positive. The segments of the positive 
streamlines beyond the envelope are mapped in the re- 
gion downstream (negative x direction) of the curve 


J = for y positive. In a similar manner, the nega- 
tive streamlines Y = —K cut the characteristic (2/3) 
—s = —6 and form an envelope at the charac- 
teristic (2/3) (—2)’* = 6. The segments of the 


streamlines between the characteristics are mapped 
into the region between the Mach line COC’ and the 
line J = wo, AOA,’ in Fig. 3, while the streamlines be- 
yond the envelope are mapped into the region down- 
stream of the curve J = ©. This flow has the prop- 
erty that the fluid particles are always accelerated 


a 


¥ ie 1 dw ] *¥ io 
e wdyt+ w + e dy 


q ' ds p 


along the streamline; and the curvature of the stream- 
lines are continuous, finite, and nonzero (except for the 
center streamline). 

FLOW WITH A DISCONTINUITY IN THE 
CURVATURE OF THE WALL 


(VII) NozzL_e 


It is possible with the addition of a solution that 
vanishes on the first characteristic to obtain a channel 
in which the curvature of the streamlines and the ve- 
locity gradients are discontinuous along the first char- 
acteristic COC’ in Fig. 3. In the region outside the 
cusp of the characteristics 6 = +(2/3) , the 
streamfunction is defined by the equation 


(—S) 
v = (1/w;) (0+ 7) + 0-17) r= 
V (4/9)s* + 6 


For those segments of the streamlines which terminate 
on both characteristics, we use 





| 
«tf —_——L ee | 1 | 1 l _— 
-16 <-12 -08 -04 4" -04 08 12 
Fic. 2. Streamlines in the hodograph plane for » = 1. Lines 
C and C’ denote the characteristics, 0 = (2/3) (—s and 6 = 
— (2/3) (—s)*/2, respectively 
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Fic. 3. Streamlines in the physical plane for «4 = 1. The 
other lines are as follows: DOD’, sonic line; COC’, Mach line; 
BOB’, line of angle 6 = 0; AOA’, Mach line on which Jacobian 
is infinite 


wv = (1/o,) Real [(6 + ia) ’* + (@ — ia)'’*] +- 
Im} c/o[(0 + ia)’* — (0 — ia)*]} 
where a = V (4/9) (—s)* — 6°. The second term 


vanishes when a goes to zero, so that the streamfunction 
is continuous along the characteristic 6 = (2/3) ( —s)' 
The curvature of all streamlines in the physical plane 
are discontinuous where they cut this first character- 
istic since the derivatives with respect to @ and s are 
discontinuous. The streamlines in this region may be 
expressed more simply in the form 


Y = (2/u) (2/3) “*(—s) * V1 + & cos [(6/3) — v] 
where 6 = cos~! }6/[(2/3) 
The angle 6 varies from 0 to x for positive values of the 
streamfunction. On _ the we 


have 6 = wand é@ = —(2/3) (—s) 


(—s)’*]} and »v = tan—'c¢, 


second characteristic, 


Thus the equa- 
tion above reduces to 


, f2\" a 6 
= (—s)’*> V1 + c? cos ~—?7)= 
WwW) >) o 
Ds ss Me : 
(2) (—s) (1 + V3 0) 
Ww) 0 


In order that the streamfunction be continuous on the 
second characteristic, the extension of the positive 
streamlines beyond the point of the tangency on the 
second characteristic must be defined by 


21+ V3) /2\" re 
- (—s)’* cos 6/3 


Ww) 3 


oO 

where 6 ranges from 7 to 37/2. The streamline pattern 
in the hodograph plane is shown in Fig. 4. The flow in 
the physical plane is similar to that in Fig. 3 except 
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that curvature of the streamlines and the velocity gra- 


dient are discontinuous along the characteristic COC’, 
The quantity C must be restricted by the condition 


that the Jacobian J = (a,"/ pq”) (s¥e? + y,") have the 


same sign in the region outside the cusp of the charac- 
teristics as in the region where the streamline segments 


Now 
7 6 06 
- $1n 7 2 
3 rays) 


touch both characteristics. 


and with 


(“.) 
ds 5 constant 


= — == § 
constant ( 7 ’ 


On the first characteristic (d0/ds);5 
and the Jacobian becomes 


wr" 2 al | = 
pq” " (ss o 


, oe, 0 - . - c 
(positive function of Ss) 


J= 


tv 

+ 0O(1) > 0 
a 
In order that / go to plus infinity as @ goes to zero, we 
require that cot v > O which is equivalent to c > 0. 
(—s) 


On the second characteristic (d6@/ds), 


constant = 


and 


J = (positive function of s) X 


}—cot [(4/3) + v] a} > 0 


This requires that cot [(7/3) + v] < 0, which is equiva- 
lent toc < 73. Combining the two inequalities, we 
see that 
0<c< V3 
in order that J be positive in the region between the 
two characteristics. 
Let us investigate the effect of these limits for c on 
the discontinuity at the sonic point along the center 


streamline. For y = 0, we have 


dx/ds = (sa"/q) (We)y=o 
In the subsonic region 
(Ve)y-0 = (2/3)'” (1 


WS ) 


and in the supersonic region 


(Ve)y-0 = (2/3) (1/ms) (1 + V3 C) 
Ass—~+ 0 we obtain ds /dx = q* [w,(0)]} (3/2)'*; and 
as s > —O, ds/dx = jq*/[w(0)}} (3/2)” [1/0 + 


V3 C)]. By letting C take on the limits 0 and 7/3, 
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we obtain the inequalities 


(4 ) (4 ) | (4 ) 
= ie 
dx s—> +0 = dx a 4 dx Pee or | 


The inequality may also be written 


| dq " dq > dq 


dx supersonic dx subsonic dx supersonic 


This is the same result obtained by Frankl using more 
general methods of analysis. 

Falkovich® has solved the problem for discontinuous 
curvature along the first characteristic, but his method 
involves the solution of a complicated differential 
equation to match the conditions on the first charac- 


teristic. 


(VIII) PROPERTIES OF THE FLOW FOR VALUES OF THE 
PARAMETER | < uw < 3/2 

The flow for 1 < yw < 3/2 has properties similar to 

The streamlines in the hodograph 

they are tangent 


the flow for u = 1. 
plane have one notable difference: 


to both the characteristics 6 = +(2/3) (—s)”. To 
show this we shall first express Eq. (23) as a function of 
6 and ¢, where ¢ = 27/(@ — r). We obtain 


. [cy (us) t'™ —_— F;(—t) a 
Bi(u)t*~"”’? Fy(—2)] 


wy = [t0/(2+ 2]! 2u/s 


The quantities a;(u) and §;(u) replace the more compli- 
cated expressions for the coefficients in Eq. (23). For 
t} < 1, we have 


wy = o~°*™ Git) 


where G(t) takes the form 


Gt) = > dt" + DY et? +n 


n=0 n=0 


After some elementary manipulation, the slope of the 
streamline becomes 


52 [f° iz 
= 2+ 2 - 
ds Wo 3 6° (| r r Z 
) 6tG(t) 1 — (2u/3) |¢ G(t) 
w (2+ 0G'(t)1/ | 2+ 2+G'(t) 


dé Vs ya 


Here we have made use of the relations 
of/00 = $[2(2+ 4) (1+ A] /at} 
and 
Ot/Os = (1/3) (s?/8?) [((2 + #)/t] 
Now for small ¢, 
(1G )/[(2 + HG'()] = cto’) - 4) 4 Q(8@-YA) 


j 


Thus 


THE AERONAUTICAL SCIENCES 


FEBRUARY, 1955 


dé 2 j [ou ; 
= — s<l—c 6+ 
ds y 3 6 | L @ 
_ Zu p(tu 3 1/3) 4 { 
3 { 
Therefore, on the characteristic 6 = (2/3) (—s)°’, we 


find that 
(d6/ds)y = —(2/3) (s?/@é) = —(—s) 


which is the same as d6/ds for the characteristic 6 = 
(2/3) (—s) It should be noted that the expression 
above is not valid for u = 1, forin this case t'?’? ~ “4? = 
tand d; = —eé . Therefore the function G (¢) contains 
no power of ¢ less than 2 (except for the term of order 
one). This can be seen by expanding Eq. (27) in 
The first derivatives for 4 = 1 are all 


At the second char- 


powers of f. 
finite on the first characteristic. 
acteristic, however, the derivatives do become infinite 
and the streamlines are tangent to the second charac- 
teristic (see Fig. 2). 

The streamlines in the hodograph plane for the upper 
half of the nozzle are shown in Fig. 5 for 1 < wp < 3/2. 
In the upper left quadrant of the hodograph plane and 
the portion of the upper right quadrant bounded by the 
s = 0 axis and the curve 0 — r = 2roré@ = (1/2 
(—s), ° Eq. (23) is used to evaluate the streamfunc- 
tion. For the region containing the characteristic for 
which y = 0, Eq. (24) is used for the streamfunction. 
Inside the two characteristics 6 = +(2/3) (—s)* the 
quantity r is imaginary; consequently, the real part of 
Eqs. (23) and (24) must be taken as the streamfunction. 
For the segments of the streamlines which touch both 
characteristics, the argument of 7 is equal to 7/2 and 
the argument of 6 — r varies from 0 on @ = (2/3) (—s 
to —7 on the characteristic 6 = —(2/3) (—s)’*. On 
6 = 0, the argument of 6 — r is —7/2. For the seg- 
ments of the streamlines beyond the second character- 
istic the argument of 7 is 37/2 and the argument of(@— 
r) varies from — 7 on the second characteristic to — 37/2 
on the line 6 = 0. 
for the latter choice of arguments and therefore gives 


Eq. (23) is seen to vanish on 6 = 0 


the correct branch of the streamfunction for the stream- 
line Y = 0. The Jacobian J = O(x, v)/0(8, s) becomes 
infinite on both characteristics but only changes sign 
at the second characteristic. 

We shall now investigate the properties of the flow 
due to the singular nature of the solution on the first 
The curvature A of the streamlines is 
We obtain 


characteristic. 
found by using Eq. (21). 


- dé 
K = ~= 
dl 
(—s)' q + Olt wit —C/s | 
20*| Po!|c) (—s) { (cw! w:) 0+ [1 — (2p/3)]} 2" -" ; 


Since wp = +e 6° 24/3)(1/3) — (44/3) 4 Q(¢—-1) we find that 
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- 2u/3) — (1/3 : : 
K = o¢" ) uw = 1, then the velocity gradient does not van sh 
The velocity gradient and curvature vary similarly 


for small @ on the characteristic 6 = (2/5) (—s) aie 
= (2/3 -5) 


along the second characteristic, @ = —(2 rl 
It is apparent that the curvature vanishes as 6 goes He _ ; sea 
* é a ie The coordinates of the sonic line can be calculated 
to zero for values of u > 1/2. The velocity gradient. ; : ane 
, ae . ath from a series expansion that converges for all finite 
along the streamline is found also by using Eq. (10). : ; : ey : ; 
: ; : : values of 6. Separating Eq. (20) into real and imagin- 
We obtain . 
ary parts, we obtain 
dq ds dq uy /3) — (2/3 
= c oer ais ) 
} * ax wy l s 
dl dz ds a a cos § — YW, sin 6 
sf = dé q pq 
on the characteristic 6 = (2/3) (—s) *. Thus both IE 
the velocity gradient along each streamline and the dy a"Ws . 
; b : . somskilla: tite nis 4 sail es eel = sin 6 + WY, cos 6 
streamline curvature vanish as we approach the axis dé g pg 


of symmetry along the Mach line in the physical plane 


corresponding to 6 = (2/3) (—s) *, when » > 1. (If By differentiating Eq. (23) with respect to 6, we find that 


wy -{1 — (24/3)] r7'™ sili [A Fy\(x) + Bx *Fo(x)] + ri eS) 14 F'(x) + (B/3)x “/9) Fi(x) 4 
Bx Fs’ (x x(r + 6)/? 
where 
x = (r — 6)/2r 
When s = 0, then x = 0, and WY, reduces to 
wiWe —{l — (2u 3) 14 tie (29) 
Similarly 
l ) ») 
Wy - i . 
wy 7 y- ( = ), as ae AF\(x) + 
WwW] } 3 
1 3 [ +. B /3 , | aS? 
By *F s(x od i ine AF,'(x) + ae Sansa Fi(x) + By *Fy'(x = 
) a oo? 
Now near s = O but 6 ¥ 0, 
x = [(r — 6)/2r] = (1/2) — 1/[2 V1 + (4/9) (87/8 
x = [(1/9) (s*/87)] + O(s%) 
Therefore, for s = 0, we find that 
y, = (a’ w2)A@' “0/9)) .. [(B/3) (9 gts 30 
Substituting Eqs. (29) and (30) into Eq. (28) and integrating leads to 
w,’ 9 > » 307" — ] )" ) { , » = a l :" 
x ( “ee » & + ve Be) “4S > 
q* n=0 (6 — 2u + 67) (2n)! q* 0 (4 - 2u + On) (2n 
9 ap2n—1 
a I , ) = a1 (—] 
| - )( , ania » 3 
3 wp*q* 0 (9 — 2u + Gn) (2n 4 ! 
a — 30°" F1( 1 
\ a— ee 4 
q* n=0 (9 — 2u + Gn) (2n + 1)! 
, 7 = gently _y 
() B (9)*'"9t- %/) > — 
q* n=0 (7 — 2u + Gn) (2n + 1)! 
2 i ae 30°"(—1 
bom = ( Ag -@a/) > 32 
3 w,p*q* 0 (3 — 2u + Gn) (2n)! 
Now the slope of the sonic line near 6 = 0 is 


dy ly /dx [ 2 l — Goa w,B rr oe 
a | 7 (: - ‘) \ 46 aa!) 4. Og") — (9)*/ig't—20)/3 4 gg! — 4) | 
an dé / de { 5 aip*q*/ / 3q* | 
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d y 
dx 


Since A < 0, B > 0, and 1 < uw < 3/2, then the slope is positive for 6 > 0. 
term of Eq. (31) and of Eq. (32) we see that « and y increase with 6. 
Now the curvature of the sonic line near the origin is of 


the incoming flow which is in the negative x direction. 


the order of 


AERONAUTICAL 


= — {[(3 — 2u)A]/[w(9)'"B]} 0 
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+ 0(6@°*) 


From consideration of the lowest order 
The sonic line is then seen to be concave to 


dx / I+ (“) | d (*") dé } / [ (“) ((9 7 (2/3) + (2u/3 
= | + — ( - ss 
dy?/ dy do \dy/ dy '/ y' dy 


The curvature of the sonic line vanishes where it crosses the axis of symmetry of the nozzle for u > 1. 


(IX) THE FLOW THROUGH A NOZZLE WITH STRAIGHT 
Sonic LINE 3/2 < p< 7/4 


9 


For values of wu greater than 3 
in Eq. (23) could be used to construct flows through 


2, the solutions given 


nozzles in which the sonic line is straight. However, 


we shall establish solutions that are more convenient. 
We shall assume the function y’ can be expressed in the 
form 

vy’ = 6s “F[—(9/4) (62/s*)] 
If this expression is substituted into the Tricomi equa- 


tion, the following differential equation results. 


nil — n) FP” + 


3 f 2 ( 1) 
E a (; ¥ | il mat F=0 


where 7 = —(9/4) (0°/s*). This is the hypergeometric 


equation whose solution is 


F ire pp 8. .) 
_ *s3°" = 4 53 
»>/y 


Unlike the solution given in Eq. (23), the value of 3/2 
for u does not lead to a trivial solution but yields 


e335 2 3 9 @ 
fk .=5 = 
s 6 2 Z i s? 


This hypergeometric function can be summed; for, 


y’ = 


with 
(9/4) [02/(—s)3] = & 
and 
rf((3/2) + a) = [(1/2) + nlrf(1/72) + 2) 


we have 


7h. 4 @ I r[(5/6) + ne" 
F(S. : : y) = 2d ee = ; 
3 a — n=0 T(5/6)(2u + 1)n! 


The quantity on the right-hand side can be recognized 


Str 


Finally, for the solution to the Tricomi equation we 


have 


) 


230 
y’ = } (1 — e)~° de 


/70 


after neglecting a constant multiplier. If the trigono- 


metric substitution € = sin ¢@ 1s made we obtain 
®sin [36/2 
yy’ = } sec ‘6 do 
/7 0 
For subsonic flow, € = if where / is real. Then 
730 /2 
yp’ = } (1 a t ) dt — 
~/ VU 


®tan 30,2 
/ sec “¢ do 34) 


The negative sign is chosen in Eq. (54) so that the 
streamlines in the physical plane form a converging- 
the 
streamfunction may be expressed as an incomplete 


diverging stream. For supersonic region, the 


Beta function 


y’ = -2 'B,(1/6, 1/6) 
where 
| 30 
x — 
y } ( —S§) 
ex 
B,(p, q) = } x?—' (1 — x)" 7! dx 
70 


For values of uw other than 3/2, the quantity s “ be- 
comes complex on passing into the supersonic region. 


For s negative and 7 < 1 we choose the branch 

vy’ = 0(—s)-* cos n ur F(n) (35) 
where v is an odd integer such that cos num is negative. 
With this choice we obtain the hodograph pattern shown 


This represents a nozzle in which the ve- 
If we 


in Fig. 6. 
locity is accelerated to supersonic velocities. 
choose another branch for the transition from the sonic 
point back into the subsonic region we obtain the pat- 
tern shown in Fig. 7 where the lower branch is obtained 
from 

y’ = 6s" cos 2n urF(n 36) 
where ” is an integer such that cos 277 < 0. This 
gives a channel that has a subsonic exhaust and in 
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which the sonic velocity is attained on a line across the 
throat of the channel. 

The origin in the hodograph plane is a branch point 
of the streamfunction. Approaching the origin along 
the different lines in the hodograph plane we obtain 
different values of y. To illustrate the point, we write 


the streamfunction as 


= (9 Jan) (=a 


For up > 3/2, 7 goes to O, and s *" goes to infinity 


as 6 and s go to zero, so that y remains finite. Thus 
near the origin the streamline y = A is the path 
figs °* —n) = 36 Ss ef [2 w,(0) | = 
K + O(s** 
For s < 0, the streamfunction 
l woe - - 
y (1 — &) dé (34 
Wy JJ 
becomes tangent to the characteristic, (38/2) |@/(—s) 


= |. The streamfunction may be continued in two 
ways. In the region outside this characteristic the 
integrand becomes complex. The streamfunction in 
this region becomes 
7] 
| = eae 
y = (1 — &?) dé + 
W@W) « 0 
cos [(57/6) (2n + 1)] iki F : 
(€ | dé 
Ww) 7 | 
tae 
. vl) + 
Ww] 
cos [(2/3) (n — 5/2)] (” i 
(g- — ] dé 
wy J 1 
This streamfunction has two branches. For n = 1, 
the cosine vanishes. For » = 0, cos [(7/3) (n — 5/2)] 
< 0 and the flow continues expanding. For n = 2, 


cos [(7/3) (n — 5/2) 
giving a combination nozzle and supersonic diffusor with 
If the streamlines are continued back 


into the region 30/[2(—s)’*] < 1, the streamfunction 


> O and the flow is decelerated, 
a second throat. 


becomes 


For n = 3 we obtain the orginal streamfunction, which 
vanishes for @ = 0. For 2 = 2 we obtain cos (27/3), 
which is negative, giving a nozzle and diffusor com- 
bination. For n = 0, we obtain the expanding flow, 
although, less expanding than that with » = 2 for the 
continuation into 3 6/2(—s)°*> 1. 

For values of w other than 3/2, the analytic con- 
tinuation of the hypergeometric function into the re- 
gion around the characteristic 36/2(—s)’° = 1 is given 
by 


H E 
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In order that y be finite on the characteristic, we see 


that u< 7/4. Since (7/6 (2u/3) < 1, the derivative 


of }1 + (9/4) (@*/s?) } 


”’ becomes infinite on the 


characteristic, and the streamlines are tangent to the 


characteristics 1 + [(9/4) (@°/s*)| = O. 


(X) PROPERTIES OF THE FLOW 


The physica! coordinates of the flow may be found by 
integrating Eqs. (8) and (9) with the streamfunctions 
given above. However, it may be more convenient 
to express Eq. (2) in terms of y and ¢’ by using Eq. 


(5). This yields 


( wa,’ qe" ye] + 
[(idp/p)e"] (37 


dz = = [ (wy, gje"dy’ = 


To find ¢’ we note that, by Eq. (6), 
y,’ = —0s""~* (uF + 3nF’) = 
or 


ge = (2/9)s *" 


> (uF + 3nF’) (0n/00 


Integrating we obtain 


») » 
2 jea7v7 1 o 3 
9 u] F ( 9 ’ =? 5? ") 
» ” 2 - 


Following a similar procedure for nu = 3/2, we have 


. gy i] 
y; = u~* 
$ (—s) 
: l On 
Ye = 2 (—s)’* (1 — ») 50 
and finally, 
gy’ _ —3(—s)' (1 = 7 


We note that ¢’ is identically zero on the sonic line of 


the physical plane. Thus, since s = 0, 6 = O on the 
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sonic line in the physical plane, then d8 = 0 and dg’ = 
0; and consequently, Eq. (37) reduces t« 


(dz),=9 = (idW)/p* 


Thus the sonic line is straight and the streamline yx, 


say, cuts it at 
x= 0 
y = vx/p* 


The streamlines then can be constructed by integrating 
Eq. (37) numerically from the sonic line with y = Wx a 
constant and with the values of 6, s, g, p, and y’ deter- 
mined by ¥(9, s) = Wr. 

Let us investigate the curvature of the streamlines 
where they cut the vertical sonic line. To use Eq. (21) 
for the curvature, we need the slope of the streamlines 
ats = 0 


(*) v, Os *—" (uF + 3nF’) + w/0s~" F(n) 
¥ 


namely, 


ds a Wo s”* (F + 2nF’) 


To approach the streamline px at the origin, we must let 
6 and s go to zero along the path 6 = As* for up > 3/2, 
K being a constant. We find by Eq. (21) and the rela- 


tion above for d@/ds that 


(5) ; 
= Q(s” )—->O 
dl ¢=WK 


Thus the streamlines have zero curvature where they 


as sO 


intersect the sonic line. 

In a similar way, it can be shown that the acceleration 
of the fluid vanishes at the sonic line. From Eq. (19) 
we have 

ds dq ds qd 
az = i) l 10 | 2 a 
= = - 41/po| (d0/ds)y? + s]} 
dl dl dq aw,” , 
Therefore the acceleration of the fluid along a stream- 


line at the sonic line is 


_ dq : dq dl 
lim = lim = 
0 @ »—~0 a at 
: dq ; , 
lim g = hm 0(s* *) = 0 
nO d s— (0) 
and 
dq 
—>U ass — 0 
dt 
Here, the limit was taken along the path 6 = As“ as 


in the foregoing explanation. 
It is of some interest to find the curvature of the line 
of constant velocity where they cross the axis of the 


nozzle. For this purpose, it is convenient to expand 
y in powers of 6. Thus, 
6 ww +1). : 
y= _ ra 6° + 0(6) 
wS” 60,5" *" 


Differentiating with respect to @ and s, we obtain for the 
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two derivatives 


l u(u + 1) 
ee a 8 


Mu l 1 dw 
¥, = (- qG-s =) 6 + 0(6°) 
w 5" @)~ s” ds 


Now the sign of the curvature at the axis of the con- 
stant velocity lines is given by the sign of d?x/dy 


[(d/d0) (dx/dy) | (d6/dy). 
dx _ (=) (3) _ par*p, cos A + YW sin 6 


From Eq. (28), 


dy dé / dé pw", sin 6 — pW, cos 0 
3S) 
dx | pw, *(d0 ds)y — tan? 
dy [ pw," (d0/ds),| tan 8 + 1 
Using 


Pe 
‘S 
+ 

© | 


(“) (*:) A 
= — = 6+ 0(6 
ds y Wo ds 


differentiating dy/dx with respect to 6 and dividing by 
dy wp, sin 6 vo 
— _ _— cos @ 
dé g pg 
(1/pgais") + O(8) 


we obtain 


(S 

dy 6=0 
fors > 0. 
the lines of constant velocity at the axis are concave 


“ Piw~ : da ee 
pqwis rT Wp et l pl 0 
Ss ds 


Since the flow is in the negative x direction, 
to the incoming flow. This concavity disappears as we 
approach the sonic point s = 0 if uw > 1. 

Using a similar analysis for s < 0, we obtain 


d*x 


- = — pgu,(—s)" X 
dy" 


ps M da 
pw) — po / + 1 ]cos nu > 0 


L (—sS) as 


Thus the streamlines are concave to the incoming flow 
in the supersonic region as well as the subsonic (see 


Fig. 7). 


(XI) CONSTRUCTION OF PERFECT NOZZLES BY MEANS 
OF THE HODOGRAPH METHOD 


In the preceding pages, we have investigated solu- 
tions of the Tricomi equation having the required 
singularities to represent the transonic flow region of a 
Laval nozzle. In addition to one of the fundamental 
solutions discussed previously, we may also add other 
solutions that are regular in the region of the s, @ plane 
in which we are interested. If we wish to design an 
inlet that starts from a uniform parallel subsonic flow, 
then we must introduce a function with an arctangent 
singularity at some point so, # in the hodograph plane, 
representing the angle and velocity of flow at infinity. 
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SOME 


A development of such a function is sketched briefly by 
usable 
shall 


Carrier and Ehlers. However, a more form 
was obtained by Falkovich® 


By the substitution 


which we derive 


simply. 


the Tricomi equation transforms to 
/ — 
Wet = i Veo = 


Falkovich introduced the orthogonal transformation 


[(1/3t)p.’] = O 


(0 + it) + ito]/[(@ + it) — ito] = e 


Since this is the bilinear transformation, the lines a = 
constant and 8 = constant are circles in the (¢, 6) plane. 
Under this transformation the differential equation 


becomes 


l 
ta 7 Ves <a xX 
3 sinh a (cosh a + cos B) 
[ Oy’ ; Oy’ 
Mis cosh @ cos 2B) + sinh a sin 8 = 0) 
a Oa 3 


Since we desire a solution with an arctangent singu- 


larity, we assume a solution of the form 


With the further substitution that cosh a = 1, 


(cosh a + cos B) * Bx(a 


the dif 
ferential equation for x reduces to the hypergeometric 


equation 


7 dx l 
dr ? 6 dr 144 


7 l 7 l 
F( , an &. 
i ts T 


Thus with 7 and 6 replaced by their equivalent relations 


l l l \ 


12’ 12? 2’ J 


in terms of s, 6, and so, we obtain for the streamfunction, 


using the second solution of x, 


1 jf4 ine { } 
y ef 9 Sy G- + 9 (S° = Se { x 
| 7 5 | ‘|? 
I - 1: |) ge 4 (5/2 — 5,°/?) | 
12° 12 { 9 


( | E Vs 
g2 + s3 + Soo |} 7] X 
( f 


[ t so 
l (39) 
g@- + 4 (s 9 J 


By the proper choice of the branch of the above many- 


6/3 
tan 
— ¢,3) 


valued function, it is possible to make all the stream- 
lines in the hodograph plane converge to the point @ = 
0, s s) and thereby prescribe a uniform flow at 
with a velocity given by so with zero angle of inclination. 
Thus the nature of our channel may be changed by 
varying our choice of Sp. 

In addition to the two singular hodograph solutions, 


we may consider solutions that are regular. We must 


SOLUTIONS OF THE 
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require, however, that y = O for 6 = O and that y be 


an odd function of @ in order to obtain a symmetrical 


channel with a straight axis. Some possible solutions 


are those obtained by separation of variables 


namely, 


l 2 
y, = § J ~ as sinh aé 
Ww] ) 
l ; 2 : 
y = tafe —~ sin af (40) 
Ww) ” 
V3 = 0/0 
V4 = 0s/w 


These six solutions are analytic everywhere in the 

We may also construct a family of solutions 
of the Tricomi equations by noting that, if yy’ is a 
solution, then yo’ = — S v'd0 is a solution if we add 
an appropriate function of s. Starting with y’ ) 
and y’ = @s, we obtain the following family of solutions 


6 plane. 


by integrating with respect to @ and adding a term in 


s whenever necessary. 


y,’ = 0 

y.’ = @ — [(1/3)s 

y;' = s*0 A 

yp," —s'9? + (64/2 s8/15 

Ws’ —- —[(539 3] + (69/10 6/15)6 

Wy! _ $1) 
po’ = s§2? — [(1/6)s 

y ‘= sO? + [(1/2)s46 

y;’ = so? — 50 1/21)s? 

v5) = —[(s6*)/5] + [(s58%)/3 1 /21)s‘0 


Only the odd numbered solutions are suitable for com- 
puting flows in a symmetrical channel. These solutions 


were first given by Tricomi in the form 


—, © 7 


where 7 is a positive integer, F is the hypergeometric 
function, and r = V[(4 We 
construct a solution of the hodograph equation which 


9)s*| + 6°. may thus 
represents a flow through a channel by taking a com 
bination of three types of solutions [Eqs. (22), or (25), 
(39), and (40), This gives us several param 


eters by which we can vary the channel shape or the 


or (41 


axial velocity distribution. Integrating along s con 


stant, the position coordinates of the point (s, #) are 


Thus, if y (9, s 


position coordinate 2, and y(6, s 


additive. is one streamfunction with 


is another with posi- 


= 


tion coordinate 22, the coordinate of y; = y yo is then 


21 T 22. 
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Fic. 8. Schematic diagram for terminating channel in a uniform 
supersonic flow. 


(XII) TERMINATION OF THE NOZZLE BY A UNIFOR) 
FLow 


It is possible by means of the method of character- 
istics to terminate the nozzle in a finite distance by a 
uniform flow. The transonic portion of the nozzle may 
be constructed up to some curve g = a constant in the 
supersonic region by means of the analytic methods 
presented in the previous sections. Thus along this 
curve the quantities g, 6, and / are known. The choice 
of g must correspond to a Mach Number greater than 
about 1.10 since the method of characteristics is not 
so accurate for lower Mach Numbers. We must also 
specify that the velocity must be tangential to the center 
axis AG in Fig.8. With these initial and boundary con- 
ditions, the flow may be constructed in the quadrilateral 
region ABCD, where BC and CD are Mach lines. 
Consequently, the channel walls may be found by con- 
tinuing the selected streamline through this region. 
The extent of the region ABCD is determined by the 
choice of Mach Number for the uniform flow, which is 
the Mach Number corresponding to the axial point D. 

To obtain a uniform flow beyond the point D, we 
make use of Prandtl-Meyer flow, which consists of one 
family of straight Mach waves. From the character- 
istic CD we construct the family of straight lines, hav- 
ing the direction of the Mach angle as measured from 
the direction of flow at each point along the character- 
istic CD as shown in Fig. 8. Along straight character- 
istics in the physical plane for two dimensional flow the 
quantities g, 6, and J are constant. The boundary 
streamlines may be extended by drawing straight line 
segments from each straight Mach line, having the di- 
rection of the velocity vector which corresponds to 
the point of intersection of the straight lines on the 
Mach line CD. Drawing a smooth curve through the 
vertices will yield the channel wall. Beyond the Mach 
line DE, the flow will be uniform. The accuracy with 
which the flow is calculated in the region DEF de- 
pends upon the number of straight lines drawn. 
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ABLE | 
Tabulation of s, 1/q, g, and 1/p 
1/4 9 ] 
0.34 1. 45156 0.68891 1. 27356 
0 26 1.35002 0.74073 1. 3336 
0 18 1.26118 0.79291 1.4980 
O.1¢ 1.24065 0 80603 1.41545 
0.14 1.22071 0.81920 1.43333 
0.12 1 20135 0.83240 1.45175 
0.10 1.18252 0.84565 1.47078 
0.08 1.16420 0.85896 1.49050 
0.06 1 14636 0.87233 1.5109 
0.04 1. 12897 0.88576 1.53223 
0.02 1.11200 0.89928 1.55437 
0.00 1 09544 0.91288 l 7744 
0.02 1.07927 0.92655 1.60151 
0.04 1.06345 0.94034 1. 62666 
-0.06 1.04797 0.95423 1.65293 
0.08 1.03282 0.96822 1.68042 
0.10 1.01796 0.98236 1.70918 
0.12 1.00339 0.99662 1.73929 
0.14 0. YS8908 1.01104 1. 77084 
0.16 0.97502 1.02562 1. 80387 
0 18 0 96120 1.04037 1. 83818 
0.22 0.93420 1 O7044 1 91318 
0.27 0.90152 1.10924 2. 01693 
0 32 0 Si986 1.14961 2.13391 
TABLE 2 
Tabulation of s, s*, and 

0.34 0.20075 1.16534 

0.26 0.14508 1 24307 

0.18 0.09521 1 30080 

0.16 0.08355 1 

9.15 0.07782 1. 

0.14 0.07217 r4 

0.12 0.06109 . 3 

0.10 0.05028 a 

0.08 0.03973 Bes 

0.06 0.02944 l 

0.04 0.01939 1 

0.02 0.00958 l 

0.00 0.0000 1 

0.02 0.00936 ] 

0.04 0.01851 1 

0.06 0.02745 1.50400 

0.08 0.03619 1.52093 

0.10 0.04474 1.53786 

0.12 0.05311 1.55480 

0.14 0 06129 1.57173 

0.16 0.06930 1. 58866 

0.18 0.07714 1.60560 

0.22 0.09233 1.63946 

0.27 0.11046 1.68179 

0.32 0.12077 1.72412 


(XIII) CONCLUSION 
In the preceding discussion, we have studied a group 
of hodograph solutions which have the singularity at 


the sonic point of the form 


(Wey =o (1/s") ass—~0O 
and which yield flows through symmetrical channels. 
For u = 1, the sonic line is concave to the incoming flow 


and has finite curvature. In the neighborhood of the 
throat, the acceleration of the gas along a streamline 1s 


positive and the curvature is nonzero, except for the 


streamline of symmetry. For | < uw < 3/2, the flow 
pattern is similar to that for « = 1, with the exception 


that the sonic line has zero curvature where it crosses 
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SOME SOLUTIONS OF THE 


the straight streamline and the acceleration of the gas 


along the axes vanishes at this point. For 3/2 <u < 


7/4, the sonic line is straight and the acceleration of the 
gas along the streamlines vanish at every point of the 
sonic line. 

For » = 1, a simple procedure was set up for con- 

structing a channel in which the curvature of the boun- 
dary is discontinuous at a point in the supersonic re 
gion. A similar procedure could be used for 3/2 
p<. 
It should be noted that for u = | and 3/2, the solu- 
tions are relatively simple, and computations of the 
boundary of the throat of a nozzle are practical. Since 
the differential equation for y’ is linear, the flows may 
be modified by adding to the singular solutions the 
solutions in Eq. (41) which are antisymmetrical in the 
flow direction and are regular in the neighborhood of 
sa 
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A Method for Investigating the Influence of Flexibility of the Mounting Structure 


(Continued from page 106 


De Je De : 
S —SvOyvy TT 2SaWa T =<SCM@CcM T 
Dp se" FF 4. 2... =e ye 
S + (2&0 + ZEcuwou + D, Js* * * + 
= (37a 


+ [Kce(1 + N,) + D, -1J]s" +... =0 (87b) 


The sum of the real parts of the roots, which describe 
the damping quality of the control loop, is given by the 
negative value of the coefficient of the second highest 
term in any characteristic equation. 

In order that the increasing degree of the character- 
istic equation with increasing complexity of the servo 
loop dynamics is not necessarily accompanied by a de 
crease of the damping quality of those roots that are 
characteristic for the preceding less complicated servo 
loop, the following relations should be guaranteed. 

2&,W T DE, uWemu) > Ky (1 + N,) - Ke 
>O (38S 


9¢ ; 
-<syoy A 


al 


On the other hand, in order to keep the power require- 
ments of the servo system as low as possible, it 1s neces- 
Sary, in a good design of a control system, to keep the 


term 2£-);we,, as small as possible. This means that, 


for a satisfactory design, the relations (36) should be 
replaced by 
2E,W_ > F Ptvw 0 38 
ZeaWa > Ke Ztyoy > Sa 
and this, in turn, finally means 
K = Keo/ (Zea) < 1 39 
) 


Thus, for all practical cases, the relations (32) should 


be replaced by 
6 hy w,,’ (£ Wy & (1 Ww l wW (40) 


for (w, wy)” > 0. The destabilizing term s w,’* in the 
load feedback loop of Fig. 3 thus becomes more and 
more dominant over the stabilizing term 2£,’ w,’ in the 
saine load feedback loop. 

From these considerations it can be seen that for a 
large variety of control systems the performance quality 
of an otherwise excellent hydraulic servo system can 
noticeably deteriorate due to the elasticity of the 
mount. It is of no use to improve the actuator dy 
namics without due consideration of the compliance of 


any structural mounts 








On the Wake and Drag of Bluff Bodies’ 


ANATOL ROSH XO7# 
California Institute of Technology 


ABSTRACT 


A modification to Kirchhoff’s free streamline introduces the 


parameter k = Y1 — C,,, which allows arbitrary base pressure 
and which must depend on the dynamics of the wake. For a 
cylinder of given cross-sectional shape, the drag, Cp, and the wake 
width, d’, are functions of k only. These functions are used to 
relate Cp and the dimensionless shedding frequency, S = nd/U 

to another number, S* = nd’/U,, which is based on wake param 


It is found that S* = 0.16 for all cylinders. In another 


eters 
approach, k& is evaluated by using Karman’s solution for the 
vortex street 


SYMBOLS 


( = distance from back of cylinder to trailing edge of inter 
ference element 

Cp = drag coefficient 

Cys = base pressure coellicient 

d = cylinder diameter or breadth 

d’ = distance between free streamlines 

h = width of vortex street 

k = base pressure parameter (U,/U,, or VY 1 — ¢,, 

l = longitudinal vortex spacing 

n = vortex shedding frequency 

R = Reynolds Number based on cylinder diameter ( U’,,d/v 

R* = Reynolds Number based on wake parameters (l’,d’/v 


S = evylinder Strouhal Number (nd/U’,, ) 


S* = wake Strouhal Number (nd’/U’,) 

u = velocity of vortices relative to free-stream velocity 

U’., = free-stream velocity 

U, = velocity on free streamline at separation 

[ = circulation per vortex 

€ = fraction of shear layer vorticity which goes into indi 
vidual vortices 

\ = wake width factor, h/d’ 


(I) INTRODUCTION 


i jew PROBLEM OF THE FLOW about a bluff body is 
returned to, in the literature, again and again. 
The reason for this is clear enough —not only is the 
problem always of great practical importance, but it has 
so far not given way to theoretical treatment. In 
fact, there has been little new advance in the theory dur- 
ing the past 40 years, the contributions of Kirchhoff 
Karman still remaining the cornerstones from 
A large part of the 


and 
which almost all writers start. 
theoretical work since then has been little more than 


Presented at the Aerodynamics Session, Twenty-Second 
Annual Meeting, IAS, New York, January 25-29, 1954. 
* The experiments reported were performed in the 20- by 


The major portion 


1) 


in. low-turbulence wind tunnel at GALCIT. 
of this work was under the financial sponsorship of the National 
Advisory Committee for Aeronautics. The author is indebted 
to Dr. H. W. Liepmann, under whose program of general turbu 
lence research the present investigation was conducted 

+ Research Fellow, Department of Aeronautics 


elaboration and extension of those theories. But, as 
already pointed out by the early investigators, there 
seems little hope of arriving at a complete theory be- 
fore something is understood of the essentially non 
steady processes in the wake immediately downstream 
of the cylinder—1.e., the dynamics of the vortex forma- 
tion. It is unfortunate that more of the energy has 
not been applied in this direction. 


of the problem does indeed appear formidable, and one 


However, this part 


would like to progress a little further without having to 
tackle it directly. 

Possibly not enough attention has been paid to the 
experimental information that has been accumulated 
and which is indispensable for the realistic formulation 
of theoretical models. The work of Fage and his asso- 
ciates" °? can be singled out as a particularly fruitful 
source of experimental results that have not been fully 
exploited. 

It is the purpose of the present paper to give an over- 
all view of the aspects of the problem which must be 
considered, and combined, in any complete theory. 
These are, principally, the potential outer flow and the 
wake. It is shown that a realistic representation of 
the potential flow, in the vicinity of the cylinder, may 
be obtained by a modification of Kirchhoff’s method. 
This depends, however, on a parameter that can be 
The 
crucial region in the ‘“‘joining’’ of the potential part of 
the solution to the wake is the region of vortex formation 


determined only from a consideration of the wake. 


just behind the body. Some experimental results are 
presented to demonstrate how critically the whole prob- 
lem depends on that region. Nevertheless, even with- 
out understanding the vortex formation, the modified 
Kirchhoff theory may be combined with some con- 
siderations about the wake to obtain a solution. This 
is still semiempirical, but the empiricism is a mint 
mum, and a useful unification of different bluff body 


shapes is obtained. 


(II) MopIFIED KIRCHHOFF METHOD 


The main feature of flow past a bluff body is its 
separation from the body surface, well ahead of the 
rear stagnation point, and the formation of a large 
wake. ‘The presence of the wake alters the flow and the 
pressure distribution on the body, as compared to the 
potential solution, resulting in a deficit of pressure on 
the downstream side and an excess in front. This gives 
a net pressure drag quite distinct from that due to fric 


tion. 
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THE WAKE AND 

In the approach taken by Kirchhoff, the problem is 
to find the drag by calculating the pressure distribution. 
For this, one has the experimental observation that the 
separated boundary layer continues downstream as a 
free shear layer, which is well defined at first and forms 
a boundary of demarcation between the wake and the 
(distorted) potential outer flow. The problem con- 
sists then of finding the shape of this boundary (which 
is idealized as a surface of velocity discontinuity, or 
so-called free streamline) so that the potential outer 
flow, and thus the pressure distribution, may be calcu- 
lated. 
that the pressure everywhere in the wake is the free- 
stream pressure, p,. This then is also the pressure on 
the ‘‘base’’ of the body, downstream of the separation 
The 


magnitude of the velocity at the separation and on the 


The basic assumption in the Kirchhoff theory is 


point, as well as all along the free streamline. 


free streamline must then be the free-stream value, 
U,. Thus, in the hodograph plane, the free stream- 
line is simply a circle and permits an elegant solution. 
However, the calculated drag is considerably lower 
than what is observed, the discrepancy being directly 
due to the assumption about the base pressure, which 
is observed to be actually much lower than p,. On 
the other hand, the measured pressure is remarkably 
This, and the fact that free 
streamlines appear to be good representations of the 
free shear layers, indicates that the Kirchhoff method 
should be sound, provided more realistic (lower) values 
But then the pressure 


uniform over the base. 


are chosen for the base pressure. 
in the wake can no longer be uniform, for it must even- 
tually become p., far downstream. On the other 
hand, there need not be too much concern about the 
details of the solution far downstream, for in any case 
the free streamline representation will be valid for only 
a few diameters downstream of the body. The main 
aim should be to find a model that will give an accurate 
representation in the vicinity of the body and permit 
the pressure distribution on the body to be calculated. 

Such a solution was first obtained by Riabouchinsky, ® 
who avoids the problem of the pressure distribution 
in the wake by constructing a closed wake, or cavity, 
instead of letting it extend to infinity downstream. 
The pressure in the cavity is uniform but may have 
arbitrary values. Thus a family of cavity flows is 
obtained in which the dimensions of the cavity, as well 
as the pressure distribution and drag of the body, de- 
pend on the cavity pressure as parameter. However, 
Riabouchinsky’s solution was not exploited by him 
or anyone else for other than cavitation flows. 

In a somewhat different approach taken by the pres- 
ent writer,‘ the Kirchhoff method is simply modified 
to allow arbitrary base pressures. To demonstrate 
how this is accomplished, Kirchhoff’s example of a flat 
plate set broadside to the (two-dimensional) flow is 
sketched in Fig. 1. The essential step, following 
Helmholtz and Kirchhoff, is the mapping of the flow 
onto the hodograph plane, from where it may be 
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mapped onto the plane of the complex potential, and 
the potential flow completely determined. The vector 
from 0, in the hodograph plane, to any other point rep- 
resents the velocity at the corresponding point in the 
Thus the velocity which at 


physical plane. “a’’ has 


the free stream value Ll’, decreases along ao until it is 
zero at the stagnation point. Then it increases again, 
along the front of the plate, until it reaches the value 
U, at the separation point on the edge of the plate. 
Whereas in the Kirchhoff theory this is fixed at the free- 
stream value, here we allow it to have an arbitrary 
value U,, greater than LU’ Correspondingly, the pres 
sure there, and on the base, will have the value 


pbs = p,. + (1/2) [p(l’ 2— U,*)] 


which is lower than p,,. Written asa pressure coeffi 


cient, the base pressure is 


; P, — Pp ta 
Cys = =1- = | 
(1/2)p , 
. U (F) 


where 


From s, in the hodograph plane, the free streamline 
must return to a, the point at infinity. Now one feels 
that, since the pressure is uniform over the base, it 
should be uniform for some distance downstream in the 
wake—that is, that the pressure and the magnitude of 
the velocity along the free streamline are constant at 
first. The corresponding trace in the hodograph plane 
is a circular arc, and we simply assume that the are 
extends through a full 90° as shown, until the direction 
is parallel to the free stream. From this point, ), it 
proceeds to a without further change of direction. In 
this way we obtain, in the hodograph plane, a ‘‘notched”’ 
circular hodograph that can easily be mapped into the 
complex potential plane, and, in the physical plane, a 
wake that has a definite width, d’. The solution de 
pends on the single parameter k 
pressure. The Kirchhoff solution is the limiting case, 
k = 1, shown in the hodograph by the dotted circle. 
A typical pressure distribution for arbitrary k is sketched 


i.e., on the base 


in Fig. 1. 

It is shown in Fig. 4 that if the value of & (1.e., of 
Cys) is chosen to correspond with the measured base 
pressure, then the calculated pressure distributions along 
the free streamline and on the front of the body agree 
well with the measured distributions, so that the meas- 
ured and calculated values of drag are also in excellent 
agreement. On the other hand, such good agreement is 
not obtained when the pressure on the free streamline 
is assumed to start increasing immediately after sepa- 
ration. Therefore the assumption of a circular, notched 
hodograph appears to be a good one and may be ap- 
plied with confidence to other cylinder shapes. In addi- 
tion, it defines a definite wake width, d’, which will 
presently prove to be most useful. The Riabouchinsky 
model gives equally good agreement between measured 
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and calculated drag, for it also assumes constant pres- 
It also gives a definite 
Thus, for repre- 


sure along the free streamline. 
wake dimension for each value of k. 
senting the early part of the flow, the two solutions are 
entirely equivalent, although they differ slightly in 
numerical detail. 

The dependence of the drag (Cp) and the wake 
width (d’ d) on the parameter k, which fixes the base 
pressure, is shown in Figs. 2 and 3, which have been 
adapted from reference 4. 

The Kirchhoff method is by no means restricted to 
this case of the “bluff plate.’’ In fact, there is ex- 
tensive literature on free streamline theory in which 
applications to many different body shapes may be 
Almost always, however, the Kirchhoff as- 
that the base pressure is p 


found. 
sumption is retained 
and correspondingly the resulting values of drag are 
always too low. For instance, when applied to the 
circular cylinder,® © the theory gives a drag coefficient 
0.5, as compared to the experimental value of about 1.0. 
In the application of the theory to a circular cylinder, 
or any smooth cylinder without edges, some additional 
difficulties are encountered, for in such cases the sep- 
aration point is not fixed once and for all. It is as- 
sumed then that at separation the free streamline has 
the same curvature as the body surface. This gives a 
unique additional condition that makes the problem 
determinate. Although this does not give the ‘‘pres- 
sure bump’ that is actually required to separate the 
boundary layer, it does give a solution that is reason- 
able in the large, and which may be expected to give 
reasonable values of the drag, provided the base pres- 
sure is chosen realistically and is not restricted to the 
Kirchhoff value. In reference 4 the notched hodo- 
graph has been adapted to this purpose. Just as for 
the bluff plate, a solution is obtained in which the drag 
and the wake width depend on the parameter k. They 
have been plotted in Figs. 2 and 3, along with the re- 
sults for the bluff plate. These figures contain, in 
addition, the results for a wedge of vertex angle 90 
While in each case the calculated flow and the drag 
correspond well with experimental results, provided 
k is chosen appropriately, the theory now is not closed, 
for it does not tell how to determine the parameter k. 
This, in fact, is the position we should expect to be in, 
for we cannot expect to obtain a complete solution with- 
out considering the mechanics of the wake. As will 
be shown in the next section, the wake plays an essen- 
tial part in setting the base pressure, especially that 
part of the wake in the first few diameters downstream 


of the cylinder. 
(III) THe CoupLinG REGION 
Although the free shear layers conform well at first 


with the calculated free streamlines, this agreement 
does not hold far downstream. Instead, the well- 
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known periodic formation and shedding of vortices 
occurs, resulting, further downstream, in the periodic 
wake structure known as the Karman vortex street. 
It is this transition from the potential free streamline 
flow to the wake regime that is at present so little under- 
stood, but which must eventually hold the key to a 
complete theory. The crucial nature of this region, 
which might appropriately be called the “coupling 
region’ between wake and potential outer flow, may 
be demonstrated strikingly by some simple experi- 
ments. 

In the first of these a thin partition, or splitter plate, 
placed along the centerline of the wake downstream of 
a cylinder, is found to stop the vortex shedding. The 
partition need not extend all the way downstream, but 
only for the first four or five diameters. It appears that 
if the two shear layers on either side of the wake cannot 
“‘see’’ each other in the region where they tend to roll 
up, then there is no stabilizing mechanism to fix a 
definite, periodically alternating, vortex formation, and 
the shear layers then break down in some other manner, 
independently of each other. What is more remark- 
able, the splitter plate has an extremely strong influence 
on the drag of the cylinder. In some measurements on 
a circular cylinder, reported in reference 7, the drag co- 
efficient was changed from 1.1 to 0.7 by introducing 
a splitter plate that extended for five diameters down- 
stream of the cylinder. The corresponding change of 
base pressure coefficient, C,,, was from —1.0 to —0.5. 
In short, the prevention of vortex formation reduces 
the base suction. 

Some indication of what is happening may be ob- 
tained from Fig. 4, which gives the corresponding pres- 
sure distributions along the wake centerline, with and 
without splitter plate. This not only demonstrates 
the change effected in the base pressure, but also shows 
a pressure ‘‘valley’’ downstream of the base, especially 
pronounced when the splitter plate is absent. The 
minimum pressure occurs about one diameter Cown- 
stream of the base, just in the region where the vortices 
form. Actually, the pressure there fluctuates at the 
shedding frequency, and it is the mean value that is 
given in Fig. 4. But it seems clear from this that the 
dynamics in the region of vortex formation are such as 
to create a low pressure there and that the suction on 
the cylinder base is intimately associated with it. 

This relation between vortex formation and_ base 
pressure is further illustrated in the measurements of 
Fig. 5, which were obtained by using a much shorter 
length for the splitter plate—namely, about one diam- 
eter. This is too short to completely shield the shear 
layers from each other, but it does interfere with the 
vortex formation. For instance, the shedding fre- 
quency, at fixed free-stream velocity, is changed when 
the position of the plate is altered along the centerline. 
In Fig. 5, the shedding frequency, n, is given in the di- 
mensionless form, S = nd/U., which is called the 
Strouhal Number. Since Ll’, is fixed in this experi- 
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ment, S is directly proportional to the frequency. As 
the plate is moved downstream, the shedding frequency 
decreases appreciably, reaching a minimum value 
when the trailing edge of the plate is at four diameters 
downstream. Then, when the plate is moved just 
beyond this critical position, the shedding frequency 
abruptly increases, almost to its undisturbed value, 
corresponding to flow without interference. 

The effect of the short plate in altering the shedding 
frequency seems to be as follows. When it is not too 
far from the base, the vortices form on its downstream 
side; this downstream shift of the formation position is 
accompanied by a decrease of the frequency. On the 
other hand, when the plate is far downstream it can have 
no effect, and the vortices form in their normal way. 
At some intermediate, critical position of the plate, the 
vortices have a choice of forming on either its upstream 
or its downstream side. Since there is no possibility for 
a gradual transition from one regime to the other, there 
must be an abrupt jump at that position. 

This phenomenon is in itself interesting, but the ex- 
periment demonstrates two other important points. 
First, the interference is effective only in the region 
where the vortices form; farther downstream it has no 
important effect on the vortex shedding. Second, the 
interrelation between the vortex formation and the 
base pressure is clearly demonstrated by the measure- 
ments of base pressure, C,,, which are also given in Fig. 
5 and which were made simultaneously with the meas- 
urement of the shedding frequency. As the plate is 
moved downstream, the decrease in shedding fre- 
quency is accompanied by an increase in base pressure, 
which abruptly decreases again at the critical position. 
The increase in base pressure—-i. e., decrease in suction 
—is no doubt partly due to the removal of the mini- 
mum pressure point from the base as the region of vor- 
tex formation is forced to move downstream. In addi- 
tion, the minimum pressure itself may increase, but 
this point was not investigated. 

The measurements given here were for the circular 
cylinder, but the same effects are observed for other 
bluff cylinders. What is clear from these experiments 
is that in every case the flow in the coupling region, in 
the first few diameters downstream of the base, is 
critical for the determination of the base pressure and 
of the overall flow. The important effect of inter- 
ference elements on the vortex formation, and thus on 
the base pressure, indicates that there will have to be 
some understanding of the dynamics of these vortices. 
That problem is made difficult not only by its nonsta- 
tionary nature, but also by the fact that, for the Rey- 
nolds Numbers of interest, the vortex formation takes 
place in turbulent fluid.* 

This introduces the difficulties that are encountered 
in any turbulent flow and makes a theoretical deter- 
mination of the base pressure parameter, k, seem rather 


remote at present. However, there are other possi- 
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TABLE 1 ts 
Cp Ro 7 

Circular cylinder 1.0 ().2] 

90° wedge is 0.18 

Bluff plate be 0.14 


bilities for exploiting the results of the free streamline 
theory. Two of these are explored in the following two 


sections. 
(IV) BLUFFNESS 


One of the features of bluff body problems which 
has not received sufficient attention is the idea of 
“bluffness’’ itself, though the usage of the term im- 
plies something about the associated flow character- 
istics. However, this usage has been mainly qualita- 
tive or intuitive. With the results of the free stream- 
line theory it becomes possible to place these intuitive 
ideas on a somewhat firmer basis and to correlate them 
with some other experimental observations that are 
not so intuitive. The main features associated with 
bluffness may be enumerated as follows. 

(1) Of two bodies having the same frontal area, the 
bluffer one tends to diverge the flow more, to create a 
larger wake, and to experience a higher drag. The re- 
sults of the free streamline theory are quite in accord 
with this definition, as is evident from Figs. 2 and 3. 
In Fig. 2 the three cylinders are arranged in terms of 
circular cylinder, 90 
Their 


increasing drag, in the order: 
wedge, and flat plate (normal to the flow). 
bluffness is, by definition, increasing in that order. 
Then, in Fig. 3, the same ordering of bluffness is ob- 
tained in terms of the wake width. 

(2) It is a matter of observation that bluffer cylinders 
have lower Strouhal Numbers. For instance, for the 
three cylinders that we have singled out for discussion, 
the typical values of Strouhal Number are given in 
Table 1, together with typical values of drag coefficient, 
for comparison. Fage and Johansen? made the per- 
tinent observation that the shedding frequency is re- 
lated not to the cylinder dimension but to the width of 
the wake, being inversely proportional to it. This ex- 
plains why the bluffer cylinders, which have wider 
wakes, have lower shedding frequencies and corre- 
spondingly lower Strouhal Numbers. 

(3) In the comparisons of bluffness made in Figs. 
2 and 3 and it has been implied that the value of k is 
is the same, or nearly so, for the different cylinders. 
This is fairly well supported by the experimental re- 


sults. For different cylinders the base pressures are 
approximately equal, within about 10 per cent, pro- 
vided there is no interference in the wake. What- 


ever variations of base pressure do exist do not show 
any systematic variation with bluffness but appear, 
rather, to be related to Reynolds Number effects. Thus 
the variation between different cylinders is no more than 
what may be obtained for the same cylinder at differ- 
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ent Reynolds Numbers. The remarks here apply to 
a range of Reynolds Numbers of about 500 to 10° for 
the circular cylinder and to even higher values for the 
other two. This is the range in which C»p is ordinarily 
described as being constant, but actually this also varies 
somewhat with R, and, in fact, the variations are asso- 
ciated with those in C,,. Likewise, the Strouhal Num- 
ber may be described as constant only to the same 
degree, its variations also being associated with the vari- 
ations in C,,. Over this Reynolds Number range the 
values of Cp and S are within 10 per cent of those given 
in Table 1 and, for our present purposes, may be as- 
sumed to be constant. 

(4) The above remarks apply to cylinders with un- 
disturbed wakes. If interference elements are placed 
in the wake, as described in the above experiments, 
then significant changes are produced in C,,, and, cor- 
respondingly in Cp and S. The important result is 
that, for a given cylinder, a decrease in drag is accom- 
panied by a decrease in shedding frequency. This is at 
odds with the behavior indicated in Table | for dif- 


ferent cylinders, but it will be seen to be entirely in 


agreement with the observation that shedding fre- 
quency decreases when wake width increases. For, 
as seen from Figs. 2 and 3, a decrease in drag of a given 
cylinder corresponds to an increase of wake width. 
This seems contrary to intuition, but it may be pointed 
out that the increase in wake width is accompanied by 
a decrease in wake velocities sufficient to give a net de- 
crease in wake energy. 

(It should be noted here that the results of Figs. 2 
and 3 for the circular cylinder are not applicable when 
the separation point is on the back side of the cylinder, 
as happens at Reynolds Numbers above the critical 
value of about 10°. In such cases the wake width will 
be smaller than the cylinder dimension, and increasing 
wake width will correspond to increasing drag. In 
those cases where the separation point is on the back, 
the free streamline theory encounters some difficulties, 
which are connected with the fact that the hodograph 
plane becomes double sheeted.) 

(5) The wake structure of different bluff bodies is 
similar. That is, the shape, or bluffness, of the body 
has no characteristic effect on the wake other than to 
determine its geometrical and velocity scales. 

These several observations on the characteristics of 
bluffness may be incorporated into a simple dimen- 
sional analysis. Consider the two parallel vortex 
sheets sketched in Fig. 6, which are intended to be 
idealizations of the two free shear layers created by 
some bluff body. If these vortex layers roll up to form 
vortices, then the frequency with which this occurs 
should depend only on the characteristic velocity, U’,, 
and the dimension, d’, associated with the configuration. 
It should not depend on the body that created it. In 
other words, there should be a Strouhal Number for the 


wake, defined as 


S* = nd’/U, (2) 


which should be universal for the vortex shedding from 

‘different bluff bodies. Of course it may be expected 
that there will be some variation with Reynolds Num- 
ber, just as for the ordinary Strouhal Number. For 
the present, we propose to take this into account only 
roughly, by defining a wake Reynolds Number 


R* = Ud'/y (3) 


To take the Reynolds Number effects into account more 
precisely will complicate matters considerably. For 
instance, the shear layers set up by different bluff bodies 
will not be the same, even at the same R*, for they 
may have had different histories up to the point where 
they start to roll up. Thus they may be thicker or 
thinner, may have become turbulent earlier or later, 
etc., and these features will have an important in- 
fluence on the vortex dynamics. They will therefore 
be extremely important for the vortex dynamics and 
the base pressure. However, our considerations are 
only kinematic, and, for that, the secondary Reynolds 
Number effects, the details of the free streamline shape, 
etc., may be less important. 

To confirm that such a universal wake Strouhal Num- 
ber exists, some measurements were made on the three 
different cylinders, with and without wake interference, 
as follows. The wake Strouhal Number is related to 
the ordinary Strouhal Number by 

nd d’ U > = 
S* = — — = (4) 

l d kd 
To determine S*, two measurements are needed. 5S is 
determined in the usual way from a measurement of 
the shedding frequency. & is determined from a meas- 
urement of the base pressure, using the relation (cf. 


Eq. 1) 
k=V1-—- Cc (5) 


These are the only measurements needed, for the other 
factor, d’/d, then is given by the free streamline results 
in Fig. 3. The relation (5) may be regarded as the 
essential one in the coupling between the potential flow 
and the wake, for it determines, for a given cylinder, 
the wake velocity and dimension in terms of the base 
pressure. The results, plotted in Fig. 7, show that 
there is indeed a universal relation, S* = S*(R*). 
Possibly the scatter is too large to show any functional 
trend, as has been attempted with the solid line, and 
one should simply assume that S* is constant at about 
the value 0.16. The discrepancies at the lower values 
of R* are probably more than experimental and may be 
related to some of the secondary Reynolds Numbers 
referred to above. At those low values, the free shear 
layers may be laminar for a considerable length, before 
transition occurs, and this transition point may be 
quite different for different cylinders, even at the same 




















130 JOURNAL OF THE AERONAUTIC 
ere? 
4 8 12 16 20 24 26 32 36 40 44 46 
0.20 
eee e 
7 
. e eo o 
o1e}—* 7 
+ 
b 
Ms De an Ge & 
0.16 ‘cal SS Ml o+ 
ad ° 
014 4. 
© ® FLAT PLATE 
6 s 90° WEDGE 
o @ CIRCULAR CYLINDER 
0.12 o FLAT PLATE 
WITH INTERFERENCE 
° CIRCULAR CYLINDER 
ne WITH INTERFERENCE 


Wake Strouhal Number. 


FG. 7. 


It will be noted that the parameter S* incorporates 
the observation made by Fage and Johansen—namely, 
that the shedding frequency depends on the wake di- 
mension rather than the cylinder dimension. S* is 
more general in that it includes also the wake velocity, 
U’,. This was not necessary in Fage’s treatment, since 
all the cylinders considered had essentially the same 
base pressure and thus the same LU’. By including 
U,, the parameter S* can be extended to a larger class 
of flows, in particular those for which the base pressure 
changes appreciably, as for instance with wake inter- 
ference. Still another important point of difference 
is that the wake width is not measured, as in Fage’s 
case, but is calculated from the free streamline theory, 
using the measured base pressure. This is especially 
convenient when the results are used in the reverse 
direction. For instance, it is possible to infer the drag 
of a given cylinder simply from a measurement of its 
shedding frequency, for then the factor 


k(d/d’) = S/S* = S/0.16 (6) 


is readily computed and & is determined with the help 
of Fig. 3. The drag coefficient is then obtained from 
Fig. 2. (In case S* varies a little with R*, the calcu- 


lation may still be made easily by iteration.) 


(V) ON CLOSING THE KARMAN SOLUTION 


There is another way in which the potential free 
streamline flow may be joined to the wake. It will be 
recalled that in Karman’s treatment of the 
street there are two undetermined parameters, which 
These two param- 


vortex 


are usually found from experiment. 
eters, which are essentially a dimension and a velocity, 
are just the ones that are needed to relate the vortex 
street, or wake, to whatever body produces it, as ex- 
plained in the last section. Now since the free stream- 
line theory also gives a velocity and dimension (U’, and 
d’), which vary with k, it would appear to be straight- 
forward to determine k uniquely by combining the two 
theories. For instance, the drag must be the same, 
whether calculated from the free streamline theory or 
from Karman’s drag formula for the vortex street. 

The possibility of joining the free streamline flow to 
the vortex street was first proposed by Heisenberg.’ 
He retained the Kirchhoff assumption that the base 
pressure is at the free-stream value, and this forced a 
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redundancy on his model, for the drag is then fixed at 
the Kirchhoff value. With the free streamline theory, 
in which the base pressure is variable, a more realistic 
formulation is possible. 

Unfortunately, the lack of understanding about the 
coupling region again precludes a strictly theoretical 
solution. For instance, it becomes necessary to know 
how much of the vorticity in the shear layers appears 
in individual vortices farther downstream. One might 
be inclined to assume that all of it does, but Fage and 
Johansen’ showed that such is not the case, and they 
estimated that only about one-half of the vorticity 
goes into individual vortices, the rest being canceled 
by mixing with opposite vorticity in the coupling re- 
gion. This introduces an empirical factor into the 
solution, but it will nevertheless be of interest to carry 
it through. The details are as follows. 

Along a shear layer (Fig. 6) the circulation per unit 
length is U’,, and the average velocity is (1/2)U,, so 
that the rate at which circulation is transported down- 
stream is (1/2) U,*. In the vortex street, where the 
vorticity is concentrated into individual vortices, each 
having circulation I’, the rate of transport of circula- 
tion in each row is mI’, where m is the shedding fre- 
quency. But according to the observation by Fage, this 
is only a fraction, «, of the transport along the shear 
layer—that is 


= (7) 


ek? U 7/2 


nv eU 7/2 


The vortices move past a given point with the velocity 
U. — u, where uw is their velocity relative to the free 
stream, while their spacing along a row is /. There- 
fore, the frequency with which they pass the point is 
(U. — u)/l, which must be equal to the shedding fre- 
quency, x. Introducing this into Eq. (7) and rearrang- 
ing it in dimensionless form, gives 


(1 — (u/U,)] [(T/(U.1)] = €(k?/2) (8) 
Then, by introducing one of Karman’s stability param- 
eters, T'/(u/) = 2 V 2, it is possible to eliminate [ from 


Eq. (8), and to obtain 
Uu l ek? 
—- = iz - y' — 

U. 2 V2 


This may be regarded as the expression that establishes 
the relation between velocities in the wake with those 


(9) 


in the free streamline flow. 
Now Karman’s formula for the drag, which he ob- 


tained from a momentum calculation on the vortex 


street, may be written in the form 


(h/d) { [5.65 (u/U,,)|—[2.25 (u/U,)*]} (0) 


Cp 
Taking h/d to the left-hand side and using Eq. (9) to 
replace u/U',, results in the expression 


9 


2 ek? 
+ 0.9635 —- 
V2 


ek” 


V2 


1+ @/1 (11) 


Cp h 1.70 
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On the other hand, the drag may also be obtained from 
the free streamline theory (Fig. 2) and, when combined 
with the wake width d’ ‘d (Fig. 3), gives the parameter 
Cp(d/d’) as a function of k. We now need an expres- 
sion that establishes the relation between the dimension, 


and for this we may simply write 
h = dd’ (12) 


where A is a factor that is unknown and which again 
expresses the lack of knowledge of the region of vortex 
formation. The solution is then finally determined 


from the equation 
Cp(d/h) = (1/d)Cp(d/d’) (13) 


The left-hand side is a function of k and ¢ and comes 
from the theory of the vortex street. The right-hand 
side is a function of k and \ and may be calculated from 
the free streamline theory for each cylinder. The 
values of k which satisfy both sides and which are the 
solutions are at the intersections of the two functions. 
In Fig. 8, the left-hand side has been plotted for several 
values of «. The right-hand side has been plotted for 
the three different cylinders, with A = 1. A most use- 
ful result is that all three fall on a single curve. (The 
branching at the higher values of k may be due to the 
approximate methods that had to be used for calcu- 
lating the wedge and the circular cylinder. An 11', 
wedge also falls on this universal curve.) 

To actually find a solution, it is necessary to know e, 
which can at present be found only empirically. But 
the important result is that this needs to be done for 
only one cylinder, since the others fall on the same 
curve. In other words, it tells that k is the same for 
for all the cylinders, and this, it was noted earlier, checks 
well with experiment. Variations in k due to vari- 
ations in Reynolds Number or due to wake interference 
would then correspond to variations in e. Actually the 
other empirical factor, A, in Eq. (15) should also be 
included in the discussion. But it will not change the 
fact that all the cylinders lie on the same curve; it will 
merely shift the position of the curve. The value 
\ = 1, which is the only one plotted, seems to be as 
good as can be inferred from the experimental infor- 
mation for cylinders without wake interference. The 
vortex spacing / itself is not well defined experimentally, 
tending to increase downstream. 

Once a value of k is determined, the theory will also 
give the shedding frequency. For instance, the uni- 
versal Strouhal Number is 


a nd’ eU,—-ud (1 _ yO)G) 
~ / U, U 1]\k 


If it is again assumed that \ = 1, then the factor d’// 
may be replaced by Karman’s stability ratio, h / = 
0.281, so that 


S* = (0.281 k) [1 — (u/U,)j (14) 


u/U. may be replaced from Eq. (9). 
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Fic. 8. Wake solutions. 


For cylinders without wake interference, the base 
pressure coefficient for a wide range of Reynolds Num- 
bers lies between the values —0.S and —1.0, tending 
more nearly to the latter value at high Reynolds Num- 
bers. A suitable average value for k then is 1.4. On 


Fig. S this gives an intersection at Cp(d/d’) = 0.96 
and ¢€ = 0.43. Then from Eq. (9), w/ U., = 0.18, and, 
from Eq. (14), S* = 0.164, remembering that \ = 1. 


These are all in fair agreement with the experimental 
information. The drag coefficients corresponding to 
this value of & are 1.10, 1.32, and 1.74 for the circular 
cylinder, 90° wedge, and flat plate, respectively. The 
corresponding Strouhal Numbers calculated from Eq. 
(4) are 0.206, 0.167, and 0.127. Both drag coefficients 
and Strouhal Numbers compare favorably with the 


values listed in Table 1. 


(VI) CONCLUDING REMARKS 


The experiments with the interference elements dem- 
onstrate convincingly that the critical region for the 
determination of the whole flow field of a bluff body is 
the region of vortex formation in the first four or five 
diameters behind the body. Some understanding of 
this region will be necessary for a really complete theory 
without empiricism. 

The flow in the vicinity of the body, as weli as the 
drag, is described adequately by the free streamline 
theory, provided the right value is chosen for the base 
pressure. Thus the theory depends on the parameter 
k, but it might be pointed out that this is the only em- 
pirical parameter. In addition, the experimental evi- 
dence indicates, and the considerations of Section V 
confirm, that k is the same for different bluff bodies 


under similar wake conditions. This allows, in prin- 
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ciple the solution of different bluff bodies from only a 
single measurement on one of them. 

The free streamline theory, based on the notched 
hodograph or on Riabouchinsky’s model, defines a 
wake width that varies with base pressure. This is a 
relationship that is extremely useful for correlating the 
real wake with the flow near the body, for the drag of a 
given body is then associated with a wake dimension 
and a wake velocity. It makes possible, for instance, 
a simple similarity treatment of the vortex shedding 
and the definition of a universal Strouhal Number, 
S*, with which the drag may be determined from a 
measurement of the shedding frequency, or vice versa. 
Other wake phenomena might be correlated in the same 
way. 

The results of this paper are applicable to any other 
two-dimensional bluff body shape if the free streamline 
theory is worked out for it. This sometimes presents 
computational difficulties, but often no more than for 
the cases with Kirchhoff’s assumption to which a great 
deal of effort has been devoted. Most of the consider- 
ations of this paper could also be applied to three- 
dimensional bluff bodies, especially bodies of revolution, 
if the appropriate free streamline flows could be worked 
out. Here the methods of conformal mappings are 
not available to facilitate the calculations, nor is there a 
theory of the wake corresponding to Karman’s theory of 
the vortex street. 

In the experiments presented here the maximum cyl- 
inder Reynolds Numbers were about 20,000, but the 
same treatment should apply to much higher Rey- 
nolds Numbers. For bodies like the flat plate and 
wedge, in which the separation point is fixed, this should 
be straightforward. For cases like the circular cylinder, 
in which the separation point moves to the back at 
higher Reynolds Numbers, it is necessary to work out 
streamline theory. Otherwise, the 


a suitable free 
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interaction between wake and cylinder should be 
qualitatively the same as at the lower Reynolds Num- 
bers. For instance, some recent measurements” on 
circular cylinders above the critical Reynolds Number 
show that the Strouhal Number is about 0.4 when the 
drag coefficient is 0.4. The typical value of base pres- 
sure at these Reynolds Numbers is —0.4, for which 
k = 1.2. Assuming that the universal Strouhal Num- 
ber, S* = 0.16, applies here, the wake width may be 
calculated from Eq. (4). It is found to be d’/d = 
0.5, which is reasonable, and supports the view that 
the processes in the wake are essentially unchanged at 
Reynolds Numbers above critical. 
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. British Fredric F. Ehrich* 1 
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ECENT DISCUSSION in this column!~‘ has centered on the prob- 
lem of numerical differentiation of experimental data. = I - ; ype 
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tives at the central point of an odd number of uniformly spaced 
points. These formulations derive from fitting high order poly- 
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where N is an ode ber, asin Fig. 1. / aight line defined by . : — - . 
is an odd number, as in Fig. 1. A straight li ined by since all its derivatives are equal. Then using the table 


of Cyl- y= Cx+C (2) 
53. y ? pas 
may be fitted by least squares to the points giving, after some n Xn Vn 
manipulation and simplification —2 3. 80 44.70) 
— ; 90) 9.40 
Wx) = Cx + Co = (1/N) ) Vy (3a ; ‘ 00 Mr = 
n 1 $10 60.34 
12 2 $ 9) 66.69 


y(x) = Ci = S~ ny (3b 


we may apply the equations for finding the derivatives and get 








——— A parabola defined by 
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P : = Least Square Fourth Order 
may be fitted by least squares to the points giving, after a Exact Parabola Polvnomial 
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It is seen that the fourth order polynomial is completely inade- 
quate for calculating the second derivative. 

This least square technique does not lend itself to calculation of 
derivati.es at points other than the central point. 

It should be noted that Eqs. (3a) and (3b), and Eqs. (4a,) (4b) 
and (4c) are extremely efficient and straightforward formulas for 
least fitting—i.e., calculation of the co- 


the actual squares 


efficients. 
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A Particular Solution of a Nonlinear Differential 
Equation Which Occurs in the Dynamics of 
Conservative Systems and Its Application to 
Undamped Airplane Motion 


Eugene Covert 
Aeronautical Engineer, Naval Supersonic Laboratory, 
Massachusetts Institute of Technology, Cambridge, Mass. 


September 20, 1954 


_* THE OPERATING SPEED and altitude of airplanes are in- 
creased, the constant-speed stability equation is character- 
ized by a reduction in damping and the appearance of a non- 
linear pitching moment term. The loss in 
damping is due primarily to the high speed and altitude while 
the introduction of the nonlinear spring is characteristic of flight 
Since, as the damping becomes small, 


“spring’”’ 


(i.e. 


at transonic speeds. 
the response of the airplane approaches that given by the solu- 
tion of the equation of motion for zero damping,! the solution 
for undamped motion is of some interest to the designer. The 
purpose of this note is to formulate the solution to this problem 
when the forcing function is a step function, or an impulse func- 
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tion. These solutions can be found using the classical methods 
of mathematics; however, the presentation below is easier to 
interpret. 

If @ and 6 are the variations in angle of attack and elevator 
deflection from a trim position, the governing undamped differ- 


ential equation can be written as? 


a+ ma) = m5-d6(t) (1) 
and 
6 = a a = a (2) 
erTrim Vrrim 
where 
Vrrim = forward velocity «t trim 
g = acceleration due to gravity 
0 = pitch angle, positive nose up 


and positive 6 is such as to pull the nose up. The dot denotes 
differentiation with respect to time. The spring constant [i.e 

m(a)| can be found by matching the trim pitching moment curve 
with the proper cubic in a. 
pitching moment curve due to elevator deflection 


restricted to a step function having a magnitude 69[1(t)] 


mg is proportional to the slope of the 
6(t) is here 
or the 
Dirac impulse function [6o(¢| 0)]. 

By changing the variables from a and ¢ to a and & Eq. (1) be- 


comes 
d (a& 7 
+ ma) = m56 Cl 3) 
da 2 
a 
If m(a) = aa* + ba? + ca and nla)da = Mla 
J ald 
then Eq. (3) may be now integrated once to give 
a? a0) | s; . P 
= — - M(a) = ms:6(t) da = mgadt)-d(*) dt 
- - a(O “0 


(4) 


The right-hand integra! in Eq. (4) can be evaluated for the step 
function by integrating by parts once and recognizing that the 
derivative of a step function is a Dirac function [é:(t| 0)|, thus 


t t 
m6 [ al(t) dt = m5-6) | a — aby(t}0) dt | = 
) /7 


oe ( 


c i] - 
ms:6o ja — a(O)} (5) 


When 6(¢) is an impulse the angular velocity has a jump propor- 


tional to 6) at zero time, that is 


ey 
a(t) = &0O) + lim mds 6o(¢)'O) dt (6) 
«— 0 —2 
Substituting this expression in the right-hand side of Eq. (4), 
then 
’ 11569 \ 2 
m56(t)a(t) dt = msivd(O) + 7) 
») 
0 2 


If the results of Eqs. (5) or (7) are denoted by g(a), then Eq. (4) 
can be solved for @ and integrated to give 


{" dn 
+- - =} (8) 
F V 20(n) — 2M(y) + &(0) 


(0) 
where 7 is a dummy integration variable. The sign on the in- 
tegral is chosen to correspond to the prope initial direction of the 
if a(Q) < 
O the positive sign is chosen. 
nomial, the solution can be expressed in terms of inverse Jacobian 


0 choose the negative sign and if a(0) > 
Since ./(n) is a fourth degree poly- 


motion, e.g. 


elliptic functions*® 4 and can be inverted to give solutions in terms 


of a = alot, k). oandk are constants in the solution 


To illustrate the above procedure consider the differential 


equation for a step input. 


@ + 1.25a3 — 1.75a? + 04a = 0.8691(t) (9) 
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with initial conditions 
a(O0) = a(O) = O (10) 


Substitution in Eq. (8) gives 


%a 1 
an 
= } = (11 
0 V 1.6607 — O.4n? + 1.17n* — 0.63n' 


In this case the positive sign is used with the positive 6) and the 
negative sign with the negative 6 Due to the initial conditions 
the motion always takes place between the root » = O of the 
When 6) = 0 
Fig. 1 


radical and its next positive or negative root. 
there is a double root at y = 0 which implies no motion. 
shows the maximum amplitude of the motion and the static trim 
angle as a function of 69. It will be noted that for a relatively 
small positive change in control setting there is a large change in 
trim and overshoot. This is caused by the unstable region in the 
moment curve. The lower right-hand part of the figure shows 
the behavior of the two curves on a larger scale in the vicinity of 
this control deflection. 

For small 6) the solution can be written as* 
a snot, k) 


b—¢ sn? (a, s k) 


Se = 


and for large 6) (see reference 4) 
all — cn(ot, k)| 


b + ccn(ct, k) 


a 


where a, b, c, ¢, and & are constants, depending upon the roots of 
the quartic in the radical. Once a(@) is known, then @(f) can be 
found simply from Eq. (2) by numerical integration. 

The most obvious effects of the nonlinear moment curve upon 
the undamped motion are: (1) the overshoot is greater than 
twice the trim or less than twice the trim depending upon whether 
the negative slope of the moment curve is decreasing or increas- 
ing; (2) near the instability region small changes in control de- 


FORUM 


flection have a large effect upon the airplane orientation and (3 
the frequency of the motion depends upon the size of the input 
Thus a great deal of useful information concerning the nature of 


the motion can be found from the elementary analysis of the 


undamped case 


REFERENCES 
1 Trilling, L., and Covert, E. E., On the Quasi-Steady Pitching Oscillations 
{ircraft To be published 
Dynamics of the Airplane, Division N, Vol. V 
Julius Springer, Berlin, 1935, pp. 99, 100 
3rd Rev. Ed.; Ginn and Co 


of Transont 
Jones, B. M 
Theory (W. F. Durand, Ed 
Peirce, B. O., A Short Table of Integrals 
Boston, 1929, pp. 66-72 
4 Grébner, W., and Hofreiter, N., Jntegrultafel Erster Teil, Unbe 
Vienna and Innsbruck, 1949 


Aerodynamic 


limmle 


Integrule; Springer-Verlag 


Further Comments on ‘‘The Effect of Boundary 
Layer on Sonic Flow Through an Abrupt 
Cross-Sectional Area Change’’! 


Robert S. Wick 
Senior Research Engineer, Jet Propulsion Laboratory 
Institute of Technology, Pasadena, Calif. 


September 20, 1954 


Caelitornia 


yes AND MILLs? have published experimental data on sonic 
flow through axially symmetric sudden enlargements of 
cross section which can be compared directly to the author's! re 
sults obtained for sonic flow through two-dimensional rectangular 
sections. The test section used by the author was of constant 
width, and the sudden enlargement consisted of a sudden change 


Thus, there were two corners 


in the height of the test section 
in which a base pressure could be measured, and, since these cor- 
ners were not interconnected, it is possible to have a nonsym- 
metrical condition if the pressure downstream of the expansion is 


raised high enough relative to the entrance condition P This 


FLOW DIRECTION——> 





A.) SYMMETRIC FLOW IN 
SUDDEN ENLARGEMENT 





B.) ASYMMETRIC FLOW IN 


SUDDEN ENLARGEMENT 


Hydraulic analogy results show possible asymmetry in 
flow through sudden enlargements 


Fic. 1. 
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Base pressures measured by hydraulic analogy during 
asymmetric flow through a sudden enlargement. 
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The effect of reservoir pressure on base pressure in 
Comparison of two-dimensional tests to 


Fic. 3. 
sudden enlargements. 
axially symmetric tests. 


was observed by the author in both the two-dimensional, air- 
flow apparatus and water-table hydraulic analogy tests. Two 
photographs are presented in Fig. 1 of water-table results which 
illustrate the possible asymmetry in the two-dimensional ap- 
paratus. Fig. 2 illustrates qualitatively from water-table data 
how the base pressure ratio varies as the asymmetry increases 
with decreasing reservoir pressure jy (or increasing downstream 
pressure ). 

In the axially symmetric tests of Cole and Mills, the base 
pressure over the annular surface must be constant, since all 
In Fig. 3 some data of Cole 
It must be 


corner regions are interconnected. 
and Mills are compared with those of the author. 
remembered that the base pressures reported by the author were 
measured in only one corner of the two-dimensional flow appara- 
tus. (The nomenclature in the figures is the same as that used 
in the author’s paper.) The explanation for the fact that the 
base pressure ratio continuously decreases as the reservoir pres- 
sure Po is raised until a certain constant value is reached is that 


the the downstream section is subsonic at the lower 


flow in 
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reservoir pressures, and therefore the base pressure is determined 
by conditions downstream of the corner. From a one-dimensional 


analysis [Eq. (5) of reference 1] 


PiAy _ Mz §1 + [(y — 1/2) | 
Plz M, 1 + [(y — 19/2109 
and the fact that P;/P) = 0.528 for sonic flow at the throat as 


well as the assumption that for subsonic flow in the expanded 
section P» is approximately equal to P; for subsonic flow, we ob- 


tain 
Me + [( 1)/21M! (0.528) (1.2) (.41/A2) 
iM of ( _ )/2] M7) 7? = 9 
P./Po 


This one-dimensional result indicates that as the pressure ratio 
P;/Po is lowered the Mach Number .V/2 in the expanded section 
may be expected to increase toward unity and also that the value 
of P;/Po at which this occurs is lower the smaller the expansion 
ratio A,;/A», The latter observation was found experimentally 
to be true by both Cole and Mills and the author, as is shown in 
When the flow in the expanded section becomes super- 
pressure in the corner, is no longer 


Fig. 3. 
sonic, the base pressure, i.e. 
determined by conditions in the downstream section but rather 
by conditions upstream of the sonic throat. In this flow regime 
both Cole’s and Mills’ and the author’s data indicate that the 
base pressure ratio is a function to the first order of the expansion 
ratio. The diameter of the throat of Cole’s and Mills’ apparatus 
was */s-in., whereas the dimensions of the throat used by the 
author were 0.380- by 0.500-in. excellent 
opportunity is afforded to compare the two-dimensional tests 
with the axially symmetric tests in the regime where upstream 
conditions may be expected to govern the base pressure ratio. 
The three sets of points of the author in Fig. 4+ correspond to 
different throat lengths as indicated by the ratio 1,/h,. The 
nozzle of Cole and Mills consisted of a thin nozzle plate with a 
'/,-in. radius on the entrance edge, and consequently the straight 
section of the throat must have been very short (approximately 


Consequently, an 


0.2-in. ). 
The theoretical calculations of Korst* for the two-dimensional 


case are also shown in Fig. 4. The agreement appears to be very 


good for L/h; = 3.94. However, Korst states that the trend of the 


0 1 
O.8F WICK 
L,/H, , 
CF 3946 | KORSF— 
7.88 ® a Pp 
l1gl © ra =e 
°| “COLE __ 
BA “7 & MILLS 


0.2 y 
of /’ 


/ 














| 


0.4 














0. y 
0.08 e 


O O2 























06 O08 10 


A,/A, 


The effect of expansion ratio on base pressure in 
Comparison of two-dimensional tests to 
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data with respect to the effect of upstream boundary layer is incor- 
rect as indicated by his theory. Korst attributes this to three- 
dimensional effects of the side walls, which is entirely possible 
when one considers that the available air supply limited the 
tests to an apparatus with h,/W = 0.7. 

An extension of the interior flow tests to the axially symmetric 
case such as Cole and Mills have done but with the throat sec- 
tion lengthened so that different boundary layers are present 
would be very worth while. This additional information would 
aid in understanding the effect of expansion waves on the wake 
boundary layer (see Figs. 1, 2, and 8 of reference 1). This 
effect is not present in the case of external flow, and one would 
expect the effect to be more pronounced at entrance Mach Num- 
bers in the neighborhood of unity. As stated, the validity of 
comparing external flow to internal flow is limited by a number 
of factors. The chief advantages of studying internal flow are 
that from an experimental point of view the experimental equip- 
ment is cheaper and is easier to instrument for probing and for 
obtaining static pressures and optical data. Also, it is possible 
to have complete information on the nature of the boundary 
layer at the entrance to the sudden expansion. 

As Korst points out, Eqs. (2) and (3) of the author’s paper 
were incorrectly stated. Also, Eq. (4) should read 

Py (1 + yM2?)P2/Pi — (1 + yMi?)Ai/A2 

P, 1 — Aj/A2 
The caption of Fig. 4 of the author’s paper should include the 
information that L. = 7.0, in. and H, = 0.38 in. 

Professor Korst in a private communication indicates that 
he has extended his theory to the three-dimensional case and 
that within the range calculated (A;/A2 = 0.6 to 1.0) it agrees 
with the experimental data of Cole and Mills. 
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Determining the Mode Shapes and Frequencies 
of Low Aspect Ratio Multispar Wings 


Marvin L. Stein 
Senior Research Engineer, Convair-San Diego Division of General 
Dynamics Corporation, San Diego, Calif. 


September 28, 1954 


es DYNAMIC EQUATIONS for determining the vibration modes 
and frequencies of a low-aspect ratio multispar wing assume 


the form 


Pf = w?®Mf (1) 
where P = (P;;)(i, 7 = 1,2,...,) is a matrix of elastic coeffi- 
cients and 

M, 0 
M={| % 
0 MM, 


is a diagonal matrix whose nonzero elements are concentrated 


masses The scalar parameter w? represents the square of the 
frequency while the vectors f which satisfy Eq. (1) for any given 
value of w? give the mode shapes corresponding to that frequency. 
From physical considerations ./ is positive definite while Pi; = 


P;; for all and j——i-e., P is symmetric. 


In general there will be only m < nm distinct values of w* for 


which (1) has a solution f The computational problem is to 


s actually determine the allowable frequencies or eigenvalues and 
their corresponding mode shapes or eigenvectors. An extensive 
mathematical literature deals with this problem (see, for example, 
Hestenes and Karush?) 
An equivalent formulation of Eq. (1) which is easier to handle 
computationally can be obtained by introducing the matrix 


l 
VM, 0 
l 
O VM 
] 
0 - 
VM, 
One can then consider the eigenvalue problem 
is = w* (2) 
where 
A = QPQ, (3) 


for if w? is an eigenvalue of (2) and s the corresponding eigen- 
vector, then w? and 

f Qs (4) 
Since Q and P are both 


can easily be shown to satisfy Eq. (1) 
This fact is essential to the 


symmetric, A will also be symmetric. 
mathematical theory on which the computational technique 
is based and allows for considerable simplification of this 
technique. 

A theory for solving eigenvalue problems of type (2) when A is 
symmetric has been developed in detail by Hestenes and Karush.' 
The ‘‘near-optimum”’ gradient method to be outlined here is a 
modification of the Hestenes-Karush method of ‘“‘optimum-a.”’ 
The application of this modified optimum-a technique to finding 
eigenvectors corresponding to zero eigenvalues has been pre- 
viously described by Stein.‘ The generalization to eigenvectors 
associated with nonzero eigenvalues proceeds as follows: 

(a) Choose an arbitrary initial vector 2; 

(b) Compute the Rayleigh quotient yu(s;), where 

p(s) = (s;*, As;)/(s:*, 2) 

(c) Form the gradient vector £(s,;), where 


2:) & = As; — p(s;)2 


wt 


(d) Compute the Rayleigh quotient for the gradient 
w(é:) = (&:*, Ab) /(E*, &) 
(e) Form the scalar 
a = [w(ti) — ls 
(f{) Choose 8 on the range 0 < 6 < 1 
(g) Compute the new estimate 
3i4+1 = Si — Bat 


and repeat, replacing 3; by = 

The values of the quotient u(s;) generated by the above process 
will converge monotonically downward to the lowest eigenvalue 
will converge to the corresponding 


while the sequence of vectors ¢ 
The monotonicity of » affords a good computa- 
The gradient vectors £ go to zero and give an 
as eigenvec tor and eige n 


eigenvector 
tional check 
indication of the accuracy of 2; and wu 
value extimates 


The same procedure will also serve to find the 
Howse ver, care 


next highest 
eigenvalue and its corresponding eigenvector 
must be taken to choose the initial estimate orthogonal to all 


eigenvectors y; that are already known, i.e., to choose 2 such that 


(z0*, 9x) = 0 i Be. ee 5) 
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The method, in theory, will preserve the relations (5) for all suc- 
ceeding estimates 3; However, in practice, an occasional re- 
orthogonalization will be required because of the round-off error. 

IBM computer techniques for carrying out a gradient method 
known as the method of ‘‘fixed-a’’ have been described by E. 
Yowell. 


ployed in the present method a separate description will not be 


Since these techniques are much like those to be em- 
necessary here. 
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Slender Bodies of Low Wave Drag 

R. N. Haskell and G. C. Grogan 

Convair-Fort Worth Division of General Dynamics Corporation, 
Fort Worth, Tex 

September 23, 1954 


eons OF MINIMUM wave drag slender bodies of revolution 
which have been published to date are all restricted to that 
class of shapes termed ‘‘smooth.’’ This requires the cross- 
sectional area gradient to be continuous everywhere including the 
nose of the body. Such a restriction is often incompatible with 
current design requirements. 

A particular case of importance is that of a ducted body whose 
inlet radius and radius slope are largely set by internal aerody 
namic considerations. It was to solve this problem that the 
authors developed a general method of designing low (though not 
necessarily the absolute minimum) wave drag bodies of revolu 
tion. To illustrate the method the following specific example 
will be studied 

Let the body be of unit length, and the following geometrical 


constraints be imposed: 


R(O) = radius at station x = 0 
R'(O) = radius slope at x = 0 
R(1) = radius ata = 1 

R'(1) = O = radius slope at x = 1 


The problem is to determine the body shape satisfying the 
above constraints and having low drag. 
The appropriate drag expression is (from Lighthill! 


1 fl 
D 1 - = ' 
= — S'(x)S"(t) log |x — t| dt dx — 
qd 2a JI) /0 
1 


[.S’(O)]? 2 


log 
BR(O 


S’(0) ss 
S"(x) log x dx + 
<7 


where 
S(x) = cross-sectional area at station x 
S"(x) = d2S(x)/dx? 
/ 
8 =V M?-1 
By letting the area distribution be a power series as 
n 
S(x*) = = aix'" (2) 
2 0 
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the above drag expression may be integrated giving 


n n n 
D l S(O) i 
= YD ayirs; 4 >» a 
q 7 ij 2% 2 


, ¢ 271-1 
[S’(O)]? 2 
as aro) 
where 
hi, a vs 1) oe l 7 l Saeed l 
Sty SlLtetry — 2 2 3 t—] 


The constraints for the variational problem may be expressed 


as: 
S(O) = ao = rR2(0) 
S’(0) = a, = 2rR(0)R’(0) 
SQ) = SO ay = rR%X1) 
1 0 
1" 
S(1) = > ia; = 0 
1 0 


By introducing 
minimized is: 


1 S10) < i [.S'(O)]? 
de ajasbi; 4 » ai 4 
wy] . - 4 T I 21—] 2r 
9 us 
log — + Ay Zz. a; +a a, — S(1 
BR(O) > 6 
\ Z. a T a (4 


since d and q are already known 

The above is minimized by differentiating it with respect to 
the a,’s (ranging from 7 = 2,3,...,”) and, A». The resulting 
simultaneous equations to be solved for the unknowns are: 


2boo (bos + dao) h 1 he) 
a» + a; 4 : a 
Tr T Tr 
1A; + 2X. = 0 
(bh + Bao) 2h (b 7+ D 
ade a3 T a 
7 Tv us 
1A, 3A 0 
dD ’ (b h ) ‘. & 
a, + a a 
7 Tr 
1X ir. = 0 
h be) bs,» + i) 
as + a 
us T 
2b, } 
T a + 1X + wr = | 
7 . . (5 
la, + la; 4+ + la = §(1) — S‘O) — sO), 
2ae + 3a, +... + na, = S’(0) 


Solutions for the case of 


R(O) = 0.05 
R'(0) = 0.1 
R(1) =0.1 
i} = ] 
have been determined for three values of n(n = 3, 4, and 5) and 


the results are shown in Fig. 1. 

As is to be observed the radius curve quickly damps out with 
increasing ” while the corresponding drags decrease. Compari 
son of this method for smooth bodies with minimum drag shapes 
has shown the convergence to be quite rapid and certainly ade- 


quate for engineering purposes 
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Fic. 1. Effect of degree of area distribution power series on pro- | 
file shape and wave drag of example problem. 

| 

| 


The illustrative example may be varied to include additional 








constraints such as prescribed areas and/or volumes. The drag 





(4) expression (3) is appropriate but Eqs. (4) and (5) will have to 

include the additional constraint terms. If it is desired to in- on 9 

clude additional area gradient discontinuities the additional terms 

for the drag expression may be easily developed from reference 1. 

sulting The remainder of the problem is treated in the identical manner 

: as the example f . : 1 


ect to 
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— Difference Solution of Shock Diffraction 
Problem! 








H. F. Ludloff and M. B. Friedman 
Daniel Guggenheim School of Aeronautics, New York University 
October 29,1954 





(5) Have TIME DEPENDENT, nonlinear equations of motion and the 

appropriate boundary condition have been formulated with 
the objective of obtaining a solution describing the reflection and 
diffraction of a strong shock from an arbitrary corner. The 
method of solution used is of more general interest than this 
particular problem. It is a difference procedure with the follow- 
ing characteristics: 

(1) In spite of the conical character of the problem, the equa 
tions of motion are formulated in x-, y-, t-coordinates, in which 
the equations are hyperbolic. 

) and (2) The equations are written in “divergence form” (see au- 


thors’ IAS preprint 456), so that in the limit they go over auto- 








with matically into the conservation laws. 

pari- i (3) The difference form of the equations supplies automatica}y 
lapes a certain amount of artificial viscosity which is necessary for the , . 

ade- | creation of varying shock intensities 4 4 6 


In the proposed difference form: Fic. 2. 
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(a) The four equations of motion, applied at each interior point 
(j, k) of the field, are: 


mags =| es m;t +o. - 3 (p,m, 7, EB), | ia 
aes = 7 u rt + a o 3 (p,m, 1, E), 


In 2,7 runs through 7 + 1,k;7 —1,k;7,k +1;7,k — 1. The 
four dependent variables p, m, r, E are the density, the two mass- 
flow vector components parallel and oblique to the wall, and the 
total energy. R, 7, U, and V are certain, nonlinear difference 
operators, applied to the subsequent quantities p, m, r, E. The 
superscripts refer to the time cycle. The subscripts refer to the 
location of the net points, which occur in certain combinations 
K, 6. (At/Ax) = (At/Ay) are chosen in accordance with the 
stability criterion. 

(b) The one true boundary condition on the wall, and the 
initial conditions describing the field of the plane shock are ob- 
vious. 

Some results of the IBM 701 machine computations are shown 
below. The development of the compression produced by the 
corner can be seen from the plot of the isopyenics (/ = 62.5 X 
p/po = const.) at three subsequent time cycles ¢ = 3, 19, 
and 87. After the incident shock strikes the corner, a small 
circular compression pulse forms which is subsequently distorted 
and swept upstream in the wake of the advancing shock. At 
the same time, the p-gradients at the boundaries of the com- 
pressed domain steepen and assume the expected shape of the 
reflected shock. The density distribution along the wall, the p- 
value at the stagnation point and the angle of the reflected shock 
at the intersection point with the incident shock agree with the 


interferograms. In the vicinity of the latter point, there should 
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be a completely homogeneous domain. At ¢ = 87, such a region 
just starts to develop. 

Note that the way in which the machine evolves, cycle by 
cycle, the final conical field, guided by the hyperbolic equations 
of motion, corresponds to the actual happenings in nature, except 
for an artificial slowdown of the sequence of events to be ex- 
plained presently. 

The initial field of the plane shock front may be described ex- 
clusively in terms of ideal fluid concepts; the same holds true 
On the other hand, the gradual 
i.e., the creation of 


for the final configuration. 
transition from the initial to the final field 
new entropy values, initially not present, can only take place if 
some kind of viscous action is taken into account. Therefore, 
the ideal fluid equations of motion could not describe events, if 
used in their analytic form; inadvertently, the difference form 
(A) introduces the viscous action necessary to reach the objec- 
The difference quotients appearing in (A) are all of first 
order. Now, suppose each difference quotient is expanded about 
the corresponding derivative in powers of a mesh-width param- 
The expanded equations, then, consist of the ideal fluid 


tive. 


eter. 
equations, augmented by second order and higher derivatives 
But the latter appear and behave like Navier-Stokes terms and 
Barnett terms, etc., in which the magnitude of the viscosity is 
related to the mesh width. Naturally, this artificial viscosity 
is not the correct one. However, that does not influence the 
final result which is again independent of viscosity, but only the 
speed with which the final configuration is reached. Thus, it 
happens that the difference formalism yields more than one could 
hope to obtain from an analytic solution; and at present even 
more than the experiment can provide. (See Figs. 1-3.) 
REFERENCE 
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On the Design of Two-Dimensional Supersonic 
Nozzles by the Method of Characteristics* 


Bernard Chaixt and Peter Henricit 
Institute for Aerodynamics, Federal Institute of Technology, Zurich, 
and National Bureau of Standards, Washington, D.C.t 


December 31, 1954 


(1) INTRODUCTION 


5 a DESIGN of two-dimensional supersonic nozzles for a 
parallel and uniform jet by the method of characteristics is 
subject to the difficulty that the initial conditions for the ve 
locity potential in the transonic region of the nozzle are not 
known and have to be derived from suitable assumptions or 
The two following hypotheses 
(a) the 


from independent calculations. 
can be traced back to the early days of nozzle design: 
computation is started at the throat of the nozzle under the 
assumption of parallel flow in the transonic region (see, e.g., 
Preiswerk,? Puckett®); (b) the computation is started in the diver 
gent part of the nozzle under the assumption that the flow can 
be approximated there by a suitable source flow (see, e.g., Ki 
senko,! Busemann,* and Foelsch,’? who, however, replaces the 
method of characteristics by an analytical method; also the “‘hy 
draulic’’ method of O. Reynolds, mentioned by Courant and 
Friedrichs’ rests on similar ideas). Concerning (a), Lighthill® 
and Sauer!® have shown by computations based upon the power- 
series method of Meyer! that the assumption of uniform flow in 


* The work of the second author was supported by the Office of Naval 
Research, U.S. Navy 

+ The authors are indebted to Professor J 
interest in this investigation and for making possible the experimental work 


Ackeret, Zurich, for his active 
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the transonic region does not correspond to reality and that 
the curves J = const. are in first approximation parabolic ares 
open to the nozzle inlet) with a radius of curvature comparable 
to that of the nozzle walls at the throat. Using these results, 
Harrop et al.!* constructed a set of nozzles starting from the 
arbitrarily chosen) line .Z = 1.2, although earlier Erdmann® had 
started from a line along which the velocity is parallel to the axis 
In the present note we describe a method which has been used 
in the design of a set of new nozzles at the Institute for Aerody 
namics at the Federal Institute of Technology in Zurich. This 
method has the feature in common with the one described in refer 
12 that the initial line of the characteristic net is determined 
However, not the line 1J = const 


ence 
by an analytical method. 
but the .Vachline originating from the throat of the nozzle is 
used as carrier of the initial data. This line, which we call S-line 
(or simply S) in view of its peculiar shape, has been determined 
experimentally by reconstruction of the flow from the geometry 
and the pressure distribution of previous nozzles and analytically 
by considerations in the theory of characteristics and by a power 
series expansion of the velocity potential of the flow. In the 
analytical work the basic assumption of references 1, 10, and 12, 
that the velocity distribution on the centerline is linear, is no 
longer made. The geometrical data of the S-line are expressed 
directly in terms of three geometrical parameters of the nozzle, 
namely (i) width of the throat, (ii) curvature at the throat, and 
(iii) rate of change of curvature at the throat Full details of 
the method including a discussion of the boundary-layer effects 
and comparisons with experimental results will be given in a forth- 
coming paper. Here we shall only announce without proofs 
some of the analytical results 
The following notations will be used (see also Fig. 1): 


= Prandtl-Mever angle of expansion 


a = slope of the streamlines 

7 = adiabatic exponent of the gas considered (= 1.4 for 
air ) 

O = sonic point on the nozzle axis 

& = throat of the nozzle, starting point of the S-line 

A = point of intersection of S-line and nozzle axis 

HW’ = turning point of the S-line 

MJ = point of maximum slope in the converging part of the 
nozzle 

R,r = radius of curvature 


O is also the origin of a (x, y)-coordinate system (positive x-axis 
coinciding with the nozzle axis). Small letters with a ~ (such 


as J, 7) refer to the S-line, capitals to the wall. 


(2) AERODYNAMIC PROPERTIES OF THE S-LINE 
Because S is a characteristic line,» + a = wron S. Between 
T and W the characteristics intersecting S are compression lines, 
Between HW’ and A the inter 
> 0 


thus 7” < 0 between JT and W 
secting characteristics are expansion lines and therefore 7” 
The characteristic intersecting S at WV originates at 0 


On this characteristic, a — a = w) = 0. Therefore at W, 


there. 


ww = aw = (1/2)w7 


This relation is exact, and no simplifying assumptions have been 


made in its derivation. It shows that between 7 and W the S 


line is crossed by compression lines whose total strength amounts 
to(!/,)jwr. Since w7 is not negligible,t the effect of this compres 
sion on the downstream region, if not properly corrected, can be 


considerable.{ In particular it appears that method (b) men 


* Mathematically this means that the nozzle contour in the transonic 


region is approximated by a polynomial of third degree 
t We measured w7 i degrees in one Vf 1.66-nozzle 
t It does not seem impossible that the sudden decrease in M, which was 


observed in reference 10 near the inflection point of several nozzles and 
which cannot be explained completely by boundary-layer effects, is due to 


this compression 
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Fic. 1. 


tioned in the introduction vields good results only for nozzles 
with a very flat inlet, for even if the transonic part of the nozzle 
is obtained correctly by backwards construction starting in the 
diverging part, these compression lines coming from the .S-lin¢ 


will not be recovered 


(3) GEOMETRICAL PROPERTIES OF THE S-LINI 
Since the whole characteristic net depends on the shape of the 
initial line, it is important to know S as exactly as possible. We 
characterize S geometrically by expressing coordinates, slope, 
and radius of curvature at each of the three points 7, HW’, and 41 
in terms of the three variables }'7, Ry, and 


d(1/R 
ky’ = ; 
dx ¥=X7 


This last quantity can be interpreted geometrically as follows 


one assumes for the nozzle contour the equation 


] : , kr’ |. 
}-—- \Vr= (Y¥ — X7)? 4 (Y¥ — Xy 


2R7 6 


If this were exact, one had between the maximum slope in the 
upstream part and k7’ the relation 

l l 

2 ky'Rr?* 


With the aid of the auxiliary quantities 


q=Vyt1Rr yr Oke’ 


these twelve relations (of which some are necessarily trivial 


can be expressed as follows: 


Vryti j 2y -— 15 | 
ir = } 1 + p 
™ 6 sl | 0 { 
jr = V7 
V3 j lWg —&y +15 |] 
jr = - p i + p 
Viat l ( ut) { 
VY34+I1.. j 10 VY 3¢ +(—-3 + 4y syy-45 1] 
‘fr = Vrp 1+ ) 
Vytl | Qi) { 
i 
Vv I (—20—10W/2)qt+(4-V/2)¥—-70 | 
ig = Vrp i+ p 
33 | 180) { 
: 1 (5 — 10N/2)g + (18 — 7W/2)7 + 45 | 
be eee Veet + 2 p 
V2V3 | 270 
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3/5 1) 
cor V-V0? 
“Wu = p ; 3 
va I 
J, 10g + (4 — 6/2) y + (15 — 10/2) | 
| GQ) P 
iy = 
Viryti ‘ ; —5q + 27 
é, 2 Vrp 4 + 0 p 
Ya = 0 
V3 j 10g —-8y +15 | ; 
‘ — y) ( y 
ss ee i a . | 90) Py ve 
2 ° y¥+9 
7, = Vrp y, 4 ie o! 
Vrytl 6 { 


These relations have to be interpreted as leading terms of 
power series expansions in p. In order that the approximations 
be good, the values of the curly brackets must therefore be close 
to 1. 
series expansion of the velocity 


All formulas have been derived from a general power 


potential. In that expansion 
sufficiently many terms have been taken into account, so that 
the coefficients in the above formulas are definite.* It is seen 
that the influence of g (and hence of the upstream geometry of the 
nozzle) is only of the second order. 

From the above data it is possible to interpolate the complete 


S-line with sufficient accuracy either graphically or analytically. 


(4) CONCLUDING REMARKS 


Three nozzles for J = 1.37, 1.66, and 1.27, respectively, of 
which the last was constructed completely on the basis of the 
theory presented in this note, have been tested thoroughly It 
was found that the discrepancies between the theory and the ex 
periment, expressed in terms of w, are less than one quarter of a 
degree, whereas for wr the discrepancy between first order theory 
(curly brackets neglected) aml experiment amounted to as much 


as 2 degrees. 


All characteristic nets have been constructed graphically. 
This seemed completely adequate for the present purpose. If a 
very high accuracy is desired, however, it is possible today to 
compute characteristic nets on fast electronic computers.t The 


objections raised in reference 7 against the method of character 
istics are thus certainly no longer justified. 

Finally it should be pointed out that the technique outlined 
above can be applied no matter what hypothesis has been adopted 
for the shaping of the downstream part of the nozzle. In par 
ticular, it can be used in connection with the method of con 
tinuous wall curvature as proposed by Evvard and Marcus!! and 
in reference 12. In the case of the new Zurich nozzles the walls 
were shaped by an iteration method in such a way as to render 


the decrease of the pressure gradient as uniform as possible 
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I 


my paper should be changed to minus 


of Engineering Sciences, Purdue 


N A RECENT PAPER,” Seide has pointed out that the plus sign 


in front of a term, namely, o, (x, 0) in the first of Eqs. (6) of 
Seide has rightfully 
shown that this sign mistake is responsible for finding asymmetri 
cal buckling shapes in place of symmetrical or antisymmetrical 
shapes that were the results of previous studies for simply sup 
ported sandwich columns. I have since then checked the whol 
computation once more and found a few other typographical 
errors and some additional results. I, therefore, wish to make 
these corrections in the order of the equation numbers of my 
paper 

(i) (+) sign in front of a, (x, 0) in Eq. (6) should be (—) 
(ii) As a result of (i), (+) signs in front of Vo, y (m, 0) and 
V4, 
(iii) Replace Eq. (14) by 


(im, 0) in Eq. (7) and in Eq. (11) should be (— ) 
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Buckling stress of a sandwich column E,/2G < p+ 


(u/v)'/2 (with yield correction ) 


Fic. 8 


(E;/12) (mat;/L)? 4 
(G/2) [ua (im, O) F we (m, O)| Ft 


[E,/(Q1 — wv)| (L/mat;) [v1 Cm, 0) + v(m, 0)] 


(14) 
(iv) Use (sinh 22;’t — 7’ sin 22’t) in place of 
(sinh? z,’t — 7’ sin? s2’t) 
in v7,’ (m, 0) of Eq. (17). 
(v) Use (sinh? 2)’ — 7’ sin? 2’t)~! in place of 
(sinh? 2,’ — j’ sin? 22’t) 
in v2’ (m, 0) of Eq. (17). 


(vi) Replace Eq. (18) by 


Ocr - J 


l 


(E,/12) (mat;/L)? 4+ 


G-u; (m, 0) + E,(1 — pv)! -(L/mat;) v(m, 0) 


G-u2.(m, 0) — EF, (1 — pv) (L/mat;) vem, 0) 


(vii) Omit Section (5a) and (5b) of the paper. 

It is curiously interesting that with this (— ) sign correction the 
buckling formulas, Eqs. (14) and (18), are changed radically. 
Eq. (14) does not contain a radical sign, and Eq. (18) does not 
have terms other than the ones contained in a radical. Conse- 
quently, by taking upper or lower signs we obtain antisymmetri- 
cal or symmetrical buckling of previous writers. Of course, the 
argument of Sections (5a) and (5b) of my paper is valid for sand- 
wich columns having unidentical faces. Even a slightest dif 
ference in the thickness or in the elastic properties will give 
asymmetrical buckling pattern. This is also true for other types 
of edge conditions. Eqs. (18) give buckling stress under pure 
bending moments at the edges of the sandwich column. The 
(—) sign of the radical in this equation has no physical signifi- 
cance, since it gives negative buckling stress. 
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New calculations were carried out on an IBM Card Program 
Calculator based on the data shown on Figures given in reference 
I We found that these curves hardly need modification, for the 
errors were less than a few per cent in the minimum critical loads 
This is probably due to the fact that 
soft. 


obtained previously. 


there were all rather Thus, the resulting 


We have, 


however, calculated new type minimum critical stresses that were 


cores chosen 
curves shown there can be used without alteration 
and have extended the range of 


not included in reference 1 


parameters so that the present curves contain the minimum 


critical stresses of previous computations as a point on each 
We, therefore, give (see Figs. 1-8) a set of curves for 
vs. (mzt/L) for various core and face combinations which 
On these 


curve. 
Ocr min 
we hope will be useful in engineering applications 
curves the numbers next to the marked point denote the value of 
(mat;/L) at which o-¢, min. occur. Values of oc; min. are corrected 
for yielding in accordance with the Karman-Engesser formula 
represented by Fig. 15 of reference 1. 

Seide also uses some continuity conditions |[Eqs. (5) and (6) 
of reference 2] replacing the usual simplifications indicated in 
his paper by Eqs. (7) and (8) of the most of previous authors, 
I have had an occasion before to estimate the effect of these sim- 


plifications and had found that for 
E.t./Ejt; KX 1 
this effect is negligibly small. Here E,and E; are moduli of elas- 
ticity, and ¢, and ¢t; are the thickness of the core and faces, respec- 
More- 
namely, « = 0 at the contacting sur- 


tively. This condition is, of course, essential in practice 
over, this approximation 
faces—is consistent with the spirit of the linear theory 

Seide states also that, at least for the problem of buckling 
under compressive edge loads, the elaborate solution of a boundary 
For the 


value problem is unnecessary. This, of course, is false. 


edge conditions Eq. (2) of his paper—i.e., o: (0, y) = o2(L, y) 
for the core—is incorrect and only represents an approximation 
for sandwich columns having soft cores. 

Goodier® gave an analysis that takes this effect into account in 
As soon as o, (0, y) # oz (L, vy) ¥ 0, we cannot use 
Thus we 


one case. 
the solutions of the type of Eqs. (14) of Seide’s paper. 
need to solve a mixed type boundary value problem to which 
In my paper, I used conditions (2) only as a 
Actually I need not use it, since a com- 


Seide objects. 
simplifying assumption. 
plete solution of the boundary value problem taking into ac- 
count the edge values of o, at x = 0, L were given in reference 13 
of my paper, which is being published presently.* 

Finally, the author is indebted to Mr. J. W. Dunkin for carry- 
ing out the computations on the IBM Card Program Calculator. 
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